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We show tliat the fundamental time reversal invariant (TRI) insulator exists in 4+ 1 dimensions, 
where the effective field theory is described by the 4 + f dimensional Chern-Simons theory and the 
topological properties of the electronic structure is classified by the second Chern number. These 
topological properties are the natural generalizations of the time reversal breaking (TRB) quantum 
ffall insulator in 2 + f dimensions. The TRI quantum spin Hall insulator in 2 + f dimensions and the 
topological insulator in 3 + f dimension can be obtained as descendants from the fundamental TRI 
insulator in 4 + I dimensions through a dimensional reduction procedure. The effective topological 
field theory, and the Z-2. topological classification for the TRI insulators in 2 + f and 3 + f dimensions 
are naturally obtained from this procedure. All physically measurable topological response functions 
of the TRI insulators are completely described by the effective topological field theory. Our effective 
topological field theory predicts a number of novel and measurable phenomena, the most striking 
of which is the topological magneto-electric effect, where an electric field generates a magnetic field 
in the same direction, with an universal constant of proportionality quantized in odd multiples of 
the fine structure constant a. = /he. Finally, we present a general classification of all topological 
insulators in various dimensions, and describe them in terms of a unified topological Chern-Simons 
field theory in phase space. 
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I. INTRODUCTION 

Most states or phases of condensed matter can be de- 
scribed by local order parameters and the associated 
broken symmetries. However, the quantum Hall (QH) 
state^'^''^'* gives the first example of topological states of 
matter which have topological quantum numbers differ- 
ent from ordinary states of matter, and are described in 
the low energy limit by topological field theories. Soon 
after the discovery of the integer QH effect, the quanti- 
zation of Hall conductance in units of /h was shown 
to be a general property of two-dimensional time rever- 
sal breaking (TRB) band insulators^. The integral of the 
curvature of the Berry's phase gauge field defined over the 
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magnetic Brillouin zone (BZ) was shown to be a topologi- 
cal invariant called the first Chern number, which is phys- 
ically measured as the quanta of the Hall conductance. 
In the presence of many-body interactions and disorder, 
the Berry curvature and the first Chern number can be 
defined over the space of twisted boundary conditions^. 
In the long wave length limit, both the integer and the 
fractional QH effect can be described by the topological 
Chern-Simons field theorj^ in 2 -f 1 dimensions. This 
effective topological field theory captures all physically 
measurable topological effects, including the quantiza- 
tion of the Hall conductance, the fractional charge, and 
the statistics of quasi-particles*. 

Insulators in 1 + 1 dimensions can also have unique 
topological effects. Solitons in charge density wave 
insulators can have fractional charge or spin-charge 
separation^. The electric polarization of these insula- 
tors can be expressed in terms of the integral of the 
Berry's phase gauge field in momentum space^^iii. Dur- 
ing an adiabatic pumping cycle, the change of electric 
polarization, or the net charge pumped across the ID in- 
sulator, is given by the integral of the Berry curvature 
over the hybrid space of momentum and the adiabatic 
pumping parameter. This integral is quantized to be a 
topological integer—. Both the charge of the soliton and 
the adiabatic pumping current can be obtained from the 
Goldstone-Wilczek formula^'^. 

In this paper we shall show that the topological effects 
in the 1-1-1 dimensional insulator can be obtained from 
the QH effect of the 2-1-1 dimensional TRB insulator by 
a procedure called dimensional reduction. In this pro- 
cedure one of the momenta is replaced by an adiabatic 
parameter, or field, and the Goldstone-Wilczek formula, 
and thus, all topological effects of the 1-1-1 dimensional 
insulators, can be derived from the 2-1-1 dimensional 
QH effect. The procedure of dimensional reduction can 
be generalized to the higher dimensional TRI insulators 
and beyond, which is the key result of this paper. 

In recent years, the QH effect of the 2-1-1 dimensional 
TRB insulators has been generalized to TRI insulators 
in various dimensions. The first example of a topologi- 
cally non-trivial TRI state in condensed matter context 
was the 4D generalization of the QH effect (4DQH) pro- 
posed in Ref. [IJ. The effective theory of this model 
is given by the Chern-Simons topological field theory in 
4-1-1 dimension&i^. The quantum spin Hall (QSH) ef- 
fect has been proposed in 2 + 1 dimensional TRI quan- 
tum modelsi^iii. The QSH insulator state has a gap 
for all bulk excitations, but has topologically protected 
gapless edge states, where opposite spin states counter- 
propagat o^^'^^1^^ . Recently the QSH state has been 
theoretically predicted^^ and experimentally observed in 
HgTe quantum wells^. TRI topological insulators have 
also been classified in 3 -t- 1 dimensions^^'^'^''^'*-. These 3D 
states all carry spin Hall current in the insulating stat^^ . 

The topological properties of the 4-1-1 dimensional TRI 
insulator can be described by the second Chern number 
defined over four dimensional momentum space. On the 



other hand, TRI insulators in 2 -I- 1 and 3-1-1 dimensions 
are described by a Z2 topological invariant defined over 
momentum spac o^^'^^i^^'^'^i^^'^''i^^i^^''^° . In the presence 
of interactions and disorder, the momentum space Z2 
invariant is not well defined, however, one can define a 
more general Z2 topological invariant in terms of spin- 
charge separation associated with a tt fiux'^^'^^. One open 
question in this field concerns the relationship between 
the classification of the 4+1 dimensional TRI insulator 
by the second Chern number and the classification of the 
3 4-1 and 2-1-1 dimensional TRI insulators by the 
number. 

The effective theory of the 4-1-1 dimensional TRI in- 
sulator is given by the topological Chern-Simons field 
theoryi^i^. While the 2 -|- 1 dimensional Chern-Simons 
theory describes a linear topological response to an ex- 
ternal J7(l) gauge fieldii^, the 4-1-1 dimensional Chern- 
Simons theory describes a nonlinear topological response 
to an external J7(l) gauge field. The key outstanding the- 
oretical problem in this field is the search for the topolog- 
ical field theory describing the TRI insulators in 2-1-1 and 
3-1-1 dimensions, from which all measurable topological 
effects can be derived. 

In this paper, we solve this outstanding problem by 
constructing topological field theories for the 2-1-1 and 
3 -f 1 dimensional TRI insulators using the procedure 
of dimensional reduction. We show that the 4-1-1 di- 
mensional topological insulator is the fundamental state 
from which all lower dimensional TRI insulators can be 
derived. This procedure is analogous to the dimensional 
reduction from the 2-1-1 dimensional TRB topological 
insulator to the 1 -f 1 dimensional insulators. There is 
a deep reason why the fundamental TRB topological 
insulator exists in 2 -I- 1 dimensions, while the funda- 
mental TRI topological insulator exists in 4 -I- 1 dimen- 
sions. The reason goes back to the Wigner-von Neumann 
classification^ of level crossings in TRB unitary quan- 
tum systems and the TRI symplectic quantum systems. 
Generically three parameters need to be tuned to ob- 
tain a level crossing in a TRB unitary system, while five 
parameters need to be tuned to obtain a level crossing 
in a TRI symplectic system. These level crossing singu- 
larities give rise to the non-trivial topological curvatures 
on the 2D and 4D parameter surfaces which enclose the 
singularities. Fundamental topological insulators are ob- 
tained in space dimensions where all these parameters 
are momentum variables. Once the fundamental TRI 
topological insulator is identified in 4-1- 1 dimensions, the 
lower dimensional versions of TRI topological insulators 
can be easily obtained by dimensional reduction. In this 
procedure, one or two momentum variables of the 4-1-1 di- 
mensional topological insulator are replaced by adiabatic 
parameters or fields, and the 4-1-1 dimensional Chern- 
Simons topological field theory is reduced to topologi- 
cal field theories involving both the external U[X) gauge 
field and the adiabatic fields. For the 3-1-1 TRI insula- 
tors, the topological field theory is given by that of the 
"axion Lagrangian" , or the 3 -t- 1 dimensional d vacuum 
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term, familiar in the context of quantum chromodynam- 
ics (QCD), where the adiabatic field plays the role of the 
axion field or the 9 angle. From these topological field 
theories, all physically measurable topological effects of 
the 3 + 1 and the 2 + 1 dimensional TRI insulators can be 
derived. We predict a number of novel topological effects 
in this paper, the most striking of which is the topo- 
logical magneto-electric (TME) effect, where an electric 
field induces a magnetic field in the same direction, with 
a universal constant of proportionality quantized in odd 
multiples of the fine structure constant a = /he. We 
also present an experimental proposal to measure this 
novel effect in terms of Faraday rotation. Our dimen- 
sional reduction procedure also naturally produces the 
Z2 classification of the 3 + 1 and the 2 + 1 dimensional 
TRI topological insulators in terms of the integer sec- 
ond Chern class of the 4+1 dimensional TRI topological 
insulators. 

The remaining parts of the paper are organized as fol- 
lows. In Sec. II we review the physical consequences of 
the first Chern number, namely the (2 + l)-d QH effect 
and (1 + l)-d fractional charge and topological pumping 
effects. We begin with the (2+l)-d time reversal breaking 
insulators and study the topological transport properties. 
We then present a dimensional reduction procedure that 
allows us to consider related topological phenomena in 
(l + l)-d and (0 + l)-d. Subsequently, we define a Z2 clas- 
sification of these lower dimensional descendants which 
relies on the presence of a discrete particle-hole symme- 
try. This will serve as a review and a warm-up exercise 
for the more complicated phenomena we consider in the 
later sections. In Sees. Ill, IV, and V we discuss conse- 
quences of a non-trivial second Chern number beginning 
with a parent (4+l)-d topological insulator in Sec. III. In 
Sees. IV and V we continue studying the consequences of 
the second Chern number but in the physically realistic 
(3 + l)-d and (2 + l)-d models which are the descendants 
of the initial (4+l)-d system. We present effective actions 
describing all of the physical systems and their responses 
to applied electromagnetic fields. This provides the first 
effective field theory for the TRI topological insulators in 
(3 + l)-d and (2 + l)-d. For these two descendants of the 
(4+ l)-d theory, we show that the Z2 classification of the 
decedents are obtained from the 2nd Chern number clas- 
sification of the parent TRI insulator. Finally, in Sec. 
VI we unify all of the results into families of topologi- 
cal effective actions defined in a phase space formalism. 
From this we construct a family tree of all topological 
insulators, some of which are only defined in higher di- 
mensions, and with topological Z2 classifications which 
repeat every 8 dimensions. 

This paper contains many new results on topological 
insulators, but it can also be read by advanced students 
as a pedagogical and self-contained introduction of topol- 
ogy applied to condensed matter physics. Physical mod- 
els are presented in the familiar tight-binding forms, and 
all topological results can be derived by exact and explicit 
calculations, using techniques such as response theory al- 



ready familiar in condensed matter physics. During the 
course of reading this paper, we suggest the readers to 
consult Appendix A which covers all of our conventions. 



II. TRB TOPOLOGICAL INSULATORS IN 2 + 1 
DIMENSIONS AND ITS DIMENSIONAL 
REDUCTION 

In this section, we review the physics of the TRB 
topological insulators in 2 + 1 dimensions. We shall use 
the example of a translationally invariant tight-binding 
modeP-S- which realizes the QH effect without Landau lev- 
els. We discuss the procedure of dimensional reduction, 
from which all topological effects of the 1 + 1 dimensional 
insulators can be obtained. This section serves as a sim- 
ple pedagogical example for the more complex case of the 
TRI insulators presented in Sec. Ill and IV. 



A. The first Chern number and topological 
response function in (2 + l)-d 

In general, the tight-binding Hamiltonian of a (2 + l)-d 
band insulator can be expressed as 



H = 



hap f, . 



(1) 



with TO, n the lattice sites and a, (3 — 1, 2, ..N the band 
indices for a TV-band system. With translation symmetry 
^rnn = i^m — ^n), the Hamiltonian can be diagonal- 
ized in a Bloch wavefunction basis: 



H 



^h''^ (k) 



Ck/3 



(2) 



The minimal coupling to an external electro-magnetic 
field is given by h'^n^'^'"" where Amn is a gauge 

potential defined on a lattice link with sites m, n at the 
end. To linear order, the Hamiltonian coupled to the 
electro-magnetic field is obtained as 



H 



^4Mk)ck+5: 

k k,q 



I dhjk) 



with the band indices omitted. The DC response of the 
system to external field A*(q) can be obtained by the 
standard Kubo formula: 



lim —QiAij + iS) 



Qij i^^m ) 



k.n 

•J,(k)G(k,ic^„)) 



^tr(J,(k)G(k,i(tj„ + i^,n)) 



(3) 



with the DC current Ji(k) = dh{\i)/dki, i,j = x,y, 
Green's function G'(k, «tt;„) = [iw„ — ft.(k)]~ , and fl the 
area of the system. When the system is a band insulator 
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with M fully-occupied bands, the longitudinal conduc- 
tance vanishes, i.e. a^x = 0, as expected, while a^y has 
the form shown in Ref. 



e2 1 /■ I 
(^xy — ~J^'i2^ j / dkyfxy (k) 

.,, , ,j ^ day{k.) dax{k) 
with/.,(k) = 

aj(k) = X! ^"^1 ^ I"'*) ' * = 



(4) 



occ 



x,y. 



Physically, ai(k) is the U{1) component of the Berry's 
phase gauge field (adiabatic connection) in momentum 
space. The quantization of the first Chern number 



Ci = J dkx J dkyfxy (k) e ^ 



(5) 



is satisfied for any continuous states |ak) defined on the 
BZ. 

Due to charge conservation, the QH response ji = 
an^-^-' Ej also induces another response equation: 

J, = ane'^E, (6) 
dp „ . _ ^ dB 

ITT = -V ■ J = -(ThV X E = CTff^ 

ot ot 
=^ p(S)-po = '7hS (7) 

where po = = 0) is the charge density in the ground 
state. Equations ^ and ([7]) can be combined together 
in a covariant way: 



ZTT 



(8) 



where /x, z/, r — 0,1,2 are temporal and spatial indices. 
Here and below we will take the units e = = 1 so that 
eV/i= l/27r. 

The response equations ([8]) can be described by the 
topological Chern-Simons field theory of the external 
field A^: 

= i^/^'^ J dtA^e^'^^d.Ar, (9) 

in the sense that SScs/6Af^ — jf^ recovers the response 
equations ([8]). Such an effective action is topologically 
invariant, in agreement with the topological nature of 
the first Chern number. All topological responses of the 
QH state are contained in the Chern-Simons theory^. 



B. Example: two band models 

To make the physical picture clearer, the simplest case 
of a two band model can be studied as an example^^. 
The Hamiltonian of a two-band model can be generally 
written as 



/i(k) =^d,(k)a'^ + e(k)I 



(10) 



where I is the 2x2 identity matrix and cr" are the three 
Pauli matrices. Here we assume that the represent a 
spin or pseudo-spin degree of freedom. If it is a real spin 
then the are thus odd under time reversal. If If the 
(ia(k) are odd in k then the Hamiltonian is time-reversal 
invariant. However, if any of the da contain a constant 
term then the model has explicit time-reversal symmetry 
breaking. If the cr" are a pseudo-spin then one has to 
be more careful. Since, in this case, — 1 then only 
is odd under time-reversal (because it is imaginary) 
while (7^,(7^ are even. The identity matrix is even un- 
der time-reversal and e(k) must be even in k to preserve 
time-reversal. The energy spectrum is easily obtained: 
E±ik) - 6(k) ± V^JW)- When Earfa(k) > for 
all k in the BZ, the two bands never touch each other. 
If we also require that maxk(i?_(k)) < mink(£'+(k)), 
so that the gap is not closed indirectly, then a gap al- 
ways exists between the two bands of the system. In 
the single particle Hamiltonian /i(k), the vector d(k) 
acts as a "Zeeman field" applied to a "pseudospin" ai 
of a two level system. The occupied band satisfies 
(d(k)-cr) | — ,k) = — |d(k)| | — ,k), which thus corresponds 
to the spinor with spin polarization in the — d(k) direc- 
tion. Thus the Berry's phase gained by | — ,k) during 
an adiabatic evolution along some path C in k-space is 
equal to the Berry's phase a spin-1/2 particle gains dur- 
ing the adiabatic rotation of the magnetic field along the 
path d(C). This is known to be half of the solid angle 
subtended by d(C), as shown in Fig[T] Consequently, 
the first Chern number Ci is determined by the winding 
number of d(k) around the origin^^i^^: 



1 

An 



From the response equations we know that a non-zero 
Ci implies a quantized Hall response. The Hall effect 
can only occur in a system with time-reversal symmetry 
breaking so if Ci 7^ then time-reversal symmetry is bro- 
ken. Historically, the first example of such a two-band 
model with a non-zero Chern number was a honeycomb 
lattice model with imaginary next-nearest-neighbor hop- 
ping proposed by Haldane^. 

To be concrete, we shall study a particular two band 
model introduced in Ref.'^^, which is given by 

/i(k) = {sin kx) cr X + {siii-ky) a y 

+ {m + cos kx + cos ky)az, (12) 

This Hamiltonian corresponds to the form with 
e(k) = and d(k) = (sin fca,, sin fcj,, m + cos fca; + cos fey). 
The Chern number of this system is'^^ 



1, < TO< 2 
Ci = { -1, -2 < TO < 
0, otherwise. 



(13) 



In the continuum limit, this model reduces to the 2 -|- 1 
dimensional massive Dirac Hamiltonian 



0=1 



h(k) = kxCJx+kyay + {m+2)az 



m + 2 fcj 



V 

m - 2 



BZ of k 



S2 of d(k) 



FIG. 1: Illustration of the Berry's phase curvature in a two- 
band model. The Berry's phase A ■ dr around a path C in 
the BZ is half of the solid angle subtended by the image path 
d{C) on the sphere 82- 
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In a real space, this model can be expressed in tight- 
binding form as 



H = 



-|-m^cJjtTzC„ 



'-Cn+y + h.C 



(14) 



Physically, such a model describes the quantum anoma- 
lous Hall effect realized with both strong spin-orbit cou- 
pling (cTj, and ay terms) and ferromagnetic polarization 
{az term). Initially this model was introduced for its sim- 
plicity in Rcf. 35, however, recently, it was shown that it 
can be physically realized in IIgi_xMnxTe/Cdi_xMnxTe 
quantum wells with a proper amount of Mn spin 
polarization^^. 



C. Dimensional reduction 

To see how topological effects of 1 + 1 dimensional insu- 
lators can be derived from the first Chcrn number and the 
QH effect through the procedure of dimensional reduc- 
tion, we start by studying the QH system on a cylinder. 
An essential consequence of the nontrivial topology in the 
QH system is the existence of chiral edge states. For the 
simplest case with the first Chern number Ci = 1, there 
is one branch of chiral fermions on each boundary. These 
edge states can be solved for explicitly by diagonalizing 
the Hamiltonian in a cylindrical geometry. That is, 
with periodic boundary conditions in the y-direction and 
open boundary conditions in the x-direction, as shown in 
FiglH (a). Note that with this choice ky is still a good 
quantum number. By defining the partial Fourier trans- 
formation 



Ckyaix) = 



1 



FIG. 2: (a) Illustration of a square lattice with cylindrical 
geometry and the chiral edge states on the boundary. The 
definition of x and y axis are also shown by black arrows, 
(b) One-d energy spectrum of the model in Eq. (|12p with 
m = —1.5. The red and black line stands for the left and 
right moving edge states, respectively, while the blue lines are 
bulk energy levels, (c) Illustration of the edge states evolution 
for ky — —f 2tt. The arrow shows the motion of end states 
in the space of center-of-mass position versus energy, (d) 
Polarization of the one-d system versus ky. (See text) 



with {x, y) the coordinates of square lattice sites, the 
Hamiltonian can be rewritten as 



H 



E 



Cky (x + 1) + h.C. 



+ ^2 '^ky i^) [si^ K'^y + (™ + cos ky)(Tz] Cky {x) 
Y^HMky). (15) 



In this way, the 2D system can be treated as Ly inde- 
pendent ID tight-binding chains, where Ly is the period 
of the lattice in the y-direction. The eigenvalues of the 
ID Hamiltonian HiD{ky) can be obtained numerically for 
each ky, as shown in Fig. [2] (b). An important property 
of the spectrum is the presence of edge states, which lie 
in the bulk energy gap, and are spatially localized at the 
two boundaries: a; = 0, L^. The chiral nature of the edge 
states can be seen from their energy spectrum. From Fig. 
[2] (b) we can see that the velocity v — dE/dk is always 
positive for the left edge state and negative for the right 
one. The QH effect can be easily understood in this edge 
state picture by Laughlin's gauge argument^. Consider 
a constant electric field Ey in the y-direction, which can 
be chosen as 



-Eyt, = 0. 
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The Hamiltonian is written H = ^j^ Hir,{ky + Ay) and 
the current along the x-direction is given by 

J. = 5I'^-(^y) (16) 

ky 

with Jx{ky) the current of the ID system. In this 
way, the Hall response of the 2D system is determined 
by the current response of the parameterized ID sys- 
tems HY£i{q{t)) to the temporal change of the parame- 
ter q{t) = ky + Ay(t). The gauge vector Ay corresponds 
to a flux $ = AyLy threading the cylinder. During a 
time period < t < 2tt /LyEy, the flux changes from 
to 27r. The charge that flows through the system during 
this time is given by 

AQ = dty^Mky) 
Jo I. 

^^AP.(A:,)|f (17) 

with At — 2TT/LyEy. In the second equality we use 
the relation between the current and charge polarization 
Px{ky) of the ID systems Jx{ky) — dPx{ky)/dt. In the 
adiabatic limit, the ID system stays in the ground state 
of Hiuiqit)), so that the change of polarization APx{ky) 
is given by AP^{ky) = Px{ky - ^-K/Ly) - Pxiky). Thus 
in the Lj^ — > oo limit AQ can be written as 

AQ = -f\ky^^^. (18) 
Jo (^'^y 

Therefore, the charge flow due to the Hall current gener- 
ated by the flux through the cylinder equals the charge 
flow through the 1-dimensional system i7iD(fcj,), when 
ky is cycled adiabatically from to 27r. From the QH re- 
sponse we know AQ = anAtEyLy — 2t:(Jh is quantized 
as an integer, which is easy to understand in the ID pic- 
ture. During the adiabatic change of ky from to 27r, 
the energy and position of the edge states will change, 
as shown in Fig[2] (c). Since the edge state energy is al- 
ways increasing(decreasing) with ky for a state on the left 
(right) boundary, the charge is always "pumped" to the 
left for the half-filled system, which leads to AQ = — 1 for 
each cycle. This quantization can also be explicitly shown 
by calculating the polarization Px{ky), as shown in FigH] 
(d), where the jump of Px by one leads to AQ — —1. 
In summary, we have shown that the QH effect in the 
tight-binding model of Eq. (|12p can be mapped to an 
adiabatic pumping effecti^ by diagonalizing the system 
in one direction and mapping the momentum fc to a pa- 
rameter. 

Such a dimensional reduction procedure is not re- 
stricted to specific models, and can be generalized to any 
2D insulators. For any insulator with Hamiltonian ([2]), 
we can define the corresponding ID systems 

H,^{e)^Y.''l.eKkx,e)cu^e (19) 



in which 9 replaces the j/-direction momentum ky and 
effectively takes the place of q{t). When 9 is time- 
dependent, the current response can be obtained by a 
similar Kubo formula to Eq. ([3]), except that the sum- 
mation over all {kx, ky) is replaced by that over only k^. 
More explicitly, such a linear response is defined as 

d9 

Ue) = G{9)- (20) 
0(9) = lim -Q(iu + iS;9) 

Q{iuJn]9) = - ^ iv[Jx{kx]9)GiT){kx,l{Vm+l^n)]d) 

dh{kx;9) ^ . \ 1 
.^^GiD(fc.,*-„;0)j 

Similar to Eq. @ of the 2D case, the response coefhcient 
G{k) can be expressed in terms of a Berry's phase gauge 
field as 

/dh 
-^fxe{kx,0) (21) 

^ J ^ \ d9 ~ Wx) 

with the sum rule 

j G{9)d9 = CieZ. (22) 

If we choose a proper gauge so that ag is always single- 
valued, the expression of G{9) can be further simplified 
to 

0„.|(/^„.,..,»,).^?f). ,.3, 

Physically, the loop integral 

P{9) = <j> dkxax/2TT (24) 

is nothing but the charge polarization of the ID 
systemlSii^, and the response equation ([20]) simply be- 
comes Jx = dP/dt. Since the polarization P is defined 
as the shift of the electron center-of-mass position away 
from the lattice sites, it is only well-defined modulo 1. 
Consequently, the change AP = P{9 = 27r) - P{9 ^ 0) 
through a period of adiabatic evolution is an integer equal 
to — Ci, and corresponds to the charge pumped through 
the system. Such a relation between quantized pumping 
and the first Chern number was shown by Thouless^^. 

Similar to the QH case, the current response can lead 
to a charge density response, which can be determined 
by the charge conservation condition. When the param- 
eter 9 has a smooth spatial dependence 9 = 9{x,t), the 
response equation ([^0]) still holds. From the continuity 
equation we obtain 



dp 


dJx _ 


dt 


dx 




dP{9) 


= 


dx 



d^p{9) 

dxdt 

(25) 
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in which p is defined with respect to the background 
charge. Similar to Eq. ([5]), the density and current re- 
sponse can be written together as 



d P{0{x,t)) 
dxi, 



(26) 



where ^^v ~ 0, 1 are time and space. It should be 
noted that only differentiation with respect to x,t ap- 
pears in Eq. (|26p . This means, as expected, the current 
and density response of the system do not depend on 
the parametrization. In general, when the Hamiltonian 
has smooth space and time dependence, the single par- 
ticle Hamiltonian h{k) becomes h{k,x,t) = h{k,9{x,t)), 
which has the eigenstates \a; k,x,t) with a the band in- 
dex. Then relabelling t,x,k as qa, A — 0,1,2 we can 
define the phase space Berry's phase gauge field 



Aa 



''^{a]qA\ 



d 



dqA 

^AB — QaAb — ObAa 



\a;<lA) 



(27) 




FIG. 3: Illustration of the d(A;, 9) vector for the ID Dirac 
model l|30p . The horizontal blue circle shows the orbit of d(fc) 
vector in the 3D space for k £ [0, 2n) with 6 fixed. The red 
circle shows the track of the blue circle under the variation 
of 9. The cone shows the solid angle fl{9) surrounded by the 
d(fc) curve, which is in times the polarization P{9). 



and the phase space current 

3a = -^eABC^BC- (28) 

The physical current is obtained by integration over the 
wavevector manifold: 

J,- j dkf^=- j^e^^'^T,. (29) 

where fi^v — 0, 1. This recovers Eq. (|26p . Note that we 
could have also looked at the component jk = J dkj^ 
but this current does not have a physical interpretation. 

Before moving to the next topic, we would like to apply 
this formalism to the case of the Dirac model, which re- 
produces the well-known result of fractional charge in the 
Su-Schrieffer-Heeger (SSH) model^, or equivalently the 
Jackiw-Rebbi model^. To see this, consider the follow- 
ing slightly different version of the tight-binding model 

h{k,6) — sin kax + {cos k — 1) (Tz 

+m (sin 9ay -\- cos Ouz ) (30) 

with m > 0. In the limit to <C 1, the Hamiltonian has the 
continuum limit h{k,9) ~ kax + m {sinOay + cos9<Tz), 
which is the continuum Dirac model in (1 -I- l)-d, with 
a real mass to cos ^ and an imaginary mass to sin ^. As 
discussed in Sec. Ill Bl the polarization § dk^a.^ /2t: is 
determined by the solid angle subtended by the curve 
d(fe) = (sin fc, msin^, TO cos + cos fc — 1), as shown in 
Fig. [3l In the limit to ^ 1 one can show that the solid 
angle ^{9) = 29 so that P{0) ~ 9/27:, in which case Eq. 
(pS)) reproduces the Goldstone-Wilczek formulapi^ : 

Jm = -^^ivdv9. (31) 



Specifically, a charge Q = — J'^^{d9/dx){dx/2Tr) — 
~{9{+oo) — 6'(— oo))/27r is carried by a domain wall of 
the 9 field. In particular, for an anti-phase domain 
wall, 9{+oo) — 9{—oo) = tt, we obtain fractional charge 
q = 1/2. Our phase space formula (^5)1 is a new result, 
and it provides a generalization of the Goldstone-Wilczek 
formula to the most general one-dimensional insulator. 



D. Z2 classification of particle- hole symmetric 
insulators in (1 + l)-d 

In the last subsection, we have shown how the first 
Chern number of a Berry's phase gauge field appears 
in an adiabatic pumping effect and the domain wall 
charge of one-dimensional insulators. In these cases, 
an adiabatic spatial or temporal variation of the single- 
particle Hamiltonian, through its parametric dependence 
on 9{x,t), is required to define the Chern number. In 
other words, the first Chern number is defined for a pa- 
rameterized family of Hamiltonians h(k, x, t), rather than 
for a single ID Hamiltonian h{k). In this subsection, we 
will show a different application of the first Chern num- 
ber, in which a Z2 topological classification is obtained 
for particle-hole symmetric insulators in ID. Such a re- 
lation between Chern number and Z2 topology can be 
easily generalized to the more interesting case of second 
Chern number, where a similar Z2 characterization is ob- 
tained for TRI insulators, as will be shown in Sec. IIV CI 
andlVBl 

For a one-dimensional tight-binding Hamiltonian H = 
J2mn'^ma^'^nW'^n/3, the particlc-hole transformation is 
defined by Cma ^"^'^Inp^ where the charge conjugation 
matrix C satisfies C^C = I and C*C = 1. Under periodic 
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boundary conditions the symmetry requirement is 

fe fe 
^ C^h(~k)C ^ ~h^{k). (32) 



From Eq. (|32p it is straightforward to see the symme- 
try of the energy spectrum: if E is an eigenvalue of h{0), 
so is —E. Consequently, if the dimension of h{k) is odd, 
there must be at least one zero mode with E = 0. Since 
the chemical potential is constrained to vanish by the 
traceless condition of h, such a particle-hole symmetric 
system cannot be gapped unless the dimension of h{k) is 
even. Since we are only interested in the classification of 
insulators, we will focus on the case with 27V bands per 
lattice site. 

Now consider two particle-hole symmetric insulators 
with Hamiltonians /ii(fc) and /i2(fc), respectively. In gen- 
eral, a continuous interpolation h(k,9), 9 G [0,7r] be- 
tween them can be defined so that 



h{k,0) = hl{k), h{k,TT) = h2{k) 



(33) 



Moreover, it is always possible to find a proper 
parametrization so that ft,(fc, 9) is gapped for all 9 E [0, tt]. 
In other words, the topological space of all ID insulating 
Hamiltonians h(k, 9) is connected, which is a consequence 
of the Wigner-Von Neumann theorem'^'*. 

Suppose h{k, 9) is such a "gapped interpolation" be- 
tween hi{k) and h2{k). In general, h{k,9) for 9 £ (0,7r) 
doesn't necessarily satisfy the particle-hole symmetry. 
For 9 e [tt, 2tt], define 



h{k, 9)^- {C-^h{~k, 27r - 9)Cy 



(34) 



We choose this parameterization so that if we replaced 
by a momentum wavevector then the correspond- 
ing higher dimensional Hamiltonian would be particle- 
hole symmetric. Due to the particle-hole symmetry of 
h{k,9 = 0) and h{k,9 = tt), h{k,9) is continuous for 
9 S [0, 27r], and /i(fc,27r) = h{k,0). Consequently, the 
adiabatic evolution of 9 from to 27r defines a cycle of 
adiabatic pumping in h(k,9), and a first Chern number 
can be defined in the (fc, 9) space. As discussed in Sec. 
Ill C) the Chcrn number C[h{k, 9)] can be expressed as a 
winding number of the polarization 



C[h{k,9)] = 
Pi9) = 



d9 



dP{9) 
89 



E^(k)<0 



where the summation is carried out over the occu- 
pied bands. In general, two different parameterizations 
h{k,9) and h'{k,9) can lead to different Chern num- 
bers C[h{k,9)] ^ C[h'{k,9)]. However, the symmetry 
constraint in Eq. ([34|) guarantees that the two Chern 
numbers always differ by an even integer: C[h{k, 9)] — 
C[h'{k,9)] =2n, ne Z. 



To prove this conclusion, we first study the behavior of 
P{9) under a particle-hole transformation. For an eigen- 
state \k,9;a) of the Hamiltonian h{k,9) with eigenvalue 
Ea{k,9), Eq. dSll) leads to 

hi-k, 2tt - 9)C |fc, 9-a)* = ~Eo,{k)C \k, 9; a)* (35) 

in which \k,9;a)* is the complex conjugate state: 
|fc, 9; a)* = J^m.p ((™' /^l ^' ^' ^))* I™' Z^) where m, /3 are 
the position space lattice, and orbital index respectively. 
Thus C\k,9;a)* = \—k,2TT — 9;a) is an eigenstate of 
h{-k, 271-9) with energy i;a(fc, 27r- 6*) = -Ea{k,9) and 
momentum —k. Such a mapping between eigenstates of 
h{k, 9) and h{—k, 27r — 9) is one-to-one. Thus 

^(^) = E i-^){k,0;a\dk\k,9;a) 



-Ea(fc)<0 



= f ^ E {-^)(.{-k,2n-9;a\ 
■dk\-k,27r-9;a)* 



^P{2tt 



(36) 



Since P{9) is only well-defined modulo 1, the equality 
((M)) actually means P{9) + P{2tt - 9) = mod 1. Con- 
sequently, for = or TT we have 2-k — 9 = 9 mod 2-k, 
so that P{9) = or 1/2. In other words, the polariza- 
tion P is either or 1/2 for any particle- hole symmetric 
insulator, which thus defines a classification of particle- 
hole symmetric insulators. If two systems have different 
P value, they cannot be adiabatically connected without 
breaking the particle-hole symmetry, because P (mod 1) 
is a continuous function during adiabatic deformation, 
and a P value other than and 1/2 breaks particle- hole 
symmetry. Though such an argument explains physi- 
cally why a Z2 classification is defined for particle-hole 
symmetric system, it is not so rigorous. As discussed in 
the derivation from Eq. (PT|) to Eq. the definition 

P{9) = § dkak/2-K relies on a proper gauge choice. To 
avoid any gauge dependence, a more rigorous definition 
of the Z2 classification is shown below, which only in- 
volves the gauge invariant variable dP{9)/d9 and Chern 
number Ci. 

To begin with, the symmetry ([55)1 leads to 



dP{9) 



dPi9), 



(37) 



which is independent of gauge choice since only the 
change of P{9) is involved. This equation shows that 
the change of polarization during the first half and the 
second half of the closed path 9 S [0, 2tt] are always the 
same. 

Now consider two different parameterizations h{k,9) 
and h'{k,9), satisfying h{k,0) = h'{k,0) = hi{k), 
h{k,Tr) = /i'(fc,7r) = h2{k). Denoting the polarization 
P{9) and P'{9) corresponding to h{k,9) and h'{k,9), re- 
spectively, the Chern number difference between h and 
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C[h] - C[h'] 



C[gi]+C[g2] 



FIG. 4: Illustration of the interpolation between two particle- 
hole symmetric Hamiltonians /ii(fc) and h2{k). 



h' is given by 



C[h]-C[h'] = 



'\o('-^-'-K^'u (38) 







V 09 89 
Define the new interpolations gi{k, 9) and .92(^1 9) as 



9i{k,9) = 
92{k,9) = 



h{k,9), 6le[0,7r] 

/i'(fc,27r - 6*), 6* e [7r,27r] 

h'{k,2TT-9), 9e[0,T:] 

h{k,9), 9e[TT,2n] 



(39) 



f/i(fc, 9) and g2{k, 9) are obtained by recombination of the 
two paths h{k,9) and h'{k,9), as shown in Fig. IH From 
the construction of 51 and 32, it is straightforward to see 
that 



C[9i] 
C[92] 



d9 



(dPl9) dP'{9) 



\ 89 



89 
8P'{9) 
89 



(40) 



Thus C[h] - C[h'] = C[gi] + C[g2]. On the other 
hand, from Eq. ((37)) we know C[gi\ — C[g2\, so that 
C[h] - C[h'] = 2C[gi\. Since C[gi\ G Z, we obtain that 
C[/i] — C[/i'] is even for any two interpolations h{k, 9) and 
/i'(fc, 9) between hi[k) and h2{k). Intuitively, such a con- 
clusion simply comes from the fact that the Chern num- 
ber C\h] and C[h'] can be different only if there are sin- 
gularities between these two paths, while the positions of 
the singularities in the parameter space are always sym- 
metric under particle-hole symmetry, as shown in Fig. 

SI 

Based on the discussions above, we can define the "rel- 
ative Chern parity" as 



JVi[/li(fc),/l2(fc)] = (-l)^[''('='^)l, 



(41) 



which is independent of the choice of interpolation 
h{k,9), but only determined by the Hamiltonians 
hi{k), h2{k). Moreover, for any three particle-hole sym- 
metric Hamiltonians hi{k), h2{k), h3{k), it is easy to 



prove that the Chern parity satisfies the following asso- 
ciative law; 

Ni[hi{k),h2{k)]Ni[h2{k),h3{k)]^Ni[hi{k),h3{k)]. 

(42) 

Consequently, Ni[hi{k), h2{k)] — 1 defines an equiva- 
lence relation between any two particle-hole symmetric 
Hamiltonians, which thus classifies all the particle-hole 
symmetric insulators into two classes. To define these 
two classes more explicitly, one can define a "vacuum" 
Hamiltonian as ho{k) = ho, where Hq is an arbitrary 
matrix which does not depend on k and which satisfies 
the particle-hole symmetry constraint C^hoC = — ft-o ■ 
Thus ho describes a totally local system, in which there 
is no hopping between different sites. Taking such a 
trivial system as a reference Hamiltonian, we can define 
Ni[ho{k) , h{k)] = Ni[h{k)] as a Z2 topological quantum 
number of the Hamiltonian h{k). All the Hamiltonians 
h{k) with Ni[ho{k), h{k)] = 1 are classified as Z2 trivial, 
while those with Ni[hQ{k), h{k)] — —1 are considered as 
Z2 nontrivial. (Again, this classification doesn't depend 
on the choice of "vacuum" ho, since any two vacua are 
equivalent.) 

Despite its abstract form, such a topological character- 
ization has a direct physical consequence. For a Z2 non- 
trivial Hamiltonian /ii(fc), an interpolation h{k,9) can 
be defined so that h{k,0) — ho, /i(/e,7r) — hi{k), and the 
Chern number C[h{k, 9)] is an odd integer. If we study 
the one-dimensional system h{k, 9) with open boundary 
conditions, the tight binding Hamiltonian can be rewrit- 
ten in real space as 



1 

7i 



k 



h{k,9), VI < m,n < L. 



As discussed in Sec. IlICi there are mid-gap end states 
in the energy spectrum of /imn(^) as a consequence of 
the non-zero Chern number. When the Chern num- 
ber C[h{k,9)] = 2n — 1, n G Z, there are values 
9^ e [0,27r), s = l,2,..2n- 1 for which the Hamiho- 
nian hmn{9s) has zero energy localized states on the left 
end of the ID system, and the same number of 9f values 
where zero energy states are localized on the right end, 
as shown in Fig. [Sj Due to the particle-hole symme- 
try between hmn(9) and /im„(27r — 9), zero levels always 
appear in pairs at 9 and 2'k — 9. Consequently, when 
the Chern number is odd, there must be a zero level at 
9 ~ or 9 = TT. Since 9 = corresponds to a trivial 
insulator with flat bands and no end states, the localized 
zero mode has to appear at 9 = n. In other words, one 
zero energy localized state (or an odd number of such 
states) is confined at each open boundary of a Z2 non- 
trivial particle-hole symmetric insulator. 

The existence of a zero level leads to an important 
physical consequence — a half charge on the boundary of 
the nontrivial insulator. In a periodic system when the 
chemical potential vanishes, the average electron density 



on each site is n„ 



N when there 
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(b) 



By using Eq.. (|37|) we obtain 




It In 



vacuum odd insulator vacuum 
e=0 8=ji 6=0 



(c) 



PA 



P(Hl) 




2 V/P(H2) 




FIG. 5: (a) Schematic energy spectrum of a parameterized 
Hamiltonian hmnifi) witli open boundary conditions. The 
red (blue) lines indicate the left (right) end states. The 6 val- 
ues with zero-energy left edge states are marked by the solid 
circles, (b) Illustration to show that the open boundary of a 
Z2 nontrivial insulator is equivalent to a domain wall between 
9 = n (nontrivial) and 6 = (trivial vacuum), (c) Illustra- 
tion of the charge density distribution corresponding to two 
different chemical potentials fii (red) and fi2 (blue) . The area 
below the curve p{fii) and p(/i2) is +e/2 and — e/2, respec- 
tively, which shows the half charge confined on the boundary. 



are N bands filled. In an open boundary system, define 
Pm(M) = {J2a 4na<^rna)^ - TV to be the density deviation 
with respect to N on each site. Then particle-hole sym- 
metry leads to PmUj) = ~Pm{—^J■)■ On the other hand, 
when fj, is in the bulk gap, the only difference between fi 
and — /i is the filling of the zero levels localized on each 
end |0L) and |0i?), so that 



lim {pm{p) - P,n{-P)) 



\{ma\OL)\^ 



for the sites m that are far away enough from the right 
boundary. Thus we have J2mPm(f^ ~^ ^~^) ~ ^/^ where 
the summation is done around the left boundary so that 
we do not pick up a contribution from the other end. In 
summary, a charge e/2 (—e/2) is localized on the bound- 
ary if the zero level is vacant (occupied), as shown in Fig. 

m 

The existence of such a half charge can also be under- 
stood by viewing the open boundary of a topologically 
nontrivial insulator as a domain wall between the non- 
trivial insulator and the vacuum. By defining the inter- 
polation hmn{S), such a domain wall is described by a 
spatial dependence of with 9{x — > -l-oo) = tt, 0{x — > 
—00) — 0. According to the response formula (^5]) . the 
charge carried by the domain wall is given by 



dx 



(43) 



2 /„ dno) = i^c[h{kM- 



(44) 



It should be noted that an integer charge can always be 
added by changing the filling of local states, which means 
Qd is only fixed modulo e. Consequently, a ±e/2 charge 
is carried by the domain wall if and only if the Chern 
number is odd, i.e., when the insulator is nontrivial. 



E. Z2 classification of (0 -f l)-d particle-hole 
symmetric insulators 

In the last subsection, we have shown how a Z2 clas- 
sification of (1 + l)-d particle- hole symmetric insulators 
is defined by dimensional reduction from (2 -I- l)-d sys- 
tems. Such a dimensional reduction can be repeated 
once more to study (0 + l)-d systems, that is, a single- 
site problem. In this subsection we will show that a Z2 
classification of particle-hole symmetric Hamiltonians in 
(0 -|- l)-d is also obtained by dimensional reduction. Al- 
though such a classification by itself is not as interesting 
as the higher dimensional counterparts, it does provide 
a simplest example of the "dimensional reduction chain" 
(2 + l)-d^ (1 + l)-d^ (0 + l)-d, which can be later 
generalized to its higher-dimensional counterpart (4+1)- 
d^ (3 -I- l)-d^ (2 -I- l)-d. In other words, the Z2 classifi- 
cation of the (0 -I- l)-d particle- hole symmetric insulators 
can help us to understand the classification of (2 + l)-d 
TRI insulators as it is dimensionally reduced from the 
(4 + l)-d TRI insulator. 

For a free, single-site fermion system with Hamiltonian 
matrix h, the particle-hole symmetry is defined as 



C^hC 



(45) 



Given any two particle-hole symmetric Hamiltonians hi 
and /12, we follow the same procedure as the last sub- 
section and define a continuous interpolation h{9), £ 
[0, 2tt] satisfying 

h{Q) = hi, /i(7r) = /i2, Ct/i(6l)C = -/i(27r- 0)'^, (46) 

where h{9) is gapped for all 9. The Hamiltonian h{6) is 
the dimensional reduction of a (l-l-l)-d Hamiltonian ft.(fc), 
with the wavevector k replaced by the parameter 9. The 
constraint (|46|) is identical to the particle-hole symmetry 
condition so that h{9) corresponds to a particle- hole 
symmetric (1 -|- l)-d insulator. As shown in last subsec- 
tion, h{9) is classified by the value of the "Chern parity" 
Ni[h{9)]. If Ni[h{9)] = —1, no continuous interpola- 
tion preserving particle-hole symmetry can be defined be- 
tween h{9) and the vacuum Hamiltonian h{d) = hQ,y9 G 
[0,27r]. To obtain the classification of (0 -I- l)-d Hamil- 
tonians, consider two different interpolations h{9) and 
h'{9) between hi and /i2- According to the associative 
law dH]), we know Ni[h{9)]Ni[h'{9)] = Ni[h{9) , h' (9)] , 



where Ni[h{6) , h' (9)] is the relative Chern parity be- 
tween two interpolations. In the following we will prove 
Ni[h{9) , h' {9)] = 1 for any two interpolations h and h' 
satisfying condition (|^^ . As a result, Ni[h{9)] is inde- 
pendent of the choice of interpolation between hi and 
/i2, so that NQ[hi,h2] = Ni[h{9)] can be defined as a 
function of hi and /i2- The Z2 quantity A'o defined for 
(0 -I- l)-d Hamiltonians plays exactly the same role as 
Ni[h{k) , h' (k)] in the (1 + I)-d case, from which a Z2 
classification can be defined. 

To prove Ni[h{9), h'{6)] = 1 for any two interpolations, 
first define a continuous deformation g{9, ip) between h{9) 
and h'{9), which satisfies the conditions below: 

g{9,^^Q) - h{e), g{e,v = n) = h'{9) 

g{0,(p) = hi, g{TT,ip) = h2 
C^g{e,^)C = -g(2^-0,27r-^)^. (47) 

From the discussions in last subsection it is easy to con- 
firm that such a continuous interpolation is always pos- 
sible, in which g{9, <p) is gapped for all 9 and Lp. In 
the two-dimensional parameter space 9, tp one can define 
the Berry phase gauge field and the first Chern number 
Ci[g{9,ip)\. By the definition of the Chern parity, we 
have Ni[h{9),h'{9)] = (_i)Ci[s(0,¥')]. However, the pa- 
rameterized Hamiltonian g{9, p) can be viewed in two dif- 
ferent ways: it not only defines an interpolation between 
h{9) and h'{9), but also defines an interpolation between 
(7(0, p}) = hi and giir, p) = h2- Since g{0, p) and giir, p) 
are "vacuum Hamiltonians" without any p dependence, 
they have trivial relative Chern parity, which means 
Ni[h{9),h'{9)]=Ni[g{%p),g{iT,p)]^Ni[hi,h2] = 1. 

In conclusion, from the discussion above we have 
proved that any two interpolations h{9) and h'{9) belong 
to the same Z2 class, so that the Chern parity Ni\h{9)] 
only depends on the end points hi and /12. Consequently, 
the quantity A'o[/ii,/i2] = -^i[^(^)] defines a relation be- 
tween each pair of particle- hole symmetric Hamiltonians 
hi and /12. After picking a reference Hamiltonian /iq, 
one can define all the Hamiltonians with NQ[hQ,h] = 1 
as "trivial" and iVo[/io,^] = —1 as nontrivial. The 
main difference between this classification and the one 
for (1 -I- l)-d systems is that there is no natural choice of 
the reference Hamiltonian /iq- In other words, the names 
"trivial" and "non-trivial" only have relative meaning in 
the (0 -I- l)-d case. However, the classification is still 
meaningful in the sense that any two Hamiltonians with 
iVo[ft.i, /12] = — 1 cannot be adiabatically connected with- 
out breaking particle-hole symmetry. In other words, the 
manifold of single-site particle-hole symmetric Hamilto- 
nians is disconnected, with at least two connected pieces. 

As a simple example, we study 2x2 Hamiltonians. A 
general 2x2 single-site Hamiltonian can be decomposed 
as 

3 

h = da(7° + dacr" (48) 

a=l 

where (T° = I and a^-'^'^ are the Pauli matrices. When 




FIG. 6: Illustration of the 2x2 single-site Hamiltonians. Each 
point on the sphere represents an unit vector d — d/\d\, and 
the north and south poles correspond to the particle-hole sym- 
metric Hamiltonians hi^2 ~ ±cr^, respectively. The blue path 
shows an interpolation between hi and /i2 satisfying the con- 
straint (|46p . which always encloses a solid angle fl = 2tt. 



C — , particle-hole symmetry requires C^hC — ^h'^, 
from which we obtain do — di = d2 = 0. Thus h — d^a^ , 
in which da ^ so as to make h gapped. Conse- 
quently, we can see that the two Z2 classes are simply 
(ia > and da < 0. When an adiabatic interpolation 
h{9) = dQ{9)a° + Y.ada{0)cT'' is defined from > to 
< 0, the spin vector d{9) has to rotate from the north 
pole to the south pole, and then return along the im- 
age path determined by the particle- hole symmetry (j46p , 
as shown in Fig. [S] The topological quantum number 
A^o[^i, ^2] is simply determined by the Berry's phase en- 
closed by the path da{9), which is tt when hi and /i2 are 
on different poles, and otherwise. From this example 
we can understand the Z2 classification intuitively. In 
Sec. IVBI we show that the Z2 classification of (2 -I- l)-d 
TRI insulators — the class that corresponds to the QSH 
effect — is obtained as a direct analog of the (0 + l)-d case 
discussed above. 



III. SECOND CHERN NUMBER AND ITS 
PHYSICAL CONSEQUENCES 

In this section, we shall generalize the classification 
of the (2 + l)-d TRB topological insulator in terms of 
the first Chern number and the (2 -I- l)-d Chern-Simons 
theory to the classification of the (4-f l)-d TRI topological 
insulator in terms of the second Chern number and the 
(4 -I- l)-d Chern-Simons theory. We then generalize the 
dimensional reduction chain (2+I)-d^ (f+I)-d^ (0+1)- 
d to the case of (4 + l)-d^ (3 + l)-d^ (2 + l)-d for TRI 
insulators. Many novel topological effects are predicted 
for the TRI topological insulators in (3 + l)-d and (2 + 1)- 
d. 
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det [[idt - ^o™) - KL^'^- ] (50) 



which determines the response of the fermionic system 
through the equation 



FIG. 7: The Feynman diagram that contributes to the topo- 
logical term (|52p . The loop is a fermion propagator, and the 
wavy lines are external legs corresponding to the gauge field. 



Second Chern number in (4 + l)-d non-linear 
response 



(51) 



In this subsection, we wiU show how the second Chern 
number appears as a non-hnear response coefficient of 
(4+ l)-d band insulators in an external U[\) gauge field, 
which is in exact analogy with the first Chern number as 
the Hall conductance of a (2 + l)-d system. To describe 
such a non-linear response, it is convenient to use the 
path integral formalism. The Hamiltonian of a (4 -f l)-d 
insulator coupled to a U{\) gauge field is written as 

H[A] = ^ (4„C''„e'^^""-c„0 + h.c) 

m.n 

(49) 



In the case of the (2 + l)-d insulators, the ef- 
fective action SeS contains a Chern-Simons term 
iCi/47r)Af,e>"'^d^Ar as shown in Eq. Q of Sec. 111X1 
in which the first Chern number Ci appears as the co- 
efficient. For the (4 -I- l)-d system, a similar topological 
term is in general present in the effective action, which 
is the second Chern-Simons term: 



'S'cff — 



C2 

247r2 



(52) 



The effective action of gauge field A^ is obtained by the 
following path integral: 



D[c]D[c'']e' 



-H[A]\ 



where /i, v, p,a,T — 0, 1, 2, 3, 4. As shown in Refs. |33||40|, 
l4ll . the coefficient C2 can be obtained by the one-loop 
Feynman diagram in Fig. [71 which can be expressed in 
the following symmetric form: 



C2 
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7^ 



-Tr 



G 



G 



dG- 



G 



dG- 



G 



dG_ 

dq° 



G 



dG_ 



(53) 



in which = (w, fci, fc2, /ca, ^4) is the frequency- 
momentum vector, and G{q^) = [uj + iS — h{ki)] ^ is the 
single-particle Green's function. 

Now we are going to show the relation between C2 
defined in Eq. (j53p and the non-abelian Berry's phase 
gauge field in momentum space. To make the statement 
clear, we first write down the conclusion: 

• For any (4 -I- l)-d band insulator with single parti- 
cle Hamiltonian /i(k), the non- linear response coef- 
ficient C2 defined in Eq. ([55]) is equal to the sec- 
ond Chern number of the non-abelian Berry's phase 
gauge field in the BZ, i.e.: 

C2 ^ ^ /rf"fce*^"'=^tr[/„/fc,] (54) 
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with !t[ - 



ft 



Q/3 



Q/3 



ia/3 



a/3 



(k) 



-z(a,k| — |/?,k) 



where «, j, k,£ = 1, 2, 3, 4. 

The index a in a"^ refers to the occupied bands, there- 
fore, for a general multi-band model, a"'^ is a non-abelian 
gauge field, and f"j^ is the associated non-abelian field 
strength. Here we sketch the basic idea of Eq. (|54p . 
and leave the explicit derivation to Appendix [B] The 
key point to simplify Eq. (|53p is noticing its topologi- 
cal invariance i.e. under any continuous deformation of 
the Hamiltonian ft.(k), as long as no level crossing oc- 
curs at the Fermi level, C2 remains invariant. Denote 
the eigenvalues of the single particle Hamiltonian ft,(k) 
as ea(k),a = 1,2, ...,7V with ea{k) < ea+i(k). When M 
bands are filled, one can always define a continuous de- 
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Generic Insulator 



Flat Band Model 



FIG. 8: Illustration showing that a band insulator with arbi- 
trary band structure ei(fc) can be continuously deformed to 
a fiat band model with the same eigenstates. Since no level 
crossing occurs at the Fermi level, the two Hamiltonians are 
topologically equivalent. 



formation of the energy spectrum so that (k) ~> ec for 
a < M and ea(k) — > for a > M (with ce > ec), while 
all the corresponding eigenstates |a,k) remain invariant. 
In other words, each Hamiltonian ft.(k) can be contin- 
uously deformed to some "flat band" model, as shown 
in Fig. [51 Since both Eq. ([55]) and the second Chcrn 
number in Eq. (j54p are topologically invariant, wc only 
need to demonstrate Eq. ([M]) for the flat band models, 
of which the Hamiltonians have the form 

/io(k) = ec ^ |a,k) (a,k| + es 

l<a<M P>M 

EE eGPG(k) + e£P£;(k). (55) 

Here Pg(}^) {Pe{}^)) is the projection operator to the 
occupied (un-occupied) bands. Non-abelian gauge con- 
nections can be defined in terms of these projection op- 
erators in a way similar to Ref. [12. Correspondingly, 
the single particle Green's function can also be expressed 
by the projection operators Pg, Pe^ and Eq. can 
be proved by straight-forward algebraic calculations, as 
shown in Appendix [Bl 

In summary, we have shown that for any (4-|-l)-d band 
insulator, there is a (4 -I- l)-d Chern-Simons term ([5^ 
in the effective action of the external U{1) gauge field, 
of which the coefficient is the second Chern number of 
the non-abelian Berry phase gauge field. Such a relation 
between Chern number and Chern-Simons term in the 
effective action is an exact analogy of the TKNN formula 
in (2-|-l)-d QH effect. By applying the equation of motion 
([5T|) . we obtain 



(56) 



which is the non-linear response to the external field . 
For example, consider a field configuration : 

A^ ^Q, Ay^ B,x, A, = ~E,t, A^=At=Q (57) 

where x, y, z, w are the spatial coordinates and t is time. 
The only non- vanishing components of the field curvature 



are F^y = Bz and Fzt = —E^, which according to Eq. 
generates the current 



47r2 



BzE, 



If we integrate the equation above over the x, y dimen- 
sions (with periodic boundary conditions and assuming 
Ez is does not depend on (x, y)), we obtain 



dxdyjy 
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dxdyBz Ez = 



^^Ez (58) 



where N^y = / dxdyBz /^-k is the number of flux quanta 
through the xy plane, which is always quantized to be 
an integer. This is exactly the 4DQH effect proposed in 
Ref. [ij. Thus, from this example we can understand 
a physical consequence of the second Chern number: In 
a (4 -|- l)-d insulator with second Chern number C2, a 
quantized Hall conductance C^^xyl'^T^ in the zw plane 
is induced by magnetic field with flux 2'KNxy in the per- 
pendicular {xy) plane. 

Similar to the (2-1- l)-d case, the physical consequences 
of the second Chern number can also be understood bet- 
ter by studying the surface states of an open-boundary 
system, which for the (4 -I- l)-d case is described by a 
(3 -I- l)-d theory. In the next subsection we will study 
an explicit example of a (4-1- l)-d topological insulator, 
which helps us to improve our understanding of the phys- 
ical picture of the (4 -I- l)-d topology; especially, after 
dimensional reduction to the lower-dimensional physical 
systems. 



TRI topological insulators based on lattice 
Dirac models 



In section IIIB[ we have shown that the model intro- 
duced in Ref. [33 realizes the fundamental TRB topolog- 
ical insulator in (2 -I- l)-d, and it reduces to the Dirac 
model in the continuum limit. Generalizing this con- 
struction, we propose the lattice Dirac model to be the 
realization of the fundamental TRI topological insulator 
in (4 + l)-d. Such a model has also been studied in the 
field theory literature^i^. The continuum Dirac model 
in (4 + l)-d dimensions is expressed as 

H= ( d'^x [il}'^{x)r {-id^) i;{x) + mV'^rV] (59) 



with i = 1,2,3,4 the spatial dimensions, and r^,/i = 
0,1, ..,4 the five Dirac matrices satisfying the Clifford 
algebra 



{r^,r'^} = 2Vi 



(60) 



with I the identity matrix^. 

The lattice (tight-binding) version of this model is 
written as 



H 



E 



V4 ( ^ ) ^n+i + h.C. 
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n 

or in momentum space, 



(61) 



^ sin kiV^ + + c ^ cos fc,; j r° 



(62) 



Such a Hamiltonian can be written in the compact form 
H = 5I4'^a(k)rVk (63) 



with 



da(k) = 




:y , sin kz , sin 



a five-dimensional vector. Similar to the (2 + l)-d two- 
band models we studied in Sec. IIIBl a single parti- 
cle Hamiltonian with the form /i(k) = da(k)T°- has two 
eigenvalues E± (k) ~ ± y^J2a '^aO^)i but with the key dif- 
ference that here both eigenvalues are doubly degenerate. 
When J2a^a.O^) = '^^(k) is non- vanishing in the whole 
BZ, the system is gapped at half-filling, with the two 
bands with E = -E_(k) filled. Since there are two occu- 
pied bands, an SU{2) x U{1) adiabatic connection can be 
defined^i^i^. Starting from the Hamiltonian ([S^ . one 
can determine the single particle Green's function, and 
substituting it into the expression for the second Chern 
number in Eq. (j53|l . We obtain 



C2 = ^ I d^kt'''"''''dad.,dl,dydAddd^de (64) 

which is the winding number of the mapping (ia(k) = 
da(k)/|d(k)| from the BZ to the sphere S'^ and 



a,b,c,d,e = 0,1,2,3,4. More details of this calculation 
are presented in Appendix O 



Since the winding number (j64p is equal to the sec- 
ond Chern number of the Berry's phase gauge field, it is 
topologically invariant. It is easy to calculate C2 in the 
lattice Dirac model (pT|) . Considering the lattice Dirac 
model with a fixed positive parameter c and tunable mass 
term to, C2(to) as a function of to can change only if the 
Hamiltonian is gapless, i.e., if X^a'^aC^'"^) = ^ ^'^^ some 
k. It's easy to determine that C2(to) = in the limit 
TO +00, since the unit vector (ia(k) (1,0,0,0,0) 
in that limit. Thus we only need to study the change of 
C2(to) at each quantum critical points, namely at critical 
values of to where the system becomes gapless. 

The solutions of equation J^a '^0(^1 to) = lead to five 
critical values of to and corresponding k points as listed 
below: 



-4c, k= (0,0,0,0) 
-2c, keP[(7r, 0,0,0)] 
TO = < 0, k e P [(tt, 0, TT, 0)] (65) 
2c, k e P[(7r,7r,7r,0)] 
4c, k — (tt, tt, tt, tt) 
in which P[k] stands for the set of all the wavevectors 
obtained from index permutations of wavevector k. For 
example, P[(7r, 0, 0, 0)] consists of (tt, 0,0,0), (0,7r, 0,0), 
(0, 0,7r,0) and (0,0, 0,7r). As an example, we can study 
the change of C2(to) around the critical value to = —4c. 
In the limit m -|- 4c <C 2c, the system has its minimal 
gap at k = 0, around which the (ia(k) vector has the 
approximate form da(k) ~ (6m,kx,ky,kz,kw) -|- o(|fc|), 
with 5m = m + 4c. Taking a cut-off A <C 27r in momen- 
tum space, one can divide the expression (|64p of C2 into 
low-energy and high-energy parts: 



StT^ \J\k\<A J\k\>A J 



Since there is no level-crossing in the region |k| > A, the 
jump of C2 at Sm — can only come from Cj^^. In the 
limit \6m\ < A <^ 27r, the continuum approximation of 
(ia(k) can be applied to obtain 

which can be integrated and leads to 

AC2t=^=AC(^),„^o_=l. (66) 



dydcdJddn,de = C^^\5m,K) + ^((5to, A). 

I 

From the analysis above we see that the change of the 
second Chern number is determined only by the effec- 
tive continuum model around the level crossing wavevec- 
tor(s). In this case the continuum model is just the Dirac 
model. Similar analysis can be carried out at the other 
critical to's, which leads to the following values of the 



second Chern number: 



C2{m) 




m < —4c or m > 4c 
-4c < m < -2c 
-2c < m < 
< m < 2c 
2c < m < 4c 



(67) 



A more general formula is given in Ref. l4Ct 

After obtaining the second Chern number, we can 
study the surface states of the topologically nontrivial 
phases of this model. In the same way as in Sec. Ill CI 
we can take open boundary conditions for one dimen- 
sion, say, w, and periodic boundary conditions for all 
other dimensions, so that k^^ky, k^ are still good quan- 
tum numbers. The Hamiltonian is transformed to a sum 
of ID tight-binding models: 



k.w 



k.w 



sinfc,T* 



ijjAw + 1) + h.c. 



c^cosfc, r° 



V'fe(w) 

(68) 



in which k = {kx, ky, kz), i = 1, 2, 3, and w = 1,2, .., L are 
the w coordinates of lattice sites. The single- particle en- 
ergy spectrum can be obtained as Ea{k), a = 1, 2, .AL, 
among which the mid-gap surface states are found when 
C2 0, as shown in Fig. [51 When the Chern number is 
C2, there are IC2I branches of gapless surface states with 
linear dispersion, so that the low energy effective theory 
is described by IC2I flavors of chiral fermions^: 



H = sgn(C2) 



^v^^pjipja -pi^.ip). (69) 



The factor sgn(C2) ensures that the chirality of the sur- 
face states is determined by the sign of the Chern num- 
ber. From such a surface theory we can obtain a more 
physical understanding of the nonlinear response equa- 
tion ([5^)1 to an external U{1) gauge field. Taking the 
same gauge field configuration as in Eq. (|57[) . the non- 



vanishing components of the field curvature are Fxy = 
and Fzt = —E^. Consequently, the (3-|-l)-d surface states 
are coupled to a magnetic field Ti = Bzi and an electric 
field E = Ez2. For simplicity, consider the system with 
—4c < m < —2c and C2 = 1, in which the surface theory 
is a single chiral fermion with the single particle Hamil- 
tonian 



h = Vd- Aj = V(JxPx+V<Jy {Py + BzX)+V(7z [pz ^ Ezt) . 

If Ez is small enough so that the time-dependence of 
Az{t) — —Ezt can be treated adiabatically, the single 
particle energy spectrum can be solved for a fixed Az as 



En±iPz 

Eoijpz 



±v^(p, + Azf 4 
v{pz + A^)sgn(B2 



2n\B._ 



1.2, 



(70) 




(b) 




FIG. 9: Three dimensional energy spectrum of the surface 
states for the parameters (a) c = 1, m = —3 and (b) c = 
1, m = — 1. The Dirac points where energy gap vanishes are 
marked with deepest red color. For m = — 3 there is one Dirac 
point at r point while for m = — 1 there are three of them at 
X points. 



When the size of the surface is taken as X Ly X Lz with 
periodic boundary conditions, each Landau level has the 
degeneracy N^y — LxLyBz/2'K. Similar to Laughlin's 
gauge argument for QH edge states'*, the effect of an in- 
finitesimal electric field Ez can be obtained by adiabati- 
cally shifting the momentum pz Pz + E^t. As shown 
in Fig. [H from the time t = to t = T = 2tt/L^Ez, the 
momentum is shifted as Pz —^ Pz + 2tt/Lz, so that the 
net electron number of the surface 3D system increases 
by Nxy In other words, a "generalized Hall current" 1^ 
must be flowing towards the w direction: 



N. 



T 



LxLyLzB^Ez 
4^ 



This "generalized Hall current" is the key property of the 
4DQH effect studied in Ref. [ij. 

In terms of current density we obtain — BzEz/4,tt'^ , 
which is consistent with the result of Eq. discussed 
in the last subsection. More generally the current density 
jw can be written as 

EB C2 



Jv 



Co 



47r2 
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^y-^ err 



(71) 



which is the chiral anomaly equation of massless (3-1- l)-d 
Dirac fermions'*^'^®. Since the gapless states on the 3D 
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a. 




FIG. 10: Illustration of the surface Landau level spectrum 
given by Eq. (|70p . Each level in the figure is Nxy fold degen- 
erate. The solid circles are the occupied states of the zeroth 
Landau level, and the red open circle shows the states that 
are filled when the gauge vector potential is shifted adia- 
batically from to 27r/Lz. 



edge of the 4D lattice Dirac model are chiral fermions, 
the current carries away chiral charge, leading to the 
non-conservation of chirality on the 3D edge. 



(72) 



Now consider a special "Landau" -gauge configuration 
satisfying A , - — A , ~ Vn, which is transla- 

tionally invariant in the w direction. Thus, under peri- 
odic boundary conditions the w-direction momentum kw 
is a good quantum number, and the Hamiltonian can be 
rewritten as 



H[A] = 



k-u} ,x,s 



4 



cT" - 



-I- (m -f ccos {k^ + Ast^)) r°] iix,k^ 

where x stands for the three-dimensional coordinates, 
Ag/^ = Aj^s+u); and s = 1,2,3 stands for the x,y,z di- 
rections. In this expression, the states with different k^ 
decouple from each other, and the (4-1- l)-d Hamiltonian 
H[A] reduces to a series of (3 -|- l)-d Hamiltonians. Pick 
one of these (3 + l)-d Hamiltonians with fixed kw and 
rename kw + Ag^ = 9^, we obtain the (3 -I- l)-d model 



IV. 



DIMENSIONAL REDUCTION TO (3 
TRI INSULATORS 



1)-D 



As shown in Sec. Ill C[ one can start from a (2 -I- 1)- 
d TRB topological insulator described by a Hamilto- 
nian h{kx,ky), and perform the procedure of dimen- 
sional reduction by replacing ky by a parameter 9. The 
same dimensional reduction procedure can be carried out 
for the (4 + l)-d TRI insulator with a non- vanishing 
second Chern number. From this procedure, one ob- 
tains the topological effective theory of insulators in 
(3 -I- l)-d and (2 -I- l)-d. Specifically, for TRI insula- 
tors a general Z2 topological classification is defined. 
Compared to the earlier proposals of the Z2 topological 
invarian t-'^^i^^i^'^i^^i^'^i^^i^^'^" , our approach provides a di- 
rect relationship between the topological quantum num- 
ber and the physically measurable topological response of 
the corresponding system. We discuss a number of the- 
oretical predictions, including the TME effect, and pro- 
pose experimental settings where these topological effects 
can be measured in realistic materials. 



A. Effective action of (3 + l)-d insulators 

To perform the dimensional reduction explicitly, in the 
following we show the derivation for the (4 -I- l)-d Dirac 
model (ISTj) . However, each step of the derivation is ap- 
plicable to any other insulator model, so the conclusion 
is completely generic. 

The Hamiltonian of Dirac model coupled to an 
external U{1) gauge field is given by 



4n 



:F° - iT' 



Vx-Hs + h.c. 



+ E ^5 \-^^^ ^^^"^ + (m + c cos 9s) F°] ^5- 

X,S 

(73) 

which describes a band insulator coupled to an electro- 
magnetic field As^3+s and an adiabatic parameter field 

Due to its construction, the response of the model ([75)1 
to Ag^gj^s and 9g fields is closely related to the response 
of the (4-f l)-d Dirac model ([HI]) to the f7(l) gauge field. 
To study the response properties of the (3-1- l)-d system, 
the effective action ^sDi^, 9] can be defined as 



exp 



£)[^]£)[^]e* / dt[Y.^ i,^{id^ -~A^o)i' — H[Afi] 



A Taylor expansion of S^d can be carried out around 
the field configuration As{x,t) = 0, 9(x,t) = 9q, which 
contains a non-linear response term directly derived from 
the (4 + l)-d Chern-Simons action ([5^ : 



S-go = / d^xdte^''"'^59d^A^d„Ar. 

An J 



(74) 



Compared to the Eq. the field 69 (x, t) = e{x, t) - 9o 
plays the role of A4, and the coefficient G3{9o) is de- 
termined by the same Feynman diagram ([7]), but evalu- 
ated for the three-dimensional Hamiltonian ([75)1 . Con- 
sequently, 03(^0) can be calculated and is equal to Eq. 
(j53p . but without the integration over k^: 
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(Pkdu 



G 



dG- 



G 



dG- 



G 



dG- 



G 



dG- 



G 



dG-^ 



(75) 



where ^ {uj,kx, ky, kz) . Due to the same calculation 
as Sec. nil Al and Appendix iBl G^IOq) is determined from 
the Berry phase curvature as 



(76) 



in which the Berry phase gauge field is defined 
in the four-dimensional space {kxj ky^ fc^, ^0)1 i.G.^ 



Q/3 



i (k, 6q; a 



(d/dk,) 



k,do;P) and a\ 



ilk, 60; a 



{d/d9o) k,9o;(3J. Compared to the second 

Chern number ([5i)) . we know that 03(^0) satisfies the 
sum rule 



GMdeo = C2eZ, 



(77) 



which is in exact analogy with the sum rule of the pump- 
ing coefficient Gi{9) in Eq. ([22]) of the (1 + l)-d system. 
Recall that Gi{d) can be expressed as dPi{d) /dd, where 
Pi{0) is simply the charge polarization. In comparison, 
a generalized polarization ^3(^0) can also be defined in 
(3 + l)-d so that G3(6'o) = dP^{9Q)/deo. ( Recently, a 
similar quantity has also been considered in Ref. 49 from 
the point of view of semiclassical particle dynamics.) The 
conventional electric polarization P couples linearly to 
the external electric field E, and the magnetic polariza- 
tion M couples linearly to the magnetic field B, however, 
as we shall show, P3 is a pseudo-scalar which couples non- 
linearly to the external electromagnetic field combination 
E-B. For this reason, we coin the term "magneto-electric 
polarization" for P3. To obtain P^{9q), one needs to in- 
troduce the non-Abelian Chern-Simons term: 



ABCD 



IGtt^ 



Tr 



Ibc ~ o [«S'«c] 



, (78) 



which is a vector in the four-dimensional parameter space 
q ~ {kx,ky,kz,9o) and A,B,C,D = x,y,z,9. JC^ satis- 
fies 



1 



OZTT 

GM = I d^kdAlC^. 



When the second Chern number is nonzero, there is an 
obstruction to the definition of a a, which implies that 
)Ca cannot be a single-valued continuous function in 
the whole parameter space. However, in an appropri- 
ate gauge choice, /C*,z = x,y,z can be single- valued, so 
that G3(6'o) = j(PkdeJC^ = dP3{9o)/d9o, with 



^3(^0) = / d^klC 
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fi] - o l-'^i 



(79) 



Thus, ^3(^0) is given by the integral of the non-Abelian 
Chern-Simons 3-form over momentum space. This is 
analogous to the charge polarization defined as the in- 
tegral of the adiabatic connection 1-form over a path in 
momentum space. 

As is well-known, the three-dimensional integration of 
the Chern-Simons term is only gauge-invariant modulo 
an integer. Under a gauge transformation Ui — > u~^aiU — 
iu~^diU {u € U{M) when M bands are occupied), the 
change of P3 is 
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j Ske^^^^Ti [{u-^diu) {u-^dju) {u-^dku)] , 



which is an integer. Thus Pi,{9o), just like Pi{9), is only 
defined modulo 1, and its change during a variation of 9q 
from to 27r is well-defined, and given by C2. 

The effective action ((74|) can be further simplified by 
introducing G3 = dPj,/d9. Integration by parts of Sj,y) 
leads to 



ftn — — 



-5- / cPxdte^""'-' A^{dPy,ld9)dJ9d^Ar. 
Ait J 

{dPz/d9)d^59 can be written as d^P^.^ where P'i{x,t) = 
P3{9{x, t)) has space-time dependence determined by the 
9 field. Such an expression is only meaningful when the 
space-time dependence of 9 field is smooth and adiabatic, 
so that locally 9 can still be considered as a parameter. 
In summary, the effective action is finally written as 



<S'3D = -r- 



J d^xdte^'"'^P3ix,t)d^A,d^Ar. 



(80) 



This effective topological action for the (3 -I- l)-d insu- 
lator is one of the central results of this paper. As we 
shall see later, many physical consequences can be di- 
rectly derived from it. It should be emphasized that this 
effective action is well-defined for an arbitrary (3 + l)-d 
insulator Hamiltonian h{k, x, t) in which the dependence 
on X, t is adiabatic. We obtained this effective theory by 
the dimensional reduction from a (4 + l)-d system; and 
we presented it this way since we believe that this deriva- 
tion is both elegant and unifying. However, for readers 
who are not interested in the relationship to higher di- 
mensional physics, a self-contained derivation can also be 
carried out directly in (3 -I- l)-d, as we explained earlier, 
by integrating out the fermions in the presence of the 
Afj_{x,t) and the 9{x,t) external fields. 
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This effective action is known in the field theory ht- 
erature as axion electrodynamica^i^'^, where the adia- 
batic field P3 plays the role of the axion field^i^l. When 
the P3 field becomes a constant parameter independent 
of space and time, this effective action is referred to as 
the topological term for the 9 vacuum^i^. The axion 
field has not yet been experimentally identified, and it 
remains as a deep mystery in particle physics. Our work 
shows that the same physics can occur in a condensed 
matter system, where the adiabatic "axion" field Psix, t) 
has a direct physical interpretation and can be accessed 
and controlled experimentally. 



From the discussion above it is clear that 3D TRI 
topological insulators realize a non-trivial solitonic back- 
ground 9 field. In Ref. the authors suggest several 
candidate materials which could be 3D topological insu- 
lators. These 3D materials are topologically non-trivial 
because of band inversion mechanism similar to that of 
the HgTe quantum wells^. Ref. 113 closely studied the 
strained, bulk HgTe. We will keep this system in mind 
since it has a simple physical interpretation, and its es- 
sential physics can be described by the Dirac model pre- 
sented earlier. We can consider the trivial vacuum out- 
side the material to have a constant axion field with the 
value 9 = and the interior of a 3D topological in- 
sulator to have a 9 = tt background field. The value 
9 = TT does not violate time-reversal (or CP in high- 
energy language) . HgTe is a zero-gap semiconductor and 
has no topologically protected features. However, when 
strained, the system develops a bulk insulating gap be- 
tween the p-wave light-hole "conduction band" and the 
p-wave heavy-hole "valence band" around the F-point. 
To study the topological features we must also include 
the s-wave band which in a conventional material like 
GaAs would be a conduction band. Because of the strong 
spin-orbit coupling in HgTe the band structure is actually 
inverted the s-wave band becomes a valence band. For 
a moment we will ignore the heavy-hole band and only 
consider the light-hole and s-wave band^. The effective 
Hamiltonian of these two bands is a massive Dirac Hamil- 
tonian, but with a negative mass. The negative mass in- 
dicates a phase shift of tt in the vacuum angle 9 from its 
original unshifted value in the trivial vacuum. The ax- 
ion domain wall structure at the surface of the topologi- 
cal insulator traps fermion zero modes which are simply 
the topologically protected surface states. If we include 
the effects of the heavy-hole band the dispersion of the 
bulk bands and surface states are quantitatively modi- 
fied. However, as long as the crystal is strained enough 
to maintain the bulk gap the topological phenomena will 
be unaffected and the boundary of the 3D topological in- 
sulator can still be described as an axion domain wall. 
Thus, this material in condensed matter physics provide 
a direct realization of axion electrodynamics. 



B. Physical Consequences of the Effective Action 

5*313 

In this subsection we present the general physical 
consequences of the effective topological action for 
(3 + l)-d insulators coupled to a P3 polarization, and in 
subsection II V CI we focus on its consequences for TRI in- 
sulators. Since the effective action is quadratic in A^, 
it describes a linear response to the external electromag- 
netic fields which depends on the spatial and temporal 
gradients of P3 . Taking a variation of S^u [A, 9] we obtain 
the response equation: 

r = ^e^^'^'^dMaAr. (81) 
Ztt 

The physical consequences Eq. (|8T|) can be understood 
by studying the following two cases. 

(1) Hall effect induced by spatial gradient of P3. 

Consider a system in which P3 = Psiz) only depends 
on z. For example, this can be realized by the lattice 
Dirac model ([73|l with 9 = 9{z). (This type of domain 
wall has also been considered in Ref. fssi). In this case 
Eq. ([5T|) becomes 

Ztt 

which describes a QH effect in the xy plane with the 
Hall conductivity axy = dzP3/27T, as shown in Fig. [11] 
(a). For a uniform electric field in the a;-direction, 
the Hall current density is jy = {dzPd./'2,n)Ex- Thus the 
integration over z in a finite range gives the 2D current 
density in the xy plane: 

In other words, the net Hall conductance of the region 
zi < z < Z2 is 

'^Iv^ dP3/27T, (82) 

J zi 

which only depends on the change of P3 in this region, 
and is not sensitive to any details of the function P3(z). 
Analogously in the (1 -I- 1) — d case, if we perform the 
spatial integration of Eq. we obtain the total charge 
induced by the charge polarization P: 

Q = - I dP/2TT. (83) 

By comparing these two equations, we see that the rela- 
tion between P3 and Hall conductance in (3 -I- l)-d insu- 
lators is the same as the relation between charge polar- 
ization P and the total charge in the (1 -I- l)-d case. As 
a specific domain wall between two homogeneous 

materials with different P3 will carry Hall conductance 
(Jh — AP3/27r, while the fractional charge carried by a 
domain wall in (1 -t- l)-d is given by Q = — AP/27r. 
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FIG. 11: (a) Illustration of the Hall effect induced by a spatial 
gradient of P3. The colors represent different values of P3, 
which decrease from at the bottom to 1/2 on top. The blue 
arrow shows the direction of a uniform electric field Ex and 
the white arrows show the Hall current density induced, given 
by the formula jy — {dzP2,/2TT)Ex- (b) Energy spectrum of 
the (3 + l)-d lattice Dirac model (|73p in a magnetic field 
towards the z direction. The boundary conditions are periodic 
in the x and t/-directions and open on the z-direction. The 
red and black curves show the surface states on the top and 
bottom surfaces, respectively, each of which is a B^L^Ly /2n 
fold degenerate Landau Level. The parameters of model (|73|) 
are chosen to be m = —3, c = 1. 



(2) Topological Magneto-electric effect(TME) 
induced by temporal gradient of P3. 

When P3 = Pj,{t) is spatially uniform, but time- 
dependent, Eq. ([5T|) becomes 



27r ' 



"-'^djAk, i,j,k = x,y,z. 



In other words, we have 



J = 



dtP3 
2tt 



B. 



(84) 



Since the charge polarization P satisfies j — dtP, 
in a static uniform magnetic field B we have dtP — 
-dt [PzB /2t:^, so that 



P = -—(P3 + const.) 



(85) 



Such an equation describes the charge polarization in- 
duced by a magnetic field, which is a magneto-electric 
effect. Compared to similar effects in multiferroic 
materials^i^, the magneto-electric effect obtained here 
is of topological origin, and only determined by the 
magneto-electric polarization P3. 

Similar to the (1 + l)-d adiabatic pumping effect, the 
response Eq. ([84| can be understood in a surface state 
picture. For example, consider the lattice Dirac model 
([75)1 with periodic boundary conditions in the x, y direc- 
tions and open boundary conditions in the z-direction. 
In the presence of a static magnetic field B^ in the z- 
direction, the single particle energy spectrum En{0) can 
be solved for at a fixed 9 value. As shown in Fig. [11] 



(b), mid-gap states appear for generic 9, which are local- 
ized on the (2 -I- l)-d boundary. It should be noticed that 
each state is A^-fold degenerate where TV = BzLxLy/2TT is 
the Landau level degeneracy. In the lattice Dirac model, 
when —4c < m < —2c so that C2 — Jg^Q^ dP^ = 1, we 
find that during a period 9 = 0-^ 2tt, N degenerate 
surface states on the bottom boundary sink below fermi 
level and N states on the top float up. Consequently, 
when 9 is adiabatically tuned from to 2tt, there are 
N electrons pumped from the top surface to the bottom 
one, which is in consistent with the result of Eq. 



AQ 



dt / dxdyj;. 



2-n 



-NCo 



Just like the relation between (2 + l)-d QH edge states 
and the mid-gap end states in the (1 -I- l)-d pumping 
effect, there is a direct relationship between the (3 -I- l)-d 
surface states of the (4 -|- l)-d lattice Dirac model and 
the adiabatic pumping discussed above. As discussed in 
Sec. IIIIBl the surface theory of a (4-1- l)-d lattice Dirac 
model with nontrivial C2 is a (3 + l)-d chiral fermion. 
As shown in Fig. 1101 the energy spectrum in a magnetic 
field Bz has a chiral dependence on the wavevector pz- 
During the dimensional reduction procedure, Pz (in the 
notation of Sec. IIIIBp is replaced by the parameter 9, 
so that the chiral energy spectrum E{pz) changes to the 
"chiral" 9 dependence of En{9) in Fig. [TT] (b). In other 
words, the adiabatic pumping in the magnetic field in the 
(3-|- l)-d system is a dinicnsionally reduced version of the 
chiral anomaly on the surface of a (4 + l)-d topological 
insulator. 

The TME leads to a striking consequence if magnetic 
monopoles are present. For a uniform P3, Eq. leads 
to 



ZTT 



B. 



Suppose we consider a compact U{1) electromagnetic 
field on a lattice, where the monopole density pm = 
V • B 121: can be non- vanishing, then we obtain 



dtPe = {dtPz) P-n 



(86) 



Therefore, when P3 is adiabatically changed from zero to 
0/27r, the magnetic monopole will acquire a charge of 



ZTT 



(87) 



Such a relation was first derived by Witten in the context 
of the topological term obtained from QCD^. 



Z-i topological classification of time-reversal 
invariant insulators 



In Sec. Ill Dl we have seen how a Z-i topological clas- 
sification is obtained for (1-1- l)-d particle-hole symmet- 
ric insulators. The key point for that case is to show 
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that any interpolation between two particle-hole sym- 
metric insulators hi (fc) and /12 (k) carries the same parity 
of Chern number, so that the "relative Chern parity" is 
well-defined for each two Hamiltonians with particle-hole 
symmetry. In this section, wc will show that the same 
approach can be applied to (3-f l)-d insulators, where the 
time-reversal symmetry plays the same role as particle- 
hole symmetry does in (1 -I- l)-d. 

For a Hamiltonian H = J2m.n^niah'^n'^nf3, the time- 
reversal transformation is an anti-unitary operation de- 
fined by Cma ^ T'^^Cm0, where the time-reversal matrix 
T satisfies T^T = I and T*T = -I. In fc-space time- 
reversal symmetry requires 



The condition T*T = —I is essential, and leads to 
Kramers's degeneracy. Now wc will follow the same ap- 
proach as Sec. HID I and show how to define a Z2 invariant 
for the TRI insulators in (3 + l)-d. For any two TRI band 
insulators hi(k) and /i2(fc), an interpolation h{k,6) can 
be defined, satisfying 



h{k,0) = hl{k), h{k,TT) = h2{k) 

T'^h{-k,-e)T = h^{k, e), 



(89) 



and h{k,6) is gapped for any 6 g [0, 27r]. Since the 
interpolation is periodic in 0, a second Chern number 
C2[h{k, 6)] of the Berry phase gauge field can be defined 
in the {k,9) space. In the same way as in Sec. IIIDi we 
will demonstrate below that C2[h{k,9)] — C2[h'{k,0)] — 



mod 2 for any two interpolations h and h' . First of all, 
two new interpolations gi^2(k,9) can be defined by Eq. 
P9[) . which we repeat here for convenience: 



h{k,e), 6ie[0,7r] 

h'{k,2TT-e), e [tt,2tt] 

h'{k,2TT-9), 6ie[0,7r] 

h{k,e), e e [7r,27r] ' 



By their definition, gi and g2 satisfy C2[h] — C2[h'] — 
C2L91] +C2[g2] and T^gi{~k,~9)T = g^{k,e). To 
demonstrate C2[gi] — C2[g2], consider an eigenstate 
k,9;a\ of gi{k, 0) with eigenvalue Ea{k, 6). We have 
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g^{-k,-e)T^ k,0;aj^ = T^gi{k, 

= E^{k,0)T^ 
{^k,-0)T^ (^k,e;a^y ^ Eaik,e)T^ (^k,e;a^y 



fc, 9; a 



Thus T'^{\k, 6 -,0)1)* is an eigenstate oi g2{—k,—9) with 
the same eigenvalue Ea{k,6). Expand over the eigen- 

states I — k,— 6,13)2 of g2{—k,—9), we have 



-k,-9;(3 



(90) 



Consequently the Berry phase gauge vector of the gi and 
g2 systems satisfies 



a?/(fc,0) 



i (k,9;a 



k,9;(3 



■y.S 



-k,-9;j 



dj (Ufjs 



~k,-9;6 



^ {/„^a]f (-fc, -9){U^),^ - t ^ [/„^(fc, 9)d,u;^{k, 

7,(5 7 



(91) 



In other words, a"j{k,9) is equal to a^j {—k, —9) up to 
a gauge transformation. Consequently, the Berry phase 
curvature satisfies fi^{k,6) = Uajf]f* {—k, ~9){U^)gfj, 
which thus leads to C2[gi{k,9)] = C2[g2{k,9)]. In 
summary, we have proved C2[h{k,9)] — C2[h'{k,9)] — 
'^C2[g{k, 9)] = mod 2 for any two symmetric interpola- 
tions h and h' . Thus the "relative second Chern parity" 

iV3[/ii(fc),/l2(fc)] = {-lf-^>^i^fi)\ 

is well-defined for any two time-reversal invariant (3-1-1)- 
d insulators, independent on the choice of interpolation. 
In the same way as in (1-1- l)-d, a vacuum Hamiltonian 
/lo(fc) = ho, Vfc can be defined as a reference. All the 



a/3/ 



Hamiltonians with N3[ho,h] = — 1 are called Z2 nontriv- 
ial, while those with N^lho, h] = 1 are trivial. 

Similar to the (1 -I- l)-d case, there is a more intu- 
itive, but less rigorous, way to define the Z2 invariant N3. 
Through the derivation of Eq. (PT|) one can see that for a 
TRI Hamiltonian satisfying Eq. ([55]) . the Berry's phase 
gauge potential satisfies ai{k) = Uai{—k)W — iUdiW , 
so that the magneto-electric polarization P3 satisfies 



2R 



247r2 



(fke'^^TY [{Ud.U'^) {UdjU'^) (f/9fe[/^)] e Z. 



Consequently, there are only two inequivalcnt, TRI val- 
ues of P3, which are i-3 = and P3 = 1/2. For 
two Hamiltonians hi and /12, the second Chern number 
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C2[h{k, 9)] = 2 (P3[/i2] - Palhi]) mod 2, so the difference 
of P3 determines tlie relative Cfiern parity N^lhi, /i2] by 
N3[hi,h2] ^ (-l)nP3lhi]-P4h2])^ Since the trivial Hamil- 
tonian ho obviously has P3 = 0, we know that all the 
Hamiltonians with P3 = 1/2 are topologically non-trivial, 
while those with P3 = are trivial. 

Once the Z2 classification is obtained, the physical con- 
sequences of this topological quantum number can be 
studied by the effective theory ([80]) . as has been done in 
the last subsection. In the (1 -I- l)-d case, we have shown 
that a zero-energy localized state exists at each open 
boundary of a Z2 nontrivial particle-hole symmetric insu- 
lator, which leads to a half charge Qd — e/2( mod e) con- 
fined on the boundary. Similarly, the nontrivial (3-1-1)- 
d insulators also have topologically protected surface 
states. The easiest way to study the surface physics of the 
(3 -|- l)-d insulator is again by dimensional reduction. As 
discussed above, for any three-dimensional Hamiltonian 
/ii(fc), an interpolation h(k,9) can be defined between 
hi and the "vacuum Hamiltonian" ho. If we interpret 
6 as the fourth momentum, h{k,9) defines a (4 + 1)- 
d band insulator. Moreover, the constraint Eq. ((5^ 
on h{k, 9) requires time-reversal symmetry for the corre- 
sponding (4 + l)-d system. The Hamiltonian h{k, 9) can 
be written in a real space form and then defined on a four- 
dimensional lattice with open boundary conditions in the 
^-direction and periodic boundary conditions for all the 
other directions. As discussed in Sec. IIII Bl there will be 
|C2[/i]| flavors of (3 -I- l)-d chiral fermions on the surface 
when the second Chern number C2 [h] is nonzero. In other 
words, in the 3D BZ of the surface states there are |C2[/i]| 
nodal points (kxn,kyn,9n) , n — 1, .., |C2[/i]| where the 
energy spectrum En{kx,ky,9) is gapless and disperses 
linearly as a Dirac cone. From time-reversal symmetry 
it is easy to prove that the energy spectrum is iden- 
tical for {kx,ky,9) and {—kx,—ky,—9). Consequently, 
if {kx,ky,9) is a nodal point, so is {—kx,—ky,—9). In 
other words, time-reversal symmetry requires the chiral 
fermions to appear in pairs, except for the ones at time- 
reversal symmetric points, as shown in Fig. [T2](a). Thus, 
when the second Chern number C2[h] is odd, there must 
be an odd number of Dirac cones at the 8 symmetric 
points in the 3D BZ. Actually, the {4 + l)-d lattice Dirac 
model (|62l) provides an example of TRI insulators with 
nontrivial second Chern number, since one can define 
to be time-reversal even and F^'^'^'"' to be odd, as in con- 
ventional relativistic quantum mechanics^-. As shown 
in Fig. [21 all the nodal points of the surface states are 
located at the symmetric points F, M, R or X . 

Now we return to the surface of (3 4- l)-d insulator. 
Since hi{k) = h{k,0), ho = h{k,Tr) by definition of the 
interpolation, the surface energy spectra of hi and ho 
are given by the 9 = and 9 = n slices of the 3D surface 
spectrum. Since all 8 time-reversal symmetric points (F, 
X and M) are at = or = tt, we know that the net 
number of Dirac cones on the surface energy spectrum of 
hi and ho is odd (even) when C2[h{k,9)] is odd (even). 



(a) 



(b) 





A / 








/ i 

/ A 





M 




(3+1 )d surface of 
(4+1 )d system 



(2+1 )-d surface of 
(3+1 )d system 



FIG. 12: Illustration of the nodal points in the surface state 
energy spectrum of (4-1- l)-d and (3-f l)-d insulators, (a) The 
nodal points in the (fca,, fc^, BZ, for a 2: = const, surface of 
the (4 + l)-d system. The red (blue) points stand for nodal 
points at time-reversal symmetric (asymmetric) wavevectors. 
The dashed lines are a guide for the eyes. There are two 
pairs of asymmetric nodal points and three symmetric points 
in this example, which correspond to a bulk second Chern 
number C2 = 7. (b) The nodal points in the (k^,ky) BZ, 
for a 2: = const, surface of the (3 -I- l)-d system. According to 
the dimensional reduction procedure (see text) the 2D surface 
energy spectrum is given by the 6 = slice of the 3D surface 
spectrum in (a) . Since 6 = ir corresponds to a vacuum Hamil- 
tonian Hq, no nodal points exist in that plane. Consequently, 
the number of nodal points in the 2D BZ (5 in this example) 
has the same parity as C2 in the (4 4- l)-d system. 



However, ho is defined as the vacuum Hamiltonian, which 
is totally local without any hopping between different 
sites. Thus, there cannot be any mid-gap surface states 
for ho- Consequently, the number of 2D Dirac cones in 
the surface state spectrum of hi must be odd (even) when 
C2 [h] is odd (even) . Since the parity of C2 [h] determines 
the Z2 invariant N^lhi], we finally reach the conclusion 
that there must be an odd (even) number of (2 -I- 1)- 
d gapless Dirac fermions confined on the surface of a 
(3 + l)-d nontrivial (trivial) topological insulator. 

Compared to earlier works on Z2 invariants and surface 
states in (3-|-l)-d, one can see that the Z2 nontrivial topo- 
logical insulator defined here corresponds to the "strong 
topological insulator" of Ref. [l^- The present approach 
has the advantage of (i) demonstrating the bulk-edge re- 
lationship more explicitly, (ii) clarifying the connection 
between the second Chern number and the Z2 topological 
number and (iii) naturally providing the effective theory 
that describes the physically measurable topological re- 
sponse properties of the system. The "weak topological 
insulators" defined in Ref. [13 are not included in the 
present approach, since these Z2 invariants actually cor- 
respond to topological properties of (2 + l)-d insulators 
(QSH insulators, as will be discussed in next section), just 
as the QH effect in (3 -I- l)-d systems^ still corresponds 
to a first Chern number, but defined in a 2D projection 
of the 3D BZ. 
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D. Physical properties of Z2-nontrivial insulators 

In the last subsection we have defined the Zi topologi- 
cal quantum number for the (3-|-l)-d TRI insulators, and 
discussed the gapless Dirac fermions on the surface of a 
non-trivial insulator. Now we will study the physical re- 
sponse properties of the non-trivial insulators. Since the 
non-trivial insulator has a magneto-electric polarization 
P3 = 1/2 mod 1, according to Eq. ([50]) the effective ac- 
tion of the bulk system should be 

S3B = / d^^xdU^^^^^df^A^d^Ar. (92) 

Svr J 

in which n — P3 — 1/2 ^ Z, is the integer part of P3 . Un- 
der time-reversal symmetry, the term e^''°"^9pA,^9crA^ = 
2E • B is odd, so that for general P3, the effective ac- 
tion (j80[) breaks time-reversal symmetry. However, when 
the space-time manifold is closed {i.e., with periodic 
boundary conditions in the spatial and temporal dimen- 
sions), the term / d^xdte'^'''""^ dfj^A^daA-r is quantized to 
be STT^m, m G Z. Consequently, = (2ri + l)rmr so 
that the action e*'^^° = e™^ = (—1)™ is time-reversal 
invariant and is independent of n, the integer part of 
P3. This time-reversal property of the effective action is 
consistent with that of P3 discussed in the last subsec- 
tion. Thus, in a closed space-time manifold, the topo- 
logical action ([5^ is a consistent effective theory of the 
Z2-nontrivial TRI insulators, and the integer part of P3 
is not a physical quantity. However, the case is different 
when the system has open boundaries. On a space-time 
manifold with boundary, the value of Ssu is not quan- 
tized, which thus breaks time-reversal symmetry even for 
P3 = 1/2 or 0. In this case, the integer part n of P3 does 
enter the action e*"^^", and becomes physical; its value 
depends on the quantitative details of the boundary. 

To understand the physics in the open boundary sys- 
tem, we study a semi-infinite nontrivial insulator occupy- 
ing the space z < 0. Since the vacuum, which fills z > 0, 
is effectively a trivial insulator (with an infinitely large 
gap) , the effective action (|80|) can be written in the whole 
of R3 as 

-^30 = ^ j d^xdte^^^'^^A^dMaAr 

Since P3 = 1/2 mod 1 for z < and P3 = mod 1 for 
z > 0, we have 

d,P3 = {n + l/2)5{z), 

where n G depends on the non-topological details of 
the surface, as will be studied later. In this case, the 
effective action is reduced to a (2 + l)-d Chern-Simons 
term on the surface: 

5su.t = - ^^y^ / drcdydie^A'-P^ (93) 
Stt J 

This is consistent with the observation in Sec. lIVBI that 
a domain wall of P3 carries a QH effect. The Hall con- 
ductance of such a surface of a Z2 nontrivial insulator is 



thus aH = (n-t- l/2)/27r, which is quantized as a half odd 
integer times the quanta /h. On the other hand, from 
the discussion in the last subsection we know that there 
are always an odd number of (2-1- l)-d Dirac fermions liv- 
ing on the surface of a nontrivial insulator. Thus the half 
QH effect on the surface can be easily understood by the 
parity anomaly of massless Dirac fermions^. Here we 
need to be more careful. The Hall conductance carried 
by a Dirac fermion is well-defined only when the ferniion 
mass is non-vanishing, so that a gap is opened. With 
the continuum Hamiltonian H — k^i-a^ + kyU^ + mcr^, 
the Berry phase curvature can be calculated as in Eq. 
pT|) . The d vector is given by d = {kx,ky,m), which 
has a meron-type configuration in kx , ky space, and thus 
carries a Hall conductance^ 

1 

an = — sgn(TO) (= — sgn(m)). 

Now consider the surface of a topological insulator with 
2n+l gapless Dirac cones. From the discussion in the last 
subsection we know that without breaking time-reversal 
symmetry, at least one of these Dirac cones cannot be 
gapped. Now consider a perturbation that breaks time- 
reversal symmetry, i.e. a term which assigns a mass 
mi,i — 1,2, ..2n -|- 1 to each Dirac cone and induces a 
net Hall conductance cth — ^S^("^i)/47''- Since 

Si^i"^ sgn(TOi) is an odd integer, the Hall conductance 
we obtain is consistent with the surface Chern-Simons 
theory ([M]) . From this discussion we can understand 
that the effective action (|93p describes a surface with 
time-reversal symmetry breaking, though the bulk sys- 
tem remains time- reversal invariant. The bulk topology 
requires there to be a 1/2 quanta in the Hall conduc- 
tance, and the surface time-reversal symmetry breaking 
term determines the integer part n. This is an exact ana- 
log of the half charge on an end of the (1 -I- l)-d particle- 
hole symmetric insulator. As shown in Fig. [5l whether 
the localized state on the end of a nontrivial insulator 
is filled or vacant can only be determined by choosing a 
chemical potential ^1 > or /j,2 < 0, or equivalently, by 
breaking the particle-hole symmetry around the bound- 
ary. The charge localized at one end of the insulator 
chain is (n -t- l/2)e, in which the integer part n depends 
on the symmetry breaking term on the surface, but the 
1/2 part is guaranteed by the bulk topology. 

To show such a surface state picture more explicitly, 
again the lattice Dirac model can be taken as an example. 
Consider the (3 -I- l)-d lattice Dirac model ([75]) with a 
domain wall configuration of the 9{x) field given by 

e{x) = e{z) = I [1 - tanh(z/4e)] , (94) 

which has the asymptotic behavior 9{z —^ —00) — tt, 
6{z — > +00) = 0. The domain wall width is ^, as shown 
in Fig. [12] (a). With periodic boundary conditions in 
the X and j/-directions, the Hamiltonian can be block 
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FIG. 13: (a) Domain wall configuration 9{z) (blue line) and 
corresponding sin 6(z) (red dashed line) as defined by Eq. (|94p 
with 5 = 20. (b) Dispersion relation of the two-dimensional 
bands trapped on the domain wall for the TRl Hamiltonian 
ho- (c) The mass of surface Dirac fermion (blue line) and the 
half-quantized Hall conductance carried by the domain wall 
(red line) in the Hamiltonian h{\) — ho + Xhi as a function of 
A. The arrows on the left and right side show schematically 
the rotation of angle 9 across the domain wall for A < and 
A > 0, respectively. 



diagonalized: 



H 



^lfe„ (^) ( o'^ ) + 1) + h.c. 



+ ^ il^l^i^^iz) [{m + c cos 9{z) + c cos kx 

+ccosky) r° +smkxT^ + sinfcyr^] ipk^^kyiz) 
+ E <fc,Wsin0(z)rVfe.fe„W^i^o + i/i. (95) 

Under a time-reversal transformation, is even and 
pi, 2, 3, 4 odd. Thus the only time-reversal odd term 
in Hamiltonian ([95|) is the last term, which is localized 
around the boundary due to the factor sm9{z) (See Fig. 

(a)). Hereby we denote H = Hq + Hi with Hi the 
last term and Hq all the other TRI terms, and define 
ho, hi as the single particle Hamiltonian corresponding 
to Hq,Hi, respectively. Then the Hamiltonian iJp de- 
scribes a time-reversal invariant interface between two 
insulators = and 9 = tt. For —4c < m < —2c, the pa- 
rameterized Hamiltonian H{9), 9 e [0, 27r] has a Chern 
number C2 = 1, so that the 6 = and 9 = tt system has 
relative Chern parity —1. It's easy to show that Hamil- 
tonian ((73|) for 6* = 0, — Ac < m < —2c is adiabatically 
connected to the m — > — 00 limit. Thus we know that 
9 = Q and 9 = n correspond to Z2 trivial and nontriv- 
ial insulators, respectively. Consequently, on the domain 
wall at z = there should be an odd number of gapless 
Dirac cones for Hq. As shown in Fig. [T3] (b), numeri- 
cal diagonalization of ho shows one single Dirac cone at 
{kx,ky) = (0,0) on the surface. To understand the ef- 



fect of hi term, notice that {F"*, /iq} = 0, with { } being 
the anti-commutator. The effective Hamiltonian of the 
Dirac cone can always be written as /isurf = k^Ux ^y^y 
in a proper basis, and it should also anti-commutc with 
F'* since the bulk Hamiltonian does. Since the only term 
that anti-commutes with /isurf in the 2x2 Hilbert space is 
(Tz, we know that the effect of F"* term is to induce a mass 
term mUz in the lattice Dirac model. More accurately, 
the amplitude and the sign of m can be determined by 
standard perturbation theory. Given the two zero-energy 
surface states j/c = 0, a) , a = 1, 2, the representation of 
the matrices Ux and Uy in the effective theory /igurf can 
be determined by 



<f} = (fc = 0,a| 



dho 
dki 



\k = 0,/3) , i = x,y. 



k=0 



Then the is given by cr^ = —ia^a^, so that the mass 
m is determined by 

m = i ^ cr^^ (fc = 0, (3\ hi \k = 0,a). 

If we consider the parameterized Hamiltonian h = ho + 
Xhi, then the mass of the surface Dirac fermion is pro- 
portional to A for A ^ 0. As shown in Fig. [T3] (c), 
the mass is positive for A < 0, which leads to a sur- 
face Hall conductance an = — sgn(A)/47r. On the other 
hand, the surface Hall conductance can also be calcu- 
lated by the effective theory through Eq. (|82p . For 
A = 1, the phase field 9 winds from tt to 0, which leads to 
<JH = J^^dPs{z)/27r = J^dP3{9)/27: = -l/47r (since 
dPs = C2/2 = 1/2). The Hamiltonian for A = -1 
can be considered to be the same lattice Dirac Hamil- 
tonian H{9) with 9{z) replaced by —9{z). This keeps ho 
invariant but reverses the sign of hi. Consequently, the 
winding of the 9 field is from — tt to 0, which leads to 
a Hah conductance an = S°_^dP^{9)/2 = l/47r. The 
winding of 9 in the two cases is shown schematically in 
Fig. [11(c). 

In summary, from this example we learn that the effect 
of a time-reversal symmetry breaking term on the surface 
is to assign a mass to the Dirac fermions which deter- 
mines the winding direction of P3 through the domain 
wall. Once each Dirac cone on the surface gains a mass, 
the whole system is gapped and the Berry phase curva- 
ture is well-defined, so that the winding number of P3 
through the domain wall is determined. Physically, the 
time-reversal symmetry breaking term on the surface can 
come from magnetic fields or magnetic moments localized 
on the surface; it could also arise from the spontaneous 
breaking of time reversal symmetry on the surface due 
to interactions. Once such a "T-breaking surface field" 
(denoted by M) is apphed, the effective action ([M)) is 
well-defined for open boundaries, and describes the elec- 
tromagnetic response of the Z2 nontrivial insulator. Ac- 
tually, the T-breaking field should be considered to be an 
external field applied to the TRI system, such that the 
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topological action (|80|) describes a nonlinear response of 
the system to the combination of M and electromagnetic 
field A^. For a Z2 nontrivial insulator occupying a spa- 
tial region V with boundary 9V, the spatial gradient of 
P3 is given by 



VP3(x) = ( g[M{x)] + i 



dn{y)5\x^y) (96) 



dV 



where g[M{x)\ g Z is the integer part of the winding 
number determined by the T-breaking field M{x), and n 
is the normal vector of the surface. Under such a con- 
figuration of VP3, the effective action ((80)) is reduced to 
the surface Chern-Simons action 

Ssurf - ^ J^^dn^ (giMix)] + e^^^'A.d.Ar. (97) 

Since g[M{x)] can only take discrete values, in general 
the surface of a nontrivial insulator consists of several do- 
mains with different Hall conductance. To obtain more 
realistic predictions of the effective theory ([97|) . in the 
rest of this subsection we will study a specific case — the 
interface between a ferromagnetic insulator and a Z2 non- 
trivial insulator, where the surface time-reversal symme- 
try breaking is generated by the magnetization of the FM 
material. Several specific experimental proposals will be 
discussed. 

(1) Magnetization-induced QH effect. 

Consider the ferromagnet-topological insulator het- 
erostructure shown in Fig. [M) The magnetization of the 
two FM layers can be parallel or antiparallel, and the 
standard six-terminal measurement can be performed to 
measure the in-plane Hall conductance. The net Hall 
conductance is given by the summation of the contribu- 
tions of the top and bottom surfaces. When the topolog- 
ical insulator is uniform, an outward pointing magneti- 
zation vector, {i.e., towards the direction of the surface 
normal vector fi), will have the same effect, no matter to 
which surface it is applied. Suppose the Hall current on 
the top surface induced by an electric field E = 
is jt = rif X E/47r, then on the bottom surface the 
same formula applies, such that jf, = ftf, x E/47r. Since 
fit = —fib = z, the current jt = — jb, as shown in Fig. 
1141 Consequently, the antiparallel magnetization leads 
to a vanishing net Hall conductance, while the parallel 
magnetization leads to (Th = e} jh. 

Just like the usual integer QH effect, the quantized Hall 
conductance here is carried by chiral edge states. To un- 
derstand the edge state picture, notice that for parallel 
magnetizations in Fig. [T3](a), the magnetization vector is 
outward pointing at the top surface and inward pointing 
at the bottom surface. Although the Hall conductance of 
the two surfaces are the same in the global x,y,z basis, 
they are opposite in the local basis defined with respect to 
the normal vector fi. In other words, the integer g[M{x)] 
in the surface Chern-Simons theory is for the top 
surface and —1 for the bottom surface. Consequently, the 
side surface is a domain wall between two different QH 
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FIG. 14: Illustration of the QH effect in ferromagnet- 
topological insulator heterostructure. (a) and (b) The electric 
field Ex (with direction into the paper) and the induced Hall 
current jt and j;, for parallel and antiparallel magnetization, 
respectively. In (b) the Hall current on the two surfaces are 
opposite and form a circulating current, (c) A top-down view 
of the device for Hall measurement. The grey regions are leads 
that contact to the surface of the topological insulator, (d) 
For the case of parallel magnetization, the chiral edge states 
are trapped on the side surfaces of the topological insulator. 
These carry the quantized Hall current. 



regions with Hall conductances that differ by one quan- 
tum. Just like a domain wall between v = Q and v = 1 
regions in the usual QH system, such a domain wall will 
trap a chiral fermi liquid, which in this experimental pro- 
posal is responsible for the net Hall effect. It should be 
noticed that the side surface is 2D, so that generically 
there are also other non-chiral propagating modes on the 
side surface, besides the one branch of chiral edge states. 
However, the existence of these non-chiral states does not 
change the stability of the chiral edge state, since there is 
always one more right mover than left mover. The stabil- 
ity of the QH effect is still protected by the "bulk" gap, 
which is the magnetization-induced gap Em in this case. 
Thus we will expect to observe this QH effect under the 
following two requirements: (i) temperature ksT <C Em] 
(ii) the chemical potential on the top and bottom surfaces 
remains in the gap induced by the applied magnetization. 

We would like to point out that this experimental pro- 
posal provides a direct demonstration of the half QH 
effect on the surface of a topological insulator. If the 
Oh = /h measured for parallel magnetization were con- 
tributed by one surface, then the magnetization flip of the 
other surface would have no effect on the net Hall con- 
ductance. Thus if an /h Hall conductance is observed 
for parallel magnetization, and the magnetization flip of 
either magnet leads to vanishing Hall conductance, one 
can conclude that the Hall conductance is contributed 
equally by the two surfaces. 
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(2) Topological Magneto-Electric Effect 
(TME). 

As has been discussed in Sec. IIVBI a TME effect 
is induced by P3, which is described by Eqs. ([84|) and 
([85|) . Now we consider the reahzation of this effect in 
a nontrivial topological insulator. Similar to the surface 
QH effect, a magnetization (or any other time-reversal 
symmetry breaking term) is necessary to determine the 
integer part of P3. Consider the FM-TI-FM structure 
in Fig. [13] (b). With antiparallel magnetizations, the 
current induced by an electric field on the top and 
bottom surfaces flows in opposite directions. If we con- 
sider an isolated system rather than the Hall bar with 
leads as discussed above, a circulating current is formed, 
which induces a magnetic field parallel or anti-parallel 
to the electric field. However, in the geometry shown 
in Fig. (|14p (b), dissipation occurs when the circulating 
current flows on the gapless side surface and the adia- 
batic condition of the TME effect is violated. To obtain 
the TME effect, a T-breaking gap for the side surface is 
necessary. This is satisfied in the cylindrical geometry 
shown in Fig. [15] (a). With a magnetization pointing 
out of the cylinder's surface, the surface is gapped and 
has a fixed Hall conductance an = {n + \)e^ /h. In an 
electric field parallel to the cylinder as shown in Fig. [15] 
(a), a tangential circulating current is induced, with the 
strength jt = urE. The magnetic field generated by such 
a current in the topological insulator can be obtained by 
solving Maxwell's equations: 

Bt = cr^E = -(2n+l)— E (98) 

c he 

in which CGS units are reintroduced. From this for- 
mula we can see that the magnetic field induced by an 
electric field is proportional to the electric field, where 
the response coefficient is quantized in odd multiples of 
the fine structure constant. When the magnetization of 
the side surface is reversed, the magnetic field induced 
is also reversed, as is expected from time-reversal sym- 
metry. Combined with the conventional, non-topological 
response, we obtain B = H + 47rM — {2n + l)j-Ei, or 

H = B-47rM+(2n+l)— E (99) 

he 

Similar induction between electric and magnetic fields 
also occurs when a magnetic field is applied. Consider 
the process of applying a magnetic field B parallel to the 
cylinder as shown in Fig. [15] (b). When the magnetic 
field is turned on from zero, a circulating electric field 
parallel to the side surface is generated, which then in- 
duces a Hall current j oc dB/dt parallel or anti-parallel 
to the magnetic field. Consequently, a charge density 
proportional to B is accumulated on the top and bot- 
tom surfaces, so that a magnetic field induces an electric 
field parallel to it. By solving Maxwell equations, such a 
topological contribution to the electric field is obtained 




FIG. 15: (a) Illustration of the magnetic field induced by an 
electric field in a cylinder geometry. The magnetization of the 
FM layer points outward from the side surface of the TI, and 
a circulating current is induced by the electric field, (b) Illus- 
tration of the electric field induced by a perpendicular mag- 
netic field in the same geometry. "©" and "0" indicate the 
positive and negative charge induced by the magnetic field on 
the top and bottom surface, respectively, (c) Illustration of a 
topological insulator with fully gapped surface states induced 
by a hedgehog magnetization configuration. P3 of the topo- 
logical insulator is well-defined since the integral along any 
two pathes L and L' gives the same AP3 = P-iix) — Pi(xii). 
(d) Illustration of a topological insulator with parallel magne- 
tization on the surface. A one-d domain wall between two-d 
surface regions with different Hall conductance (as shown by 
red and blue) is induced, which carries chiral edge states as 
shown by (g) and 0. In this case the P3 of the topological 
insulator cannot be determined, since two different pathes L 
and L' lead to different AP3. 

as 

e2 

Et = (2n+1)— B. (100) 
he 

Combined with the conventional, non-topological re- 
sponse, we obtain: 

D = E + 47rP - (2n+ 1)— B (101) 

he 

The conventional Maxwell's equations, supplemented by 
the constituent relations (|10ip and give the complete 
description of the electrodynamics of the 3D topological 
insulators. 

An alternative description is to use the conventional 
constituent relations D = E -I- 47rP and H = B — 47rM, 
and work with a set of Maxwell's equations modified by 
the topological term. The total action of the electro- 
magnetic field including the topological term is given by 

dsni): 

'S'tot = 'S'Maxwoll + >5'topo 
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= / d'^xdt 



16^ + 2 ''''' ^ 



(102) 



in which a = /he is the fine structure constant, and 
poi ^ pi g^j^j -pij ^ Mfc are the electric and mag- 
netic polarization vectors, respectively. The equations of 
motion arc obtained by variation of the action over 
as 

+ d.V^^" + ^e^''^--d,{P:iF,r) = -f- (103) 
47r 47r c 

These equations can also be written in the more familiar 
component form as 

V • D = 47rp + 2a(VP3 • B) 

VxH-l^ = ^j-2a (^(VPs X E) + ^ (9tP3) B 
c 



V-B = 



(104) 



where D = E + 47rP and H = B - 47rM include only 
the non-topological contributions. These are the equa- 
tions of motion of axion electrodynamic o^°'^^i^^ . By shift- 
ing the topological terms to the left-hand side, and re- 
defining D and H according to (|10ip and ([M)) . and tak- 
ing P3 — n + we recover the conventional Maxwell's 
equations, but with modified constituent relations, thus 
demonstrating the equivalence to the formulation given 
above. The quantization of the TME effect in odd units 
of the fine structure constant is a deeply profound quan- 
tization phenomenon in condensed matter physics. The 
flux quantization inside a superconductor determines the 
fundamental constant hc/e, while the quantization of 
the Hall resistance determines the fundamental constant 
h/e^. To date, there has been no other known quan- 
tization phenomenon in units of the dimensionless fine 
structure constant a = /he. 

It should be emphasized that the value of P3 in the 
topological insulator can only be determined when a mag- 
netization is applied to open a gap on the surface. As 
shown in Fig. [15] (c), by defining a path L from a refer- 
ence point xq deep in the vacuum, P3 can be determined 
by Pz{x) = Jggj^dl ■ VP3. However, this definition only 
applies when the result does not depend on the choice 
of path. If a magnetic "shell" covered the surface of 
the topological insulator, with the magnetization outgo- 
ing everywhere on the interface, then the change of P3 
across the interface is the same for different points on 
the surface, such that the bulk P3 is well-defined without 
dependence on the choice of path. In this case, Eqs. ([M]) 
and (jlOip are well-defined, and the integer part of P3 can 
change if the magnetization direction is reversed. On the 
other hand, when there are domain walls on the surface, 
the integer part of P3 is not well-defined in the bulk of 
the topological insulator, and Eqs. ([M]) and poip do not 



apply, as shown in Fig. [15] (d). Physically, the failure of 
Eqs. and (|10ip is simply due to the existence of a 
QH edge current on the domain wall, which requires the 
more general Maxwell equations Eqs. (|104p including the 
contribution of the current. This analysis also provides 
a new picture of the surface QH effect, that is, the QH 
effect on the surface is carried by the chiral edge states 
living on vortex rings of the P3 field. It is only when 
there are no vortex rings of P3 on the surface, that the 
surface is fully gapped and the electro-magnetic response 
is simply given by Eqs. and (|10ip . 

(3) Low-frequency Faraday rotation. 

The TME effect can be experimentally observed in the 
settings discussed above, by applying an electric field 
through a capacitor, and measuring the magnetic field by 
a SQUID device. Alternatively, we consider the experi- 
ment of Faraday or Kerr rotation. The modified Maxwell 
equations (|104p can be applied to another phenomenon — 
photon propagation in the system^''. It should be noted 
that the effective theory (|102l) only applies in the low- 
energy limit E Eg, where Eg is the gap of the surface 
state. Thus, to detect the topological phenomena we 
should study the low frequency photons with lo <ti Eg/h. 
Consider a FM-TI interface at z = 0, as shown in Fig. [TH] 
Normally incident, linearly-polarized light can be written 



A{z,t) 



ae 



■( — kz—ujt) _j_ ^^i{kz—LL!t) 2 > 



^i{ — k'z—ujt) 



z < 



(105) 



in which k — lo/v and fc' — lo/v' are the wavevectors 
of the photon in the z > and z < region, re- 
spectively. The VP3 terms in Eq. (|104p contribute a 
non-conventional boundary condition at z = 0. Define 
VP3 = Az5(z) (with A - 1/2 e Z), the boundary condi- 
tions are given by 



a + b = c 
z X [A;(-a + b)//i + fc'c//i'] c 



in which the dimensionless constants e, e' and /i, are 
the permittivity and permeability of the z > and z < 
materials, respectively. Denote a± = a^; ± iay and the 
same for 6±, c±, the equations above lead to 



a+ 



fc'//.*' - 2ia/\uj/c 



k/^i 



Consequently, when the incident wave a is linearly polar- 
ized, the transmission wave c is also linearly polarized, 
with the polarization plane rotated by an angle 



9topo = arctan ■ 



2aA 



U/Jl + ^T/Jj/ 



(106) 



Here we always assume that the magnetization of the FM 
material is perpendicular to the xy plane, so that H = 
/iB holds for in-plane magnetic fields. In the simplest 
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(a) _ (b) 




Experimentally, the main difficulty of measuring such 
an effect comes from the low-frequency constraint lu <ti 
Eg/h. For a typical value Eg — lOmeV we get / = 
Lo jl-K <C 2.4THz, which is in the far infrared or microwave 
region. In principle, it is possible to find a topological 
insulator with a larger gap which can support an accurate 
measurement of Faraday rotation. Similar proposals as 
above can also be worked out for the rotation of reflected 
wave (Kerr effect). 



FIG. 16: (a) Illustration of the Faraday rotation Stopo on 
the interface between a ferromagnet (or equivalently, a para- 
magnet in magnetic field) and a topological insulator, (b) 
Illustration of the total Faraday rotation angle ^ as a func- 
tion of magnetic field B in the sandwich structure as shown 
in the inset. The zero-field extrapolation of Q(B) determines 
the topological term Stopo. 



case ~ 1, e, e' ~ 1 and A = 1/2, we get ~ a ~ 
7.3 X XQ-'^T&A. 

Since the ferromagnetic material itself also induces a 
Faraday rotation, its necessary to distinguish these two 
contributions in order to measure the topological contri- 
bution (|106p . Replace the FM layers by paramagnetic 
materials with large susceptibilities, and apply an ex- 
ternal magnetic field to polarize them. In this case the 
magnetization is proportional to magnetic field, such that 
the Faraday rotation contributed by the bulk is also pro- 
portional to magnetic field. The net Faraday rotation is 
given by = ^topo + ^topo + ^'buik, which has the following 
dependence on the magnetic field: 



Q{B) = u-B 2sgn(B) arctan ■ 



(107) 



Consequently, the topological contribution can be ob- 
tained by measuring d(S) at different applied magnetic 
fields and extracting the linear extrapolation of as 
B 0+. 



V. DIMENSIONAL REDUCTION TO (2 + 1)-D 

By carrying out the same dimensional reduction pro- 
cedure once more, we can obtain the topological effective 
theory for TRI (2 -I- l)-d insulators. Additionally a Z2 
classification can be defined for (2 -t- l)-d TRI insula- 
tors, which is in exact analogy to the Z2 classification of 
(0-f l)-d particle-hole symmetric insulators. We will show 
that the (2 -I- l)-d nontrivial phase corresponds to the 
QSH insulator proposed recentl y^^d'''^°i^^ , and study the 
physical consequences of the effective theory. 



Effective action of (2 + l)-d insulators 



In Sec. IIV Al we have seen how a (3 -I- l)-d insulator 
with a parameter field d{x,t) is related to a (4 -I- l)-d 
insulator through dimensional reduction. In the same 
way, two parameter fields can be defined to obtain the 
dimensional reduction from (4 -|- l)-d to (2 -I- l)-d. In 
the following we will still take the lattice Dirac model as 
a canonical example to show the dimensional reduction 
procedure and derive the effective theory. 

Starting from the lattice Dirac model ((72)) and 
choosing a special gauge vector configuration satisfying 
K,n+i = A^+ra.n+w+i = K+z,n+i+h (^^ that the gauge 
vector is homogeneous along z and w) we obtain the 
Hamiltonian 



kz .kuj ,x s— 1,2 



/cF"-iF" 
V 2 



h.c. 



fe^,fc„,xs=l,2 



• [sin (fc^ + ylxs) F'"^ -f sin {k^ -f A^^) -f (m -I- ccos {k^ + ^xs) + ccos (fc^ -I- A^^)) F"] ip^-k,.k^ 

I 

in which x = (x, y) is the two-dimensional coordinate. As family of (2 -t- l)-d Hamiltonians: 
in the (3-1- l)-d case, the gauge fields in the z and w direc- 
tions can be replaced by parameter fields (fcz-l-Axa) Oy, 



and (fcu, -|- Ax4) </'x; resulting in the parameterized H2B[A,9,ip\ = 



v4 



cF" - 



x+s 



h.c. 



-^^t [sin0,F3+sin(^xr^ 
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+ (m + ccos6'x + ccos ip^) r°] V'x- (108) in which the coefRcient G2(6'o, (fio) is determined by the 

same correlation function as Eq. (j53p . but without the 
By integrating out the fermion fields and expanding integrations over kz,k^: 
the resulting effective action around Ag — 0, 9 — do, </? = 
(po: the same nonlinear term shown in the Feynman dia- 
gram in Fig. [7] leads to the topological term 

52D = / (fxdte^"'PA^,^Je^p5^ (109) 

ZTT J 
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(110) 



in which /X, //. jO = 0, 1, 2, i,j = l,2, q^^ — {u},kx,ky), and 
the Berry curvature is defined in the four-dimensional 
parameter space {k^, ky, 0, ip). The coefficient G2 {0o, 'Po) 
satisfies the sum rule 



G2 {0Q, Vo) dOodipo = 27rC2. 



(Ill) 



To simplify the expression further, the Chern-Simons 
form /C-^ in Eq. ^ can be introduced again. Here 
A runs over kx,ky, 6, Lp, and G2{Oo, ipo) can be written in 
terms of K.^ as 

G2 (Oo, Vo) = -27r J d^k (a,/C" -I- dylC^ + deK." + S^/C^) 

Similar to the (3 -I- l)-d case, the momentum derivative 
terms 9(3; j,)/C*^^'^^ lead to vanishing contributions if /C'^'^-' 
is single- valued, in which case G2 can be expressed as 



G2 (do, fo) = defl^ - d^flo, 



with 



= -2tt I d^kJC\ Vtg ^2ti I d^kJCf. (112) 

Notice that 



1 



167r2 

2 ( fi^p 



gUXr 



[ai,aj 



and similarly for IC^. We know that the vector fl = 
(flg,^ip) has the correct transformation properties of a 
gauge vector potential under the coordinate transforma- 
tions of the parameter space (0, </?), and also under gauge 
transformations of the wave functions. Consequently, 
when the parameters 9 = 9{x^) and Lp = ip{x'^) have 



smooth dependence on space-time coordinates, an effec- 
tive gauge vector potential 17^ can be defined in (2 + l)-d 
space-time as 

= Qed^se + n^d^s^, (113) 

the gauge curvature of which is related to G2 as 

5^fi^ - dMfi = [de^^ - d^fle) {di_i_S0d^6(p - d^S0df_,S(p) 
= G2 {d^Sed.S^ ~ djed^S^) . (114) 

Mathematically, G2 is a density of second Chern 
form in the 2D parameter space {0,(j)), and {0,ip) = 
{d{xP), p{x'^)) defines a smooth map from the (2 -I- 1)- 
d space-time manifold to the 2D parameter space. The 
curvature of the gauge potential fl^ is the puUback of G2 
to (2 + l)-d space-time. By making use of Eq. (|113p the 
effective action (|109p can be rewritten in a parameter- 
independent form: 



52D = 7^ / d^xdte^"'^A^^MT■ 
2n 



(115) 



The physical consequences of the effective theory (|115p 
can be understood by studying the response equation: 



ZTT 



(116) 



As will be shown in the next subsection, Eq. (|116p is 
the fundamental response equation for the QSH effect, 
which takes the form similar to the fundamental response 
equation for QH effect (O, with the replacement of the 
external gauge field by a effective Berry's phase gauge 
field. In this sense, our formalism provides a unifying 
theory for both effects. This type of relationship between 
different types of topological insulators will be discussed 
in more detail in section IVll 

It is worth to noting that the response equation (|116p 
can be expressed in an explicit form for the Dirac model 
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(|108p . According to Eq. ([64)) . the momentum-space 
second Chern number of the (4 + l)-d Dirac model 
^0^) = Sa'^a(^)-'^° equal to the winding number of 
d(k) on the unit sphere 5^. Correspondingly, the Hamil- 
tonian of the (2 + l)-d Dirac model (|108p with constant 
9 and cp has the form /i(k, 6*, (^) = o?a(k, 0, iy9)r'', so 
the correlation fmiction G2 defined in Eq. (jllOp can be 
obtained as 

G2{0, ^) = -J dke |d(k,0,^)|5 • 

Thus the curvature of effective gauge vector potential fi^ 
is expressed as 



An 



|d(k,0,^)|5 



(117) 



Now consider a slightly different version of lattice Dirac 
model given by 

/i(k, n) = sin /c2;r^ + sin fcj,r^ + (cos fca: + cos /cj, — 2) 
+ni ^ naE", 

a=0,3,4 

in which m > and n = (no, fis, ^4) is a 3D unit vector. 
For such a model the d vector can be decomposed as 

d(k,0,^) =do(k)+ ( 
mn 

with do(k) = (sin /cj:, sin fcj,, 0, 0, cos fca: -l-cosfcj, — 2). In 
the limit m ^ 2, the Hamiltonian has the continuum 
limit h{k, h) ~ J2a=i.2 ^a'^"' + 1^6=0,3,4 mfibV'', which is 
the continuum 4x4 Dirac model with three possible mass 
terms. In this limit the integral over k in Eq. (|117p can be 
explicitly carried out, leading to the following expression: 



d^^u - dyVl^ = -n • S^n x 9,,n, 



which is the skyrmion density of the unit vector n. Com- 
bined with Eq. (fTT6)) we obtain the response equation 
for the Dirac model in the continuum limit: 



on 



(118) 



Eq. (|118p describes a topological response where the 
charge density and current are equal to the skyrmion 
density and current, respectively. Such an equation 
can be considered as a (2 + l)-d version of Goldstone- 
Wilczek formula (|3ip . which has been studied extensively 
in the literaturei^i^i^iSliSa. Thus, through the discus- 
sion above we have shown that the topological response 
formula( |118p of (2 -I- l)-d Dirac fermions is a special ex- 
ample of the generic response equation (|116p . 

To understand the physics described by Eq. (|116p . 
consider the lattice Dirac model in Eq. (|108p with an 
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FIG. 17: (a) Illustration of the charge pumping effect in the 
(2 -I- l)-d lattice Dirac model with spatial dependent 6 — 6{jj) 
and time-dependent tp — tp[t). This effect is the dimensional 
reduction of the domain-wall QH effect in the (3-1- l)-d system 
as shown in Fig. (b), in which dttp plays the role of electric 
field Ez- (c) Illustration of the half charge trapped at the 
crossing point of the 6 and ip domain walls. This effect is 
the dimensional reduction of the charge trapped on the do- 
main wall by a magnetic field, as shown in Fig. (d). ip{x) 
corresponds to the gauge vector potential as shown by red 
arrows, the curvature of which leads to a magnetic field By in 
the yellow region with net flux 2n. Due to the half quantized 
Hall conductance of the 9 domain wall, such a magnetic flux 
induces a half charge. 



adiabatic time-evolution of Lp{t) = 27rt/T, and a spatial 
domain wall configuration of 6{x): 



1 + tanh 



as shown in Fig. [T7] (a). According to Eq. (|116p . the 
charge current along the domain wall is given by 



27r 



{dyftt 



dtfly) 



When the parameter evolves adiabatically from to 27r, 
the net charge flowing across the line a; = is given by 

/ dth = / dtdyidynt-dtny)/2TT = de " difide^ - 

d^^ilg) /2Tr, which is the integration of the second Chern 
form Tr ^ f ab fc d\ /327r^ over the parameter range 

9 e [0,7r], </? g [0, 27r], where A,B,C,D = ka:,ky,9,(p. 
According to the discussion in the (3-|-l)-d case, we know 
that a magneto-electric polarization P3 can be defined as 

Psid) ^ J d^kdipKf = - j dipVL^/2iT, 

which implies / dtl^ = — d9dgP3{9). For —4c < m < 
—2c, we have dPz{9) — 1/2, corresponding to the 
pumped charge AQ = / dtlx = 1/2. In summary, the 
physical consequence of the topological response equa- 
tion (|116p is to induce a topological charge pumping effect 



30 



during the adiabatic evolution of tp, in which the charge 
pumped in one period is proportional to the spatial gra- 
dient of the magneto-electric polarization P3 defined in 
(f) space. Specifically, a charge e/2 is pumped 
along each A6 = it domain wall of when (/? evolves from 
to 27r. 

Such a charge pumping effect can also be viewed as 
the dimensional reduction of the half QH effect on the 
9 domain wall of (3 + l)-d lattice Dirac model, which 
has been studied in Sec. IIVBI This dimensional reduc- 
tion procedure is in exact analogy with the usual reduc- 
tion from the (2 -I- l)-d QH effect to (f -I- l)-d quantized 
pumping effect studied in Sec. Ill CI Similar to the lat- 
ter case, a fractional charge effect can also be proposed 
in (2 -I- l)-d according to Eq. (jll6p . To show this ef- 
fect, one can consider the same domain wall as shown 
above, and a 27r domain wall of ip along the j/-direction 
(/3(x) = TT [1-f tanh(a;/^)], as shown in Fig. [T7] (b). The 
charge density is given by p = {dx^y — dyVLx) /27r. By 
integrating over the x-direction we obtain J VLxdx — 
J fl^dip — — 27rF3(0), such that pi£, = J dxp = dyP-i{9) 
and / dypiY) — 1/2. Thus, a half charge is localized at 
the crossing of 9 and Lp domain walls^. Such a frac- 
tional charge effect can also be understood through the 
dimensional reduction from (3-1- l)-d. The spatial depen- 
dence of ipix) corresponds to the spatial dependence of 
kz — Az {x) , which describes a magnetic field perpendicu- 
lar to the 2D domain wall in (3-|-l)-d system. When tf{x) 
has a 27r domain wall, the net flux of the corresponding 
magnetic field is 27r, which thus induces a half charge as 
shown in Fig. [TTIfb). 

In summary, we have studied the physical conse- 
quences of the topological effective action (IllSp in a spa- 
tially and/or temporally inhomogeneous insulator. In the 
rest of this section we will show how to define a Z2 topo- 
logical invariant in (2 + l)-d TRI insulators and study 
the physical properties of the Z2 nontrivial phase — QSH 
phase — by applying the effective theory (|115p . 



B. Z2 topological classification of TRI insulators 

In Sec. Ill Dl and III El we have shown how a Z2 classifi- 
cation of particle-hole invariant insulators can be defined 
in both (1-1- l)-d and (0 -I- l)-d through dimensional re- 
duction from the (2-|- l)-d QH effect. The second-Chern- 
class analogy of the (2 -I- l)-d QH effect is the (4 -I- l)-d 
QH effecti^ described by the Chern-Simons theory ([5^ . 
which then leads to the Z2 classification of TRI insula- 
tors in (3 -I- l)-d, as shown in Sec. IIV CI Following this 
line of reasoning, it is straightforward to see that a Z2 
classification can be defined for (2 -I- l)-d TRI insulators, 
as an analog of (0 -I- l)-d particle-hole symmetric insula- 
tors. In this subsection we will sketch the demonstration 
of such a topological classification without going into de- 
tail since the derivation here is exactly parallel to that in 
Sec. iHEl 

First of all, for two TRI (2 + l)-d insulators 



/ii(k), /i2(k) an adiabatic interpolation h(^,9) can be 
defined, satisfying 

/i(k,0) hi, /i(k,7r) = h2 
T^h{-k,-9)T = h^{k,9), (119) 

Since h(k, 9) corresponds to the Hamiltonian of a (3 + 1)- 
d TRI insulator, a Z2 topological quantity A^3[/i(k, 9)] — 
±1 can be defined as shown in Sec. IIV CI The key point 
to defining a Z2 invariant for the (2 -I- l)-d Hamiltonians 
hi, /i2 is to demonstrate the independence of N-i[h{k, 9)] 
on the choice of /i(k, 0). Consider two different pa- 
rameterizations h{k,9) and h'{k,9). An interpolation 
(7(k, 9, ip) can be defined between them which satisfies 

g(k,0,O) = Mk,^), g{K9,^)^h\li,9) 
5(k,0, (^) = /ii(k), g(k,7r,(^) = /i2(k) 
g^{K9,p) = T^g{-k,-9,~p)T. 

.g(k, 9, ip) corresponds to a (4 + l)-d insulator Hamil- 
tonian, for which a second Chern number C2 [g] is de- 
fined. By its definition, the "second Chern parity" N3 
of /i(k,6') and h'{\i,9) satisfies TVs [/i] TVs [/i'] = (-l)^^^. 
At the same time, g(k, 9, p) can also be considered 
as an interpolation between = and 9 = t: sys- 
tems, i.e., between g(]s.,0,p) = /ii(k) and g(k,TT,p) = 
/i2(k). Since /ii^2(k) are both independent of p, the p- 
component of the Berry's phase gauge field vanishes for 
g(k,0,p) and g{k,TT,p). Consequently, it can be shown 
that {-lf2[9] = N3[g{k,0,p)] = m[g{k,Tr,p)] = 1, 
so that N3[h]N3[h'] — 1 for any two interpolations 
h and h' . Thus, we have shown that the Z2 quan- 
tity 7V2[/ii(k),/i2(k)] = N3[h{k,9)] only depends on the 
(2 + l)-d Hamiltonians hi and /12. By defining a constant 
Hamiltonian /io(k) = ho as reference, all (2 -|- l)-d TRI 
insulators are classified by the value of N2[hQ, h(k)]. An 
insulator with N2[hQ,h] = — 1 cannot be adiabatically 
deformed to the trivial Hamiltonian Hq without breaking 
time-reversal symmetry. 

In the next subsection, the physical properties of the 
Z2 non-trivial insulator defined here will be studied. We 
will see that the Z2 non-trivial insulator defined here 
has non-trivial edge dynamics, and corresponds to the 
QSH insulator studied in the literatur o^^'^^i^°'^^ . Com- 
pared to the former definition of the Z2 topological 
classificationiSiSi2^, our definition has the advantage of 
providing a direct relationship between the topological 
quantum number and the physical response properties of 
the system. 

C. Physical properties of the Z2 nontrivial 
insulators 

Similar to the (3 -I- l)-d case, the topological proper- 
ties of a Z2 nontrivial insulator lead to non-trivial edge 
state dynamics described by the effective theory (|115p . or 
equivalently, the response equation pi6p . The edge of a 



^2 nontrivial insulator is equivalent to a domain wall be- 
tween a nontrivial insulator and a trivial insulator (since 
the vacuum can be considered as a trivial insulator with 
a large gap). Thus, in the following we will focus on the 
domain wall between a nontrivial system with Hamilto- 
nian /ii(k) and a trivial system with Hamiltonian ho. 

As discussed in the last subsection, an interpolation 
h{h, 9) can be defined between Hq and hi satisfying 
h{k,0) = ho, h{k,n) = /ii(k) and rt/i(-k, -6i)r = 
/i-^(k, 0). Since hi is nontrivial, ft,(k, 6) has to break time- 
reversal symmetry for general 9 to adiabatically connect 
hi to Hq. Making use of /i(k, 0), two different interfaces 
between hi and ho can be defined. Consider a spatially 
dependent 9 given by 



(a) 



(b) 



0{x,y) 



1 — tanh 



Then the spatially dependent Hamiltonian ft,(k, 9{y)) and 
/i(k, —6{y)) both describe a spatial domain wall between 
hi (for y <C — 5) and ho (for y 3> ^). The only dif- 
ference between these two Hamiltonians are the time- 
reversal symmetry breaking terms around the interface. 
Now consider a more complicated interface, with 



/i(k,x) 



h(k,9iy)), x<0 
/i(k,-%)), a;>0 



as shown in Fig. [TS] (a). In such a system, the time- 
reversal symmetry on the interface is broken in oppo- 
site ways for a; > and x < 0, in the sense that 
h^{k,{x,y)) = Tt/i(-k, (-a;,y))r. Now we study the 
charge localized around the point a; = 0, y = 0. For a 
loop C enclosing this point as shown in Fig. [18] (a), the 
charge in the region A enclosed by C is given by Eq. 
(fTTe]) as 



2tt 



(fx {dx^y - dyVlx) 



2-K 
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When the size of the loop is large enough compared to the 
boundary width ^, such a loop integration is equivalent 
to an integration over 9 from to 27r, which leads to Q = 
§ VtgdB /2tt — F3[/i(k, 9)]. According to the definition of a 
Z2 non-trivial insulator in the last section, P3[/i(k, 0)] = 
1/2 mod 1 for any interpolation hik.,9) between ho and 
hi. Consequently, the charge confined on the domain 
wall \s Q = {n -\- l/2)e with n an integer depending on 
the details of the interface^. 

To summarize, a time-reversal symmetry breaking 
term can be applied at the interface of trivial and non- 
trivial insulators. For a given interface described by 
Hamiltonian /i(k, y), its time-reversal partner, h'(k,y) = 
r/i"^(— k, y)T'' , describes a different connecting condition 
at the interface. If the ID interface is described by h in 
one region and by h' in another region, then the domain 
wall between these two regions will trap a half-charge as 
a consequence of the non-trivial topology. To understand 
such a domain wall better, we can consider the case with 
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FIG. 18: (a) Illustration of an interface between the vacuum 
(VAC) with Hamiltonian /io(k) and a QSH insulator (QSH) 
with Hamiltonian /ii(k). An interpolation h(k,8),0 £ [0, 27r] 
can be defined between ho and hi. On the left (right) half of 
the interface marked by red (blue), 9 has a domain wall from 
to — TT (tt). a fractional charge Q = e/2 is trapped on the 
domain wall between the red and blue interfaces, which can 
be calculated by an integration of along the loop L (see 
text), (b) Physical realization of the domain wall between 
two interfaces in (a) by an anti-phase domain wall of magnetic 
field or magnetization. The red curve shows schematically the 
charge density distribution. 



a magnetic field as the time-reversal symmetry breaking 
term on the interface. When the magnetic field has an 
anti-phase domain wall as shown in Fig. [TH] (b) , a half- 
charge must be trapped on the domain wall. 

Moreover, one can also obtain the distribution of ID 
charge density and current density on the interface by 
integrating the Eq. (|116p only along the y-direction: 
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Pidix) = — y dy {dxVly - dy^lx) 
iid{x) = J dy {dyVtt - dtVty) 



in which L is a cut-off in the y-direction, satisfying L 3> ^ 
so that the contribution to piD and jid from the region 
|y| > L is negligible. According to the definitions (|112p 
and pi3p of the effective gauge vector potential il^, we 
know that J7^(a;,y,t) — s- for a point deep in the QSH 
or VAC region, i.e., when \y\ — > 00. Thus the expression 
of ID density and current can be simplified to 

Pm{x,t) = dj;P3{x,t), jiB{x,t) = -dtP3{x,t) 

(120) 

with ^3(0:, t) — j'^j^ dyVly{x, y, t)/2-K the magneto-electric 
polarization defined in (kx,ky,y) space. Eq. (|120p is 
exactly the Goldstone-Wilczek formulcsi^ describing the 
charge fractionalization effect in the (1-1- l)-d Dirac 
model, and 2ttP3{x, t) plays the role of the phase angle of 
the Dirac mass term. When the interfaces on the left and 
right sides of the domain wall are related by time-reversal 
symmetry, the change of 27rP3 through the domain wall 
must be (2n -f l)7r, which gives the half charge on the 
domain wall. 

Such a relation between the interface and the (1 -f l)-d 
Dirac model can be understood more intuitively in terms 



FIG. 19: Illustration of the dimensional reduction from the 
surface of a (3 + l)-d topological insulator (a) to the edge of 
a (2 + l)-d QSH insulator (b). The red and blue points in 
(a) are the positions of gapless (2 + l)-d Dirac cones in the 
surface BZ. The blue line sd, ky — tt defines the edge theory 
of a QSH insulator, as shown in Fig. (b). 

of the edge effective theory^. Similar to the relation be- 
tween the edge theory of a (4 + l)-d topological insulator 
and the reduced edge theory of a (3+ l)-d Z2 non-trivial 
insulator, we can obtain the edge theory of a (2 + l)-d 
Z2 nontrivial insulator from that of a (3 + l)-d nontrivial 
insulator. The interpolation ft,(k, 6) between hi and ho 
can be viewed as the Hamiltonian of a (3-|-l)-d TRI insu- 
lator, in which 9 plays the role of k^. Consider a specific 
point of the (2 -I- l)-d boundary of the (3 -I- l)-d system, 
say y — 0. According to the discussion in Sec. IIV C| an 
odd number of (2 -I- l)-d Dirac fermions are propagat- 
ing on the boundary of the (3 -I- 1) — d system. Due to 
time-reversal symmetry there must be an odd number of 
Dirac cones on the four time-reversal symmetric points, 
as shown in Fig. [T^] (b). For a slice at y = 0, the wavevec- 
tor of the surface state is {kx,0). Consequently, when 9 
is considered to be a parameter, the surface energy spec- 
trum E{kx,9) for a given 9 describes the dispersion of 
(1 + l)-d edge states of (2 + l)-d insulators. Specifically, 
9 = corresponds to the vacuum Hamiltonian /ig, which 
cannot support any non-trivial edge states. Thus, the 
Dirac cones at the time-reversal symmetric points can 
only appear at {kx,9) = (0, tt) and (tt, tt). To have the 
minimal odd number of Dirac cones in the (kx,9) BZ, 
there must be one gapless Dirac cone at (0,7r) or (tTjTt), 
but not both, as shown in Fig. [121 In other words, the 
edge effective theory of a (2 + l)-d nontrivial insulator is 
given by a gapless Dirac theory 

/dk 
—V ak^jliJka (121) 
a=±l 

where a = ±1 means left and right movers, respectively. 
This edge theory agrees with the former descriptions of 
QSH edge states^^'^^'^^ and shows the equivalence of the 
Z2 nontrivial insulator defined in this section and the 
QSH insulator. 

Once the edge theory is obtained, it is easy to under- 
stand the charge fractionalization proposed above. Due 
to Kramers's degeneracy, any TRI perturbation cannot 



FIG. 20: Illustration of quantized charge pumping on the 
boundary of a QSH insulator induced by a rotating magnetic 
field. During the time period t G [0, T] the magnetic field is 
rotated from B to —B, and a half charge Q = J dt J_^^ dyjx = 
e/2 is pumped along x direction. L is taken as a cut-off with 
jx vanishingly small for \y\ > L. 

open a gap at the edge. Only when a magnetic field, or 
other time-reversal symmetry breaking term is applied, 
can a mass term rrixcrx + friy'^v be generated in the edge 
theory. Time-reversal symmetry also guarantees that the 
mass induced by opposite magnetic fields is exactly op- 
posite. This implies that an anti-phase domain wall of 
the magnetic-field corresponds to a sign-change of the 
mass of the Dirac fermion. Thus, the edge state theory 
is described by the well-known Jackiw-Rebbi model^ or 
equivalently, the Su-Schrieffer-Heeger model^. The study 
of the fractional charge in the edge theory approach and 
its experimental consequences have been presented in 
Ref. [6^. Thus, we have seen that the effective theory 
we obtained from dimensional reduction agrees with the 
edge theory analysis, just hke in the (3 -I- l)-d case. The 
effective theory correctly describes the half-charge asso- 
ciated with a magnetic domain wall, which is a direct 
physical manifestation of the Z2 quantum number. 

A quantized charge pumping effect always accompa- 
nies a fractional charge effect and can be realized when a 
time- dependent T-breaking field is applied at the edge. If 
the system is described by a time-dependent Hamiltonian 
/i(k,x, t) which satisfies /i(k, x, t = 0) = h{'k,9{y)) and 
h{k,x,t = T) = h{k,~9{y)) with 9{y) the domain waU 
configuration discussed above, then the charge pumped 
through the interface during the time t G [0, T] is given 
by 

Qpump = ^ dtjiB = - [PsiT) - P3(0)) = - + • 

In the example of an applied magnetic field, such a pump- 
ing process implies that a half charge is pumped when a 
magnetic field rotates adiabatically from B to — B, as 
shown in Fig. 1201 The experimental proposal of such a 
charge pumping effect is also discussed in Ref. 69. 

Besides providing a quantized response property of 
the QSH insulators, the fractional charge and charge 
pumping effects proposed here are a (1 -I- l)-d version 
of electro-magnetic duality. In (3 -I- l)-d, the electro- 
magnetic duality gives rise to the Witten effect^^, where 
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a magnetic monopole carries a charge Q/2tt and becomes 
a "dyon" when a topological 8-term is introduced in the 
Lagrangia n^^i^° . Such an effect can occur in a (3 + l)-d 
topological insulator where the charge of the dyon is 1/2, 
as studied in Sec. IIVBI In comparison, the magnetic do- 
main wall on the boundary of (2 + l)-d QSH insulator 
can be considered as a topological point defect of mag- 
netic field in (1 -I- l)-d, which also carries a half-charge. 
In this sense, if we consider the magnetic domain wall 
as a dynamical degree of freedom of the system, [e.g., 
when the magnetic domain wall is generated by a ferro- 
magnetic stripe on top of the 2D QSH system) it can be 
considered as the (1 -I- l)-d manifestation of dyons. 

Interestingly, such an analogy can also be generalized 
to (2-t-l)-d, where the topological defect of magnetic field 
is a flux tube. Recently it has been shown that a tt fiux 
tube threaded into a QSH insulator carries either charge 
±e, spin or charge 0, spin 1/2, where the spin (1/2) is 
generically defined as a Kramers' singlet (doublet) under 
time-reversal symmetry i^i^ In other words, the tt fiux 
tube becomes a dyon-like object and realizes spin-charge 
separation in (2 -I- l)-d. Such a spin-charge separation 
phenomenon also provides an alternative definition of the 
Z2 topological insulators in (2 -I- l)-d^. 

In summary, we have shown how a Z2 classification of 
(2 -|- l)-d TRI insulators is obtained, and how the physi- 
cal properties of the Z2 non-trivial insulator are described 
by the effective theory derived from dimensional reduc- 
tion. Together with the discussion of (3-|-l)-d topological 
insulators in Sec. IIVI we have seen that the nontriv- 
ial topology and its consequences in both (3 + l)-d and 
(2 -|- l)-d TRI systems have their origin in the nontriv- 
ial second Chern number in (4 -f l)-d. The dimensional 
reduction series (4 + l)-d ^ (3 + l)-d ^ (2 + l)-d is 
in exact analogy of the lower dimensional one (2 -f 1)- 
d ^ (1 + l)-d ^ (0 + l)-d. In next section, we will 
develop the unified framework of dimensional reduction 
in generic dimensions, which contains the two series as 
simplest examples. 



VI. UNIFIED THEORY OF TOPOLOGICAL 
INSULATORS 

A. Phase space Chern-Simons theories 

Up to now, we have systematically studied several re- 
lated topological phenomena, including the (2 -I- l)-d QH 
insulator with nontrivial first Chern number, the (4+l)-d 
topological insulator with nontrivial second Chern num- 
ber, and their dimensional reductions. Comparing Sec. 
ITTl with Sees. Illlirvl one can easily see the exact anal- 
ogy between the two series of topological insulators: the 
(2 + l)-d and (4 + l)-d fundamental topological insula- 
tors are characterized by an integer — the first and second 
Chern number, respectively. Under a discrete symmetry 
(particle-hole symmetry for the (2-|-l)-d family and time- 
reversal symmetry for the (4-f l)-d family), a Z2 topologi- 



cal classification can be defined for the lower dimensional 
descendent systems; the physical properties of which can 
be described by effective theories obtained from the di- 
mensional reduction procedure. The main facts about 
these topological phenomena are summarized in Table [H 
In this section, we will show that the effective theories 
for all these systems share a universal form when written 
in phase space. 

As a simple example, we first consider the effective the- 
ory ([9|) of the QH effect. Expanding the expression of the 
first Chern number explicitly, Eq. ([9]) can be expressed 
as 

S=^J ^e''Tr[a,aj] J dtd^xe^^^PA^d^Ap (122) 

in which i,j are indices 1,2 in momentum space and 
fi,v,p = 0,1,2 are space-time indices. Here and below, 
and a.i stand for the external electromagnetic gauge 
field in real space and Berry's phase gauge field in mo- 
mentum space, respectively. The trace is carried over all 
occupied energy levels. If we define the phase-space coor- 
dinate as q = {t,x,y,kx,ky) and the phase-space gauge 
potentials A = (^0,^1,^2,0,0), a = (0, 0, 0, oi, a2), 
then the action above is equivalent to the following sec- 
ond Chern-Simons term: 

^2+1 = ^ / d^'qe^^'^^'^AAdBAcTr [doaE] (123) 

where all capital roman indices e.g. A, B,C . . . run over 
the appropriate phase space coordinates. Since in this 
system A and a are always orthogonal to each other, 
such a reformulation seems trivial. However, it turns out 
to be helpful when considering the dimensional reduction 
procedure. As discussed earlier, dimensional reduction of 
the (2 + l)-d system to (1 + l)-d is defined by replacing, 
say, ky -|- by a parameter 9{x,t), which in general is 
space-time dependent. Four changes are induced by this 
substitution: 

1. Ay is replaced by 9{x,t) — 9q with 0o playing the 
role of ky 

2. d/dky is replaced by d/dO 

3. The Berry's phase gauge field a'^^ = 

—i (a; k| dky |/3; k) is replaced by a^^ = 
-i {a;ka:,e\de \l3;kx,0) 

4. The integrations /dy and J dky/2TT are removed 
from the effective action. 

By making these substitutions, the effective action of the 
(1 -I- l)-d system can be obtained. To help understand the 
general dimensional reduction procedure, we show the 
derivation of the (1-1- l)-d effective theory explicitly. For 
simplicity, one can start from Eq. (|122p . Note that the 
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Dimension 


Topological quantum 
number 


Effective theory 


Symmetry 
requirement 


Physical properties 


Family 1 


2+1 


1"* Chern number Ci £ Z 




Not required 


QH effect; axial anomaly 
on the boundary 


1+1 


Chern parity 




C 


Half charge on the 
boundary 


0+1 


1""' Chern parity 
A^o = Ni\ho,h] e Z-i 


C = Tr[ao] 


C 


Not applicable 


Family 2 


4+1 


T"^ Chern number C2 G Z 




Not required 


4DQH effect; chiral 
anomaly on the boundary 


3+1 


2""^ Chern parity 

% = (-l)C2[h0,hl £ ^2 


£ = j:;;P3t^'"'''d^A^d„Ar 


T 


Half QH effect on the 
boundary, TME effect 


2+1 


2"^* Chern parity 
A^2 = Nz\hn,h\ G ^2 




T 


QSH effect; half charge at 
anti-phase domain wall on 
the boundary 



TABLE I: Summary on the properties of the topological insulators. In the effective Lagrangians, the indices q,/3.. = 0,1, 
= 0,1,2, ^i,v.. = 0,1,2,3 and a,h.. = 0,1,. .,4. The vector ai (a^, et al.) stands for the gauge vector of the external 
electro- magnetic field, and Ai that for the Berry phase gauge field. C and T stands for particle-hole symmetry and time-reversal 
symmetry, respectively. The quantities Pi, P3, are defined by Eqs. (f24)) . ([79j and (|112p . See text of Sec. IVII for explanations 
of the effective actions. 



gauge field depends only on the (1 + l)-d coordinates 
(t, x) after dimensional reduction. Consequently, in the 
Chern-Simons form e^'^'^ A^^d^Ar the terms containing dy 
are identically zero, so that 

e^'^A^d.Ar = 2 {Atd^Ay - A^dtAy) - 2 {Atd^9 - A,dt9) 

after an integration by parts. The effective action (|122p 
after making all four substitutions above is expressed as 

Si+i = ^ J dkj:TT [dk^ag - dettkj 

dtdx {Atd,J ~ A^dt0) . 

With a smooth space-time dependent 9{x,t) field, the 
eigenstates \a;kx, 6) can be considered as space-time de- 
pendent "local" eigenstates, whose space-time depen- 
dence originates only from that of 9{x,t). In this way, 
each state \a;kx,9{x,t)) = \a;kx,x,t) is defined in the 
full phase space {t,x,kx), and the Berry's phase gauge 
field can gain real-space components defined as a"^ — 
—i{akx,x,t\d^\(3,kx,x,t) — aed^9, in which ^ — t,x. 
Similarly, the space-time derivative of a^^ is given by 
dfiO-k^ = deak^dfid- By making use of these observations, 
the effective action can be simplified to 

Si+i = ^ J dkx J dtdx {AtTr [dk^Qx - dxttkj 

-AxTT[dk^at~dtakJ}. (124) 

By generalizing the definition of the gauge vector poten- 
tial {At, Ax) to the phase space vector {At,Ax,0), the 
equation above can be expressed as the mixed Chern- 
Simons term in the phase space: 



which describes an inhomogeneous (1 + l)-d insulator. 
Note that momentum derivatives acting on {At, A^) van- 
ish. This effective action agrees our discussion in Sec. 
Ill CI as one can see by taking ^^^j- and comparing the 
resulting response equations with Eq. ([28|l . 

It should be noted that such a phase space formalism 
is only applicable when the space-time variation of 9 is 
smooth and can be approximated by a constant in the 
neighborhood of a space-time point (i, x, y). More quan- 
titatively, the characteristic frequency lu and wavevector 
k of 9{x,y,t) should satisfy 



hu, hvk <^ E, 



9' 



(126) 



where u is a typical velocity scale of the system. For 
example, in the lattice Dirac model p2|) . v is the speed 
of light (which is normalized to be 1 in Eq. p2p ). Under 
condition (|126p the space-time variation of the 9 field 
does not generate excitations across the gap, and the 
system can be viewed locally as a band insulator with a 
"local Hamiltonian" h{kx',x,t). 

Carrying out such a procedure once more to the ac- 
tion (|125p one can obtain the (0 + l)-d action. We 
will show the derivation explicitly. First one must take 
ax — 0, dxttk^ = in Eq. p24p since nothing can de- 
pend on the spatial x-coordinate after it is dimensionally 
reduced. Next we replace Ax by the parameter (f)[t) — (po, 
which leads to 



So+i = - I dt{(j>- 0o) Tr [d^at - dta^] 
dt {dt(j>d^ + dt) . 



d'qe^^'^AATvidBac] 



(125) 



where the integration J dxdkx /2t: has been removed, and 
at, Gk^ in (l + l)-d are relabelled as dt, d^ for later conve- 
nience and to obtain the second equality, an integration 
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by parts is carried out. It should be noted that 9^ comes 
from d/dky, which actuaUy means d/d(j)Q since it is (j)Q 
that is replacing ky. Compared to the (2 + l)-d ^ (1 + 1)- 
d dimensional reduction, the difference here is that the 
wavefunctions are, in general, already time-dependent in 
(1 + l)-d. This comes from the dimensional reduction 
from (2 + l)-d. Consequently, the Berry's phase gauge 
potential in the (0 + l)-d system consists of two terms as 
shown below: 



d_ 

dt \d4>y , 



1/3; t, 0) 



in which (d/dt)^ means to take the f-derivative while 
keeping constant. Both of these terms are necessary 
for the correct topological response and similar terms (in- 
cluding spatially dependent ones) will be present in all 
higher dimensions when more than one reduction is car- 
ried out. Combining the two equations above we finally 
obtain 



dt Tr [at\ 



(127) 



which has the form of a "zero-th" Chern-Simons term, 
and describes the "(0 -I- l)-d insulator" discussed in Sec. 
IIIE( i.e., a single-site fermion system. In this case the 
only gauge invariant quantity remaining is the Berry's 
phase the single-site system obtains during a closed path 
of adiabatic evolution. 

For the second family of topological insulators we dis- 
cussed, the effective theory ([52]) can be expressed in the 
following phase space form: 



1^ / cPqe^^^^^^^^^'AAdsAcdoA, 
IZtt^ J 



1927 

•Tr [DpaGDHai] (128) 

in which the covariant derivative Db =83+ ias is in- 
troduced for the non-abelian Berry's phase gauge field. 
The dimensional reduction to (3-1- l)-d can be performed 
similarly to the (2 -I- l)-d case. Denoting the 4-th spatial 
dimension as w, which is the dimension to be reduced, 
then any term with dw vanishes, and so does aw Con- 
sequently, the only non-vanishing terms in the effective 
action (|128p are those with Ayj, which now is replaced 
by the parameter 0{x, y, z, t). On the other hand, one of 
F,G,H,I in the form DpacDuaj has to be k^, which 
is now replaced by 9. In summary the theory can be 
rewritten as 

53+1 = ^ ! qe'^'-^^ e'^^ A^dMrO 
Odtt'' J 

•Tr [DeaiDjak + cycL] 

where fi,!^.. — 0,1,2,3 and i,j,k — 1,2,3 are space-time 
and momentum indices of (3-1- l)-d system, and cycl. de- 
notes the other three terms obtained from cyclicly per- 
muting 9 and i,j,k. The integration J dk^dw /2t: has 
been removed, and a prefactor 3 appears due to the fact 



that there are three A as in the effective action (jl28l) . In 
the same way as in (2 -I- l)-d to (1 -|- l)-d case, df^ODe 
can be replaced by I?^ , so that the effective theory of the 
(3 + l)-d insulator is finally obtained: 



53+1 



AB..G 



327r3 



AaObAcTt [DDaEDpaG] 



(129) 



According to the definition of P3 in Eq. ([79]) we know 
that 



which shows the equivalence of the action (|129p to the 
action ((80|) we derived earlier. 

Now from the two examples of (2 -|- l)-d — > (1 -I- l)-d 
^ (0 l)-d and (4 + l)-d ^ (3 l)-d, one can easily 
obtain the general rule of dimensional reduction to the 
phase-space Chern-Simons theories. For a (d+l) dimen- 
sional system, the phase space dimension is 2d -I- 1, and 
the dimensional reduction of corresponding phase-space 
Chern-Simons theory is defined as 

1. Remove a term OaAb from the action and cor- 
respondingly replace the {2d -f l)-dimensional 
anti-symmetric tensor by the one in {2d — 1) 
dimensions. 



2. Remove the integration J dxddk^/ 2Tr from the 
action when x^, kd are the spatial and momentum 
indices to be reduced. 



3. Multiply the action by a factor n when the power 
of external gauge field Aa in the original action is 
n. 

Following these rules, the effective action for (2 + l)-d 
TRI insulator can be easily obtained by one more step of 
dimensional reduction from Eq. (|129p : 



2+1 



d^q 
8^' 



ABODE 



AATi[DBacDDaE] (130) 



where the coefficient is determined by 87r^ = 327r'^/(2 • 
27r). By considering the space-time and momentum in- 
dices separately, one can easily confirm that Eq. (|130p is 
equivalent to Eq. (jlisp we obtained earlier. 

In summary, we have shown that all the known topo- 
logical insulators are described by a Chern-Simons effec- 
tive theory in phase space, and the topological theories 
in different dimensions can be related by the dimensional 
reduction procedure. It is straightforward to generalize 
this formalism to arbitrary dimensions. As shown in Ap- 
pendix [B] the relation between the non-linear response 
function (j53[) and the corresponding Chern number in 
momentum space can be generically proven for any odd 
space-time dimension. The effective theory of such a 
{2n + l)-d topological insulator is given by^ 
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T 




generation 

FIG. 21: The family tree of topological insulators. The black points on the "trunk" (i.e., 0-th generation) stand for the 
fundamental topological insulators in odd space-time dimensions characterized by a nontrivial Chern number in momentum 
space. The blue and red markers show the descendants of the (4n — l)-d and (4n + l)-d insulators, respectively. Physical 
effects associated with some of the prominent topological insulators are indicated in parenthesis. A Z2 classification is defined 
for each blue circle (square) under (pseudo) particle-hole symmetry C (C), and for each red circle (square) under (pseudo) 
time-reversal symmetry T (T). The definitions of C, C, T, T are given in Sec. I VI Bl For all the physically realizable systems 
with spatial dimensions d < 3, the names of topological insulators are labelled, where CIn (Tin) stands for a particle-hole 
symmetric (TRI) topological insulator in n + 1 dimension. The black circles and squares stand for other topological phenomena 
obtained from dimensional reduction, which are also described by the phase space Chern-Simons theories but do not correspond 
to Z2 topological insulators. The phase space Chern-Simons theory CS2„+i (as defined in Eq. (|134|l ') corresponding to each 
topological phenomenon is also specified on the figure. 



^" [ fj2n+l tiitJ.2..tJ.2„ + l A r) A r) A 

with the n-th Chern number in momentum space defined as 

^" ^ n!2"(27r)" / ' ^''"'^^ [•^*i*2-^^3^*--'^»2„-i»2„] ■ (132) 
Thus the (4n -f l)-d phase space formula for this effective action can be written as 

FoUowing the general rules of dimensional reduction procedure discussed above, one can obtain the effective actions 
for lower-dimensional topological insulators as "descendants" of the topological theory p33p . From the examples 
discussed above it can be seen that the number of Berry's phase gauge vectors in the effective action remains 
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invariant during dimensional reduction, while the number of external gauge field insertions A^^ decreases by one at 
each step of dimensional reduction. After n + 1 reductions we obtain an effective action of an ((n — 1) + l)-d system 
that contains no Ai. Just like in the (0 + l)-d case, such an effective action does not result in any response equation 
of the system, but only describes the Berry's phase the system obtains during adiabatic evolution. Obviously the 
dimensional reduction cannot be carried out again on such an {{n— 1) + l)-d system. Thus, the (2n + l)-d topological 
insulator with a nontrivial n-th Chern number only has {n + 1) "descendants" under dimensional reduction. It is 
straightforward to show that the m-th descendant (1 < m < n + 1) of the (2n + l)-d topological insulator has the 
effective action 

, , (2n+l-m\ n .4n+l-2m„ 

•Tr [DA2,^ + 2-2m^A2r^ + 3-2m■■DA4,„-2m^A4,r^ + l-2m] 

= ClS4„_2m+l (134) 

in which CS* stands for the mixed Chern-Simons action in s phase-space dimensions with t powers of the external 
Aa field. Specifically, 08° is a pure non-Abelian Chern-Simons term of the Berry's phase gauge field a a, which 
cannot be reduced to a function of /ab = Daclb ~ Dbcla alone. Thus we have seen that a whole family of topological 
phenomena are described by phase-space Chern-Simons theories with different s, t values. In a given spatial dimension 
d, all possible topological phenomena in band insulators are given by the actions C&^^j^i with all possible values of 
t. It should be noted that the external gauge field Aa is only defined in real space, meaning that Aa or ObAa are 
both vanishing if A or i? is a momentum index A, i? = d + 2, d + 3, 2d + 1. Since in CSj^^j^ there are t Aa's and 
at least t — 1 partial derivative operators acting on A^'s, the Chern-Simons action CS2(i^i vanishes if 2t — 1 > d + 1. 
Consequently, there are in total [d/2] + 2 available Chern-Simons terms in the phase space of a (d -I- l)-dimensional 
system, which are 

CS*rf+i, i = 0,I,...,[d/2] + l. 

Here [d/2] denotes the maximal integer that does not exceed d/2. For example, in (2 + l)-d there are three available 
phase space Chern-Simons terms, two of which are CS5 describing a QH insulator, and CS5 describing a QSH insulator. 
The third one is given by 

CS° - j d''qe^^''''''Tr[aAdBacdDaE 

3^3 
+-aAaBacODaE + -aAaBacaoaE 
2 5 

which contains no Aa and thus does not describe any electromagnetic response properties of the system. The 
information contained in the effective action CS5 is the Berry's phase the system obtains during adiabatic evolution, 
just like the effective action of the (0 + l)-d system CSj = J dtao- We have grouped the phase-space Chern-Simons 
theories based on the parent theories and their descendants. The relationships are summarized in a "family-tree" 
in Fig. [5T] Similar to the generalization of the (2 -I- l)-d QH insulator to any odd space-time dimension, the Z2 
topological insulators we have studied can also be generalized to higher dimensions, which will be explained in the 
next subsection. 

Before moving on to that, we would like to point out an interesting mathematical fact about the phase-space 
Chern-Simons theories. For an (n -I- 1) dimensional system with A'^ occupied bands, a U{N) gauge vector potential 
can be defined in phase space as 

A\ = XAa + ua, 

with Aa being the external gauge potential and ua being the Berry phase gauge potential. The non-Abelian Chern- 
Simons term for Aa can be expressed as 

CS2„+i(A) ^ CS2„+i [A\] = J d2n+i^^AM....A.„+,Tr [A\^dA2A\^...dA2„Ai,+, + N.T.] . (135) 

I 

Here N.T. stands for the non-Abelian terms containing commutators of Aa, which can be determined by the 
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relation of the (2d+l)-d Chern-Simons form to the {{2d+ 
1) + l)-d Chern form. For more details, c.f. Section 
11.5 of Ref. [tiI. By expanding the parameterized action 
CS2r!+i(A) over A, the following equality can be obtained: 



CS2ri+l(A) — 



[n/2] + l 

E 

t=0 

1 d' 



A*CS2„_|.i 



CS2„+i - ^^^CS2„+i(A) 



(136) 



A=0 



This implies that all possible phase space Chern-Simons 
terms can be obtained from a single "generating func- 
tional" CS2n-i-i(A). We present l|135p as the unified the- 
ory of all topological insulators. 



B. Z2 topological insulator in generic dimensions 

For the descendants of the (2 + l)-d and (4 + l)-d in- 
sulators, we have defined a Z2 classification under the 
constraint of a discrete symmetry. For the descendants 
of the (2 + l)-d QH insulator, the Z2 classification is de- 
fined for particlc-hole symmetric insulators satisfying Eq. 



C''h{-k)C = -h^{k), C^C = C*C = 1. (137) 

The key point of this classification is to show that an in- 
terpolation between two particle-hole symmetric Hamil- 
tonians hi{k) and h2{k) forms a closed path when com- 
bined with its particlc-holc transformed partner. The 
Chern number enclosed in such a closed path always has 
a certain parity which does not depend on the choice of 
the path. In the same way, a Z2 classification of particle- 
hole symmetric insulators is also defined in (0 -I- l)-d. 
For the family of (4-|- l)-d insulators it is the same story 
except that the particle-hole symmetry is replaced by 
time-reversal symmetry 



Tt/i(-k)T = /i^(k), T^T = -T*T = I. 



(138) 



Following this one can easily generalize such Z2 classifi- 
cations to higher dimensions. To do that, one first needs 
to understand what is the difference between (2-1-1) and 
(4-1-1) dimensions that requires the choice of different 
discrete symmetries. The easiest way to see such a dif- 
ference is to study the transformation of the correspond- 
ing Chern-Simons theories under particle-hole symmetry 
(C) and time-reversal symmetry (T). Under particle- 
hole symmetry, the charge density and charge current 
both change sign. The vector potential does as well, as 
required by the invariance of the minimal coupling A^j^^. 
In the same way one can obtain the time-reversal prop- 
erty of A^, as summarized below: 



C:A,. 



Af_i , T : Af^ 




(139) 



In both cases of C and T, the momentum operator —id^ 
has the same transformation property as A^. Based on 
these facts the transformation properties of the Chern- 
Simons Lagrangian (jl3ip are 



n ■ qCS 

'-^ ■ '-'2n+l 
^ ■ '^2n+l 



n+1 qCS 



(-1)"+^S'. 



2n+l 



/ 1 yi cCS 
I ^) '-'2n+l 



(140) 



Thus, we see that S'^^^^j^ is T-even but C-odd, while 
S2n-i is T-odd but C -even. In other words, a (4n-|- l)-d 
topological insulator has to break particle-hole symme- 
try but can be time-reversal invariant, just like the case 
of (4 -I- l)-d; a [{An — 2) -I- l]-d topological insulator has 
to break time-reversal symmetry but can be particle-hole 
symmetric, just like the case of (2 + l)-d. Consequently, 
for the descendants of (4n-f l)-d ((4n— 1) — d) topological 
insulators, it is only possible to define Z2 topological clas- 
sifications by the dimensional reduction procedure under 
the constraint of T (C) symmetry. 

Naively, it seems that the procedure we introduced 
to define the Z2 classification by dimensional reduction 
could be applied recursively to all the descendants of a 
{2n+ l)-d topological insulator. However, this turns out 
to be incorrect. As an example, we can study the (1 + 1)- 
d TRI insulator as a descendant of the (2 + l)-d QSH 
insulator(not the (2 -I- l)-d QH insulator). In the dimen- 
sional reduction from (3 + l)-d to (2 -f l)-d discussed in 
Sec. IV Bl we define an interpolation ft,(k, 0) between two 
(2 + l)-d TRI Hamiltonians hi(k) and ft.2(k). When the 
interpolation h(k, 9) is required to satisfy time-reversal 
symmetry (Eq. (|119p ). it corresponds to the Hamilto- 
nian of a (3 -I- l)-d topological insulator, for which a Z2 
index N3[h{'k,9)] can be defined. In Sec. IV Bl we have 
shown that N3[h{'k,9)] does not depend on the choice of 
the interpolation h{'k,9), which thus provides a criteria 
on whether /ii(k) and /i2(k) are topologically equivalent. 
If we carry out the same procedure on (1 + l)-d TRI 
insulators, it seems that for two Hamiltonians hi{k) and 
/i2(fc) a Z2 topological classification can be defined in the 
same way. To see if this is true, one can again take the 
lattice Dirac model as an example. The single-particle 
Hamiltonian of (2 -I- l)-d 4x4 lattice Dirac model with 
time-reversal symmetry is written as 

ft.2D(k) — sin /cj; + sin fcj, 

-t-r° [rn + c(cosfca; -t- cosfcy)] (141) 

which is in the topological nontrivial phase for < to < 
2|c| or — 2|c| < TO < 0. By dimensional reduction, 
/i2D(k) can be considered as the interpolation between 
two (1 -|- l)-d TRI Hamiltonians hi{k) — h2D{k,0) and 
h2{k) — /i2D(^,7r). Thus if the Z2 classification pro- 
cedure applied to (1 + l)-d systems, hi{k) and h2{k) 
should be topologically distinct when /i2D(k) is in the 
nontrivial phase. In other words, it should be impossi- 
ble to define another interpolation ft.o(A:, 9) between hi{k) 
and /i2(fc), which satisfies T^ho{k,9)T = h^{-k,-9) = 
hQ{—k,9), V0. However, such an interpolation can in fact 
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be constructed as follows: 

hoik, 9) = sin fc + sin^ 61 

+r'' (m + ccosfc^; + ccos6') (142) 

in which F''^ = iV^T^ is even under time-reversal. 
The existence of two topologically distinct interpolations 
ho{k, 6) and h2D{k, 0) shows that it is not possible to de- 
fine a Z2 classification of (1 -f l)-d TRI insulators in the 
same way as in (2 -I- l)-d and (3 + l)-d. The main reason 
for the failure is that the proof in Sec. (jVBp requires 
the parameterized Hamiltonian manifold to be simply 
connected. In other words, an interpolation /i(k, 0, 1/3) 
can always be defined for two interpolations /i(k, 0) and 
/i'(k, 0). Similar arguments do not work in the classifi- 
cation of (1 -t- l)-d Hamiltonians because it may not be 
possible to adiabatically deform one path to the other. 
In the example of the lattice Dirac model, the paths /12D 
and /iQ cannot be adiabatically connected, because the 
combined path 

aike)^{ ^2D(fc,0), 0G[O,^] 
' ' \ ha{k,2TT-e), ee [tt,2tt] 

is a (2 -|- l)-d Hamiltonian that breaks time-reversal sym- 
metry and has a nontrivial first Chern number Ci = —1. 
Consequently, the path g cannot be contracted to a point, 
and the definition of a path-independent Z2 invariant 
fails. 

From this example we have seen that the definition 
of a Z2 topological classification for the descendants of 
(2ri-|- l)-d topological insulators fails when the dimension 
is reduced to ((2n — 3) -I- l)-d, since the lower Chern num- 
ber Cn-i is defined for each closed path of ((2n — 3) + l)-d 
Hamiltonians, thus obstructing the adiabatic connection 
between two different paths. In other words, the Z2 topo- 
logical insulators as descendants of (2n + l)-d topologi- 
cal insulators can only be defined in {{2n — 1) -I- 1) and 
((2n— 2)-|-l)-d . Since a bulk topological insulator always 
corresponds to a topologically protected gapless edge the- 
ory, the validity of a Z2 classification can also be justified 
by studying the stability of edge theories. As discussed in 
Sec lIII B1 the boundary theory of a (4 -I- l)-d topological 
insulator with second Chern number C2 contains of IC2I 



flavors of chiral (Weyl) fermions. In the simplest exam- 
ple of the lattice Dirac model (|6ip with —4c < m < —2c, 
the boundary single-particle Hamiltonian is 

-ffa(4+i) = • p 

which is topologically stable since no mass term is avail- 
able for the edge system. Under dimensional reduction 
the boundary theory of a (3-1- l)-d Z2 nontrivial insulator 
is simply given by taking — in the above equation: 

Hdi3+1) = V [axPx + CTyPy) 

which is stable in the presence of time-reversal symmetry 
since no T-invariant mass terms are available. The same 
analysis shows the stability of the edge theory of a (2-1-1)- 
d topological insulator, given by Hq^2+i) = vUzPz- When 
dimensional reduction is carried out once more, we ob- 
tain the (0+ l)-d edge of the (1 -I- l)-d insulator described 
by = 0. This just describes a Kramers's pair of 

localized states on the boundary. Since such a pair of 
mid-gap states can be easily removed by a constant en- 
ergy shift without breaking time-reversal symmetry, the 
(l-l-l)-d TRI insulator does not have a topologically non- 
trivial class. Different edge state stabilities for effective 
theories in different dimensions arc illustrated in Fig. [2T] 
Such an edge theory analysis can be easily generalized 
to higher dimensions. The boundary states of a {2n + 1)- 
d topological insulator with nontrivial Chern number are 
described by a ((2n— 1) + l)-d chiral fermion theory with 
the Hamiltonian 

2n-l 

H2n-l{Y>) = V p^T' (143) 

1=1 

in which F' are 2"^^ x 2"^^ matrices forming a represen- 
tation of the so(2n — 1) Clifford algebra. The boundary 
theory of the m-th descendant of the {2n + l)-d system 
is given by simply taking pi = Q for i — 2n — m, 2n — 
TO-|-l,...,2n— 1. The symmetry properties and stability 
of the theory can be studied by studying the properties 
of the F° matrices. Here we will display the conclusions 
of the edge state analysis, with the details presented in 
Appendix |DJ 



1. The chiral Hamiltonians (jl43p in different dimensions satisfy different discrete symmetry properties, as listed 
below: 

' C^H2u-i{v)C = -i/2«-i(-p), C*C = 1, n = 4m - 3, m e N 
rtiJ2„-i(p)r = H2n-i{-v), T*T - -I, n = 4m - 2, m G N 
C^H2n~i{p)C = -H2n-i{-p), C* C = -1, n = 4m - 1 , m E N 
fti/2n-i(p)T = i72«-i(-p), f*f = l, n = 4m, meN 

2. Only the first two descendants of H2n-i, i-e., H{p) = Pi^^ ^^'^ H{p) — X^i"^"^ Pi-^*, are topologically 
stable under the constraint of the discrete symmetry in the given dimension {C,C,T or T). Consequently, 
the Z2 topologically nontrivial insulators descending from the (2n + l)-d topological insulator only exist in 
{{2n - 1) -f l)-d and ((2n -2) + l)-d. 
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The edge state analysis confirms our insight from the bulk 
picture, that is, the Zi topological classification is only 
well defined for the first two generations of descendants 
of the (2n + l)-d topological insulator. Moreover, it also 
provides more information about the discrete symmetries 
in different dimensions. In the (6 + l)-d topological in- 
sulator, and its descendants, the correct discrete sym- 
metry is C which is similar to particle-hole symmetry C 
but satisfies C*C = —I. This is necessary since a usual 
particle- hole symmetry cannot be defined for the (5-1-1)- 

d chiral fermion 'Ylii=\'P^^ ■ Such a symmetry C can be 
called a "pseudo particle- hole symmetry". Similarly, in 
the (8 -I- l)-d topological insulator, and its descendants, 
the discrete symmetry is a "pseudo time-reversal sym- 
metry" satisfying T*T = I. In Fig. [211 the dimensions 
with "true" C or T symmetry are labelled with filled cir- 
cles, and those with C or T symmetry are labelled with 
squares. 

This paper is partly inspired by work in high-energy 
physics. The study of topological insulators with non- 
trivial boundary states is analogous to the genera- 
tion of massless fermions on higher dimensional domain 
wall S ''^r'T^i. In a high energy context these surface states 
would be subsequently used to mimic chiral fermions 
in lattice gauge theories. In our case they become the 
gapless boundary liquids that generate novel transport 
properties and characterize the topological stability of 
the state. Our picture of a (4 -I- l)-d topological insula- 
tor characterized by the Chern-Simons term has a very 
special meaning in high-energy physics when reduced to 
(3 -|- l)-d. The Chcrn-Simons term becomes a ^-term 
which is related to the so-called vacuum angle. If 6* is a 
constant then this term does not contribute to the equa- 
tions of motion. In our case the 0-term has a solitonic 
structure with a domain wall at the surface of the topo- 
logical insulator. Inside the insulator Q jumps by tt which 
still preserves CP if the original vacuum does. Thus, the 
only effect on the system is a non-zero boundary term at 
the 9 domain wall. From a high-energy perspective we 
have introduced an axionic domain wall and a topologi- 
cal insulator exists in one domain while a trivial insulator 
exists in the other. 

The dimensional reduction procedure we introduced 
here is not new to physics, and was first used in the 1920's 
in an attempt to unify gravity and electro-magnetism in 
(4-1- l)-d by Kaluza and Klein^^i^i. Basically, the dimen- 
sional reduction amounts to compactifying the "extra di- 
mension" with periodic boundary conditions, e.g. on a 
circle, and shrinking its radius to zero. The compactifi- 
cation creates a tower of modes labelled by a discrete in- 
dex, but as the circle is shrunk only one low-energy mode 
from each field remains (which can be seen by taking 
the Fourier transform of a field on a circle with periodic 
boundary conditions). These "zero- modes" become the 
propagating fields in the lower dimensional space. Car- 
rying out this procedure via our method or by compact- 
ification yields the same results. The adiabatic param- 
eter we introduced is connected with the fiux threading 



the higher dimensional circle. One could imagine com- 
pactifying using higher dimensional manifolds, such as a 
sphere or something more exotic, and threading various 
fluxes through non-trivial cycles of the compact space. 
The zero-mode structures of these manifolds are more 
complicated and we will not deal with them here, but 
perhaps other interesting theories can arise. 

Additionally, we have unearthed a ladder of topolog- 
ical insulators with gapless fermionic boundary states 
whose stability depends on the presence of a discrete 
symmetry. The discrete symmetries repeat modulo 8, 
which should be no surprise since their form is derived 
from the representation theory of real Clifford algebras 
which exhibit Bott periodicity with period 8^^. Due 
to this periodicity we can analyze the types of gapless 
fermions allowed to exist at the boundaries of topologi- 
cal insulators. Any boundary can support Dirac fermions 
and any even dimensional space-time boundary can have 
Weyl (chiral) fermions as well. For Majorana fermions 
we are restricted to boundaries with spacetime dimen- 
sions {(1 + 1), (2 -h 1), (3 -h 1), (7 + 1), (8 + 1)} mod8. 
Finally, there is a special representation, the Majorana- 
Weyl fermion, which is a real fermion with a definite 
handedness which can only exist in (1-1-1) mod 8 di- 
mensions i.e. (1 + 1), (9 + 1), (17 + 1), (25 + 1) . . 
It is this type of fermion which appears at the edge of 
(j) + ip) superconductors^SiiliZ^. However, because this 
representation is missing in (5 -I- l)-d a (6 -I- l)-d "(p-fip)- 
superconductor" which obeys the "pseudo-particle hole" 
symmetry will not have single Major ana- Weyl boundary 
states. This can also be seen by the fact that ^ —1 
which means there would have to be at the minimum 
two Majorana- Weyl fermions at the boundary due to 
a Kramers'-like theorem. Beginning with the parent 
(2n+l)— d topological insulator we see that the boundary 
theories of itself and its stable descendants are massless 
fermions in ((2n-l) + l), ((2n-2) + l) and ((2n-3) + l) 
dimensions respectively. Theories of massless fermions 
often result in field theory anomalies and there is a deep 
connection between this boundary theory ladder and the 
corresponding anomaly ladde r^^i^^ . 

In summary, we have provided a unified framework 
to describe a whole family of topological insulators in 
generic dimensions. All the topological effects are de- 
scribed by phase-space Chern-Simons theories, which ei- 
ther describe the topological insulators with nontrivial 
Chern number in odd space-time dimensions, or describe 
their lower dimensional descendants through dimensional 
reduction. Z2 topologically nontrivial insulators exist in 
((2n — 1) -I- 1) and {{2n — 2) -I- 1) space-time dimensions 
and are protected by a given discrete symmetry that is 
preserved by the parent {2n -f l)-d topological insulator. 
We found that in {2n -I- l)-d there are two types of topo- 
logical insulators, one of which is characterized by the 
Chern number C„ and the other by the Z2 invariant, as 
a descendant of ((2n -I- 2) -I- l)-d topological insulator. In 
comparison, in {{2n — 1) + l)-d there is only one type 
of topological insulator, which is characterized by a Z2 



41 



invariant, as a descendant of (2n + l)-d topological in- 
sulator. There are many tantalizing connections of our 
work with well-developed sectors of high-energy physics, 
and diving deeper into these subjects is sure to benefit 
both condensed-matter and high-energy physics. 
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VII. CONCLUSION AND DISCUSSIONS 

In conclusion we have constructed the topological field 
theory of TRI insulators. We showed that the fundamen- 
tal TRI insulator naturally exists in 4-1- 1-d, and the effec- 
tive topological field theory is the Chern-Simons theory 
in 4 -t- 1-d. We introduced the concept of dimensional re- 
duction for microscopic fcrmion models, where some spa- 
tial dimensions are compactified and the associated mo- 
mentum variables are replaced by adiabatic fields. This 
method enables us to obtain the topological field theory 
for the 3D and 2D TRI insulators from the dimensional 
reduction of the 4 -I- 1-d Chern-Simons field theory. In 
particular, we obtain the "axion" field theory for the 3D 
insulator, and many experimental consequences follow di- 
rectly from this topological field theory. The most strik- 
ing prediction is the TME effect, where an electric field 
induces a magnetic field along the same direction, with 
a universal constant of proportionality quantized in odd 
multiples of the fine structure constant. The role of the 
"axion" , or the adiabatic field, is played by a magneto- 
electric polarization, whose change is quantized when an 
adiabatic process is completed. The topological field the- 
ory for the 2D TRI insulator involves two adiabatic fields, 
and this theory directly predicts the fractional charge of 
a magnetic domain wall at the edge of the QSH insulator. 
These topological effects illustrate the predictive power of 
the topological field theory constructed in this work, and 
many more experimental consequences can be obtained 
by the proper generalization of the concepts introduced 
here. 

Our work also presents the general classification of 
topological insulators in various dimensions. The funda- 
mental topological insulators are described by the topo- 
logical Chern-Simons field theory, and the effective topo- 
logical field theory of their descendants can be obtained 
by the procedure of dimensional reduction. The descen- 
dent topological insulators are generally classified by dis- 
crete symmetries like the charge conjugation and the time 
reversal symmetries. This way, the Chern number classi- 
fication of fundamental topological insulators and the Zi 
classifications of their descendants is unified. Finally, we 
present a framework in terms of the Chern-Simons field 
theory in phase space, which gives a unified theory of 
all topological insulators and contains all experimentally 
observable topological effects. 

We would like to thank Dr. B. A. Bernevig for many 
insightful discussions on this subject and collaborations 
at the early stage of this project. We would like to thank 
H. D. Chen, S. Kachru, C. X. Liu, J. Maciejko, M. Mul- 
ligan, S. Raghu, S. Shenker and J. Zaanen for helpful 



APPENDIX A: CONVENTIONS 

Due to the special importance of dimensionality in 
this work we tried to be very consistent with our dimen- 
sion and index conventions. In addition, gauge fields of 
all types appear and we have selected a convention for 
the electro-magnetic vector potential and the adiabatic 
(Berry's phase) connection. The conventions we used in 
this paper are summarized in Table |TT] and also explained 
below. 

For spacetime conventions we have chosen the form 
(n + l)-d where n can be 0. For spatial dimensions only 
we use nD with a capital D. 

For indices, space-time and frequency-momentum in- 
dices share the same convention. Greek indices from the 
middle of the alphabet such as /i, i^, p, cr, t run from to 
the spatial dimension in the current context. Examples 
being 0,1, 0,1,2 etc. Additionally, when indexing mo- 
mentum space objects is frequency and 1,2,3... are 
the momentum components in the 1,2,3,... directions. 
Latin indices from the middle of the alphabet such as 
i, j,k,£ are purely spatial indices and run from 1 onward 
to the spatial dimension of the current context. For mo- 
mentum space objects they index the spatial momenta 
e.g kx, ky, kz, kw . . . . We always use the Einstein sum- 
mation convention unless stated otherwise. Since we are 
considering flat space we make no distinction between 
raised and lowered indices. 

In some special contexts we will need two more sets of 
indices. Several of the Hamiltonian models we use can 
be written in terms of the 2 x 2 or 4 x 4 Dirac matrices. 
In these cases where you see vectors indexed by low- 
ercase Latin letters from the beginning of the alphabet 
they run from 1, 2, 3 and 0, 1, 2, 3, 4 respectively. Finally, 
for cases where the indices don't just run over coordi- 
nate or momentum space separately, but instead cover 
"phase-space" coordinates, we use capital Latin letters 
from the beginning of the alphabet such as A, B,C, D, E. 
These run over the phase space variables in the cur- 
rent context. Some examples being = {t, x, y, k^, ky), 
= (i, x,kx), or A running over {kx, ky, 6, 0). 

For orbital, band, or state labels we use Greek letters 
from the beginning of the alphabet such as a, /?, 7. 

For the electro- magnetic U{1) gauge field we use 
where A is capitalized. For the Berry's phase gauge field 
we use a* with a lower case a. Note the different index la- 
bels. The electromagnetic gauge field has a 0-component 
while the Berry's phase has no "frequency" -component. 
For the curvatures we use F^i, and /.y respectively. 
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symbol 


explanation 


(71+ l)-d 


space-time dimension (n + 1) 


nD 


spatial dimension n 


^i,i/,p,o-, ... 


space-time or frequency-momentum indices 0, 1, ..,n 


^1 i) ^1 ^1 •■• 


spatial or momentum indices 1, 2, .., n 


A,B,C,D,... 


phase space indices 0, 1, ...n, n -I- 1, ..2n -I- 1 


a, 6, c, d, ... 


indices for anticommuting F" matrices and corresponding coefficients da 


Q,/3,7,(5... 


energy band indices 




gauge vector potential and gauge curvature of external electro-magnetic field 




(generally non-Abelian) Berry phase gauge vector potential and gauge curvature in momentum space 




electro-magnetic or Berry phase gauge vector potential defined in phase space, respectively 



TABLE II: Table of conventions used in the paper. 



APPENDIX B: DERIVATION OF EQ. (HH) will demonstrate these three steps separately. 



In this appendix we will prove the conclusion ([51]) . As 
is briefly sketched in Sec. IIII Al the demonstration con- 
sists of three steps: (1) Topological invariance of Eq. 
(155)1 : (2) Any Hamiltonian h{k) is adiabatically con- 
nected to an /io(k) with the form of Eq. ((55|) : (3) For 
such an /io(k) the non-linear correlation function (j53p is 
equal to the second Chern number. In the following we 



1. Topological invariance of Eq. (|53p 

To prove the topological invariance of Eq. (|53p . we 
just need to prove any infinitesimal deformation of the 
Green's function G(k, to) leads to a vanishing variation 
of C2. Under a variation of G'(k, lo) we have 



5 {Gdf,G-^) = SGdf,G-^ + Gd^ (SG-^) = SGd^G-^ - Gd^, {G'^SGG-^) 
= -G{d^,G-^)5GG-^ -dt,{5G)G-^. 

Thus the variation of C2 is 

SG2 = -^e^-P-^ [ ^-^tT[S {Gd^G-') {Gd^G-') {GdpG-') {Gd,G-') {GdrG-')] 
15 J (27r) 

= ^e^""''" / ^^tr [{Gd^G-^SGG-^) {Gd,G-^) (GdpG-^) (Gd^G-^) {GdrG-')] 

f d^^^ [{d,5GG-') {Gd.G-') {GdpG-') {Gd.G-') {GdrG-')] 
l-^ J (27r) 

^ ^^^upar f ^!^tr [d,, (G-^SG) {d,G-^G) (dpG-^G) {d^G'^G) {drG-^G)] 
15 J (27r) 

^ ^^^.upar f ^!^a^tr [(G-MG) (d^G-^G) (dpG-^G) (d^G-^G) (drG-^G)] 
15 J (27r) 

_^^,.p.r f ^^^-1^^^ {d.G-^G) {dpG-^G) (9.G-1G) (a^G-^G)] 

1'3 J (27r) 

+tr [{G-^5G) {d^-^G) (dpG-^G) (S^G^^G) (a^G^^G)] 

+tr [{G-^SG) {d^G-^G) (dpG-^G) {d^G-^G) (a^G^^G)] 

+tr [(G^^G) (d^G-^G) (dpG-^G) (S^G^^G) dp (a^G^^G)] } 

= /■ ^^a^tr[(G-i<5G) (a^G-^G) (9,G-1G) (a^G-^G) (a^G-^G)] 

15 J (27r) 



(Bl) 
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Thus Eq. ((53)) is topologically invariant. The topologi- 
cal invariance of the second Chern number defined in Eq. 
([54| is a well-known mathematical fact. Such a topologi- 
cal invariance is quite helpful for showing the equivalence 
between the second Chern number and Eq. ([53]) . 



we obtain an adiabatic interpolation between /i(k, 0) = 
h{k) and h(k, 1) = U{k)Do{k, l)C/t(k). Since the matrix 
U{k) can be written in the eigenstates of h{k) as U{k) = 
(|l,k),|2,k),...,|7V,k)), wehave 



2. Adiabatic deformation of arbitrary /i(k) to ho(k) 

An adiabatic deformation h(k, t), t G [0, 1] can be writ- 
ten down, which connects an arbitrary gapped Hamil- 
tonian h{k) to a "maximally degenerate Hamiltonian" 
ho{k) in the form of Eq. ([55]) . Any single particle Hamil- 
tonian h{k) can be diagonalized as 

h{k) = U{k)D{k)U^{k) 

with U{k) unitary and ho{k) = 
diag [ei(k), e2(k), eAr(k)] the diagonal matrix of 
energy eigenvalues. Without loss of generality, the 
chemical potential can be defined to be zero, and the 
eigenvalues can be arranged in ascending order. For an 
insulator with M bands filled, one has 



M N 
a=l I3=M+1 

= eaPam + esPEi^k) (B5) 



In summary, we have proven that each gapped Hamil- 
tonian h{k) can be adiabatically connected to a Hamil- 
tonian with the form of Eq. (|55p . 

3. Calculation of correlation function (|53[) for hoik) 



For Hamiltonian of the form (IB5P the Green's function 
is written in the simple form: 



ei(k) < e2(k) < ... <eM(k) <0 
< eAf+i(k) < ... < eAr(k). 



For t e [0, 1], define 
E4k,t) -- 



ea(k)(l - <) + eci, l<a<M 
ea(k)(l -t) + est, M <a<N 



(B2) 



(B3) 



^ Pg jk) ^ PE{k) 

uj + iS — eG u! + iS — eE 



On the other hand, we have 



(B6) 



and Do{k,t) = dia.g[Ei{k, t), E2ik,t), En {k,t)], then 
we have 



i^o(k,0) =i^(k), Do{k,l) 



xM 



MxN-M 



As long as ec < < , -Do (k, t) remains gapped for 
t e [0, 1]. Thus by defining 



h{k,t) = U{k)Da{k,0)U\k) 



dG'^{k,uj) 
dio 

dG-^{k,uj) 
dki 



= 1 



dPG{k) dPE{k) 

-CG — ^E- 



(e_E - tG 



dki dki 



dki 



(B4) where i = 1, 2, 3, 4. Thus Eq. ([52]) can be written 



(B7) 



3 J 2n 



(2^^)^ n,m,s,t=1.2 {^^ + iS - {uJ + iS - £,„) {uj + i6 ~ Cs) {uj + iS - Ei) 



Tr 



p dPap dPg p dPc p dPG 
^ri Ok, dkj a/cfc dke 



(ee - Eg) 



(B8) 



in which ei_2 — eg^e and Pi,2(k) = Pg,e{^), respectively. From the identity Pe + Pg = ^ and P| = Pe, Pq = Pg, 
we obtain 



P, 



DPg 
dki 



dh 



dki 



^^d_PE ^ d_PGp^ 

dki dki 



(B9) 
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Consequently, PcdiPcPG = PEdiPaPE = 0, so that the trace in Eq. (jB8[) can be nonzero only when n ^ m, m ^ 
s, s t, t n. In other words, only two terms are left out of the 16 terms summed over in Eq. (|B8p : 



C2 



d kduj 

7^ 



Tr 



Pg^Pe^Pg^Pe^ 



Tr 



P 9Pg n OPg p dPg p aPc 

Pe^Pg^Pe^Pg^ 



{uj + i5 — ec)^ {ui + i5 — e^)^ 



(tj + — ec) (cj + i(5 — 



(££; - Eg) 



Carrying out the integral over lo and using identities (|B9|) again, we obtain 



C2 = / d^fce'^'^^T^ 



dPcdPc^ OPgOPg 

^E-^, — ^E 



dki dki 



dkk dki 



(BIO) 



Now we will show this is just the second Chern number. 
The Berry phase gauge field is defined by 



d 



ar(k) = -^(a,k| — |/3,k) 

in which a,f3 = 1,2, ..,M are the occupied bands. The 
U (M) gauge curvature is given by 

. f d{a,k\d\P,k) 
= -'[-^. dk--^'^'^ 

. I d (a,k 



.f d{a,k\ r. n.. 9\P,k) 
I '^^ '~dk ^ -^^ 



in which (z means the term with i,j exchanged. In 
operator form, we have 



dki 



Thus Eq. (|B10[) can be written in terms of the Berry 
phase curvature as-'*^'''- 

C2 = ^ 1 d^fce'^'^^Tr [f,,U . (Bll) 



In conclusion, the equivalence of the non-linear corre- 
lation function (|53|) and the second Chern number ((54|) 
has been proven for the specific models of the form Eq. 
([55]) . Due to the topological invariance of both Eq. ((53|) 
and the second Chern number, such a relation holds for 
any band insulator in (4 + l)-d. 

The procedure in this section can be easily generalized 
to any odd space-time dimensions. In {2n + l)-d space- 
time, a n-th Chern number C„ is defined in the BZ of a 
band-insulator, which appears as a topological response 
coefficient to an external gauge field. Equivalently, the 
coefficient of the n-th Chern-Simons term in the effective 
action of an external gauge field obtained from integrat- 
ing out the fermions^. 

APPENDIX C: THE WINDING NUMBER IN 
THE NON-LINEAR RESPONSE OF 
DIRAC-TYPE MODELS 



Starting from the symmetric form, in terms of general 
Green's functions, of Eq. ([55]) we want to calculate C2 
for the generalized lattice Dirac model H{k) — (ia(fc)r". 
The Green's function for this model can be simplified to^ 



G(k,u;) 



Inserting this into Eq. 
tion 



uj + d''{k)r'' 

{UJ-^ - d^{k)da{k)y 



(CI) 



3]) yields the following calcula- 



C2 = -yey-M j "Isty (c.2 - |rfP)5 ^^ + drnVm^ + dnTn{dd°T°){uj + d^T^id^d^Vi) 
X (CJ + drV){dkd''T'){uj + dtr*)(9;rf"r")] 

= 'Y'-'"' J (271-7 (^^ -1^1^)5 [^'^ + d„,Tn{co + d„r")r°(c. + d^rnviuj + drVnT^ico + d,r*)r«] 

I 



where ni,n,o,p,q,r,s,t,u = 0,1,2,3,4. Among all the are those of 5,7, and 9 T matrices''^ and the complete 
terms in the bracket [], the only ones with non-zero traces trace simplifies nicely to — 4e(o5s«rf*(w^ — |fip)^, which 
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simplifies C2 as 



C2 = —eijkl 



-4e 



(27r)5 (^2 _ |rf|2)3 



TT' / d^/c d^did^djd'idkd'did'^ 



\d\' 



, d^ketoqsud'dJ°dyd'id,d'd^d^. (C2) 

Thus we have proved that the winding number given in 
the expression (j64p is equal to the second Chern number 
defined in Eq. for generic (4 + l)-d Dirac models. 



APPENDIX D: STABILITY OF EDGE THEORIES 
IN GENERIC DIMENSIONS 

In this appendix, we will determine the existence or ab- 
sence of Z2 topological insulators in generic dimensions 
by studying the stability of boundary theories. First 
of all, the boundary theory of a topological insulator in 
(2n + l)-d with nontrivial n-th Chern number C„ is |C„| 
copies of chiral fermions: 



C„| 2n-l 



(Dl) 



The 2n — 1 r° matrices are 2"~^ dimensional and form 
an so(2n — 1) Clifford algebra. We will study the lattice 
Dirac model, since other systems with the same Chern 
number can be obtained by an adiabatic deformation 
from the lattice Dirac model and the topological stability 
of edge states does not depend on the adiabatic deforma- 
tion, as will be discussed below. For simplicity, in the 
following we will focus on the case C„ = 1, in which 
case the single particle Hamiltonian can be written as 
h{p) = wpaF". We have the following theorem about the 
symmetry of the Hamiltonian: 

• Theorem I. In a 2""^ dimensional representa- 
tion of an so(2n — 1) Clifford algebra generated by 
F", a = 1, 2, ..2n — 1, a unitary matrix M(„) can be 
defined such that 



(D2) 
(D3) 



Consequently the chiral Hamiltonian h{p) = vpaT"^ 
satisfies 



Aft /,(p)M(„) - (-l)"/i*(-p). 



(D4) 



In Eq. (jD3[) [n/2] means the maximal integer that does 
not exceed n/2. 

Theorem I can be proved by induction. First, the F" 
matrices for n = 2 are the Pauli matrices F" = CTq , a = 



1, 2, 3, and M(2) = i(T2 satisfies the theorem. Suppose the 
theorem is true for case n with the matrix M(„) and the 
2n — I Gamma matrices F?„ , s , then the F?„ , , x matri- 

{271—1) ' {2n+l} 

ces for the so(2n -I- 1) Clifford algebra can be generated 
by 



^ (2n+l) 



r°2„-i)®^y, a = l,2,..2n-l 

I«)T^, a^2n ■ (D5) 

I(g)Tz, a = 271 + 1 



It is straightforward to check the anticommutation rela- 
tions. Defining 



Ml 



M, 



(n+l) 



(n) 



Ml 



in) 



>iTy, n odd 
g) I, n even 



(D6) 



we find M(„_|_i) satisfies the Theorem I. 

From Eqs. (|D3[) and (|D4[) one finds different prop- 
erties of A/(„-| in different dimensions, (i) For n = 
4A: — 3, k E Z, M(„-) reverses the sign of energy, and sat- 
isfies M^*^^M(„) = I, which can be identified as particle- 
hole symmetry; (ii) for n = 4k — 2, fc € Z, M(„) pre- 
serves the energy and satisfies M^*^^M(„) = —I, which 
can be identified as the time-reversal symmetry; (iii) for 
n — Ak — 1, k G Z, M(„-) reverses the sign of energy but 
satisfies M*j^^M(„) = —I, which we call "psuedo" particle- 
hole symmetry and denote by C; (iv) for n = 4fc, k E Z, 
Af(„) preserves the energy but satisfies Afj*^^M(„) — I, 
which behaves like the time-reversal symmetry of an 
integer-spin particle, and we call "pseudo" time-reversal 
symmetry T. 

In the following, we will denote the symmetries C, C, 
T and T, defined by Eq. (|D4p , by Af-symmetry. Now we 
study what other terms can be added in the Hamiltonian 
without breaking the M symmetry. Given F° for the 
so(2n — 1) case, all the 2'^~^ x 2"^^ Hcrmitian matrices 
can be expanded in the basis 



IS" 



^m(m — 1) /2pai pa2 pdn 



1,2, 



As expected, the total number of matrices forming the 
basis is 1 + X;^=\ {Tn-i) = 2^""^- The Af(„) transfor- 
mation property of r^i^^- """ can be determined by that 
of F° as 



(n) (n) 



If we have a constant term maja2..a,„r°^°^"''™ in the 
Hamiltonian without breaking the M symmetry, the fol- 
lowing condition must be satisfied: 



Eqs. (|D7l) and ((D8|l thus require 

^ -^yn(2n+m — 3) /2 ^ -^^n 



(D8) 



(D9) 



so that (to— l)(2n— 2-f to)/2 must be odd. This condition 
can be satisfied by several possibilities: (1) when m is 
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odd, (to— l)/2 must be odd, which means m = 4fc— 1, k G 
N. (2) When to is even, n — 1 + to/2 must be odd, which 
means to = 4fc if n is even, or to = 4A: — 2 if n is odd. 
In summary, the terms available in the Hamiltonian are 
given by 

{4/c — 1 or 4k. n even , 
4fc - 1 or 4fc - 2, n odd 

in which fc € N and to, is also bounded by 1 < to < n. 
For example, the first non-trivial case is n = 2 in which 
r°,a — 1,2,3 are Pauli matrices. Since n is even, to is 
required to be 4fc — 1 or 4fc, in which the lowest value is 
TO = 3. Consequently there is no mass term available. 
When n — "i {2n —1 = 5) the F", a = 1,2, .., 5 matrices 
are the usual Dirac matrices, and the only value of to 
satisfying Eq. (jPlOp is to = 2. In other words, the terms 
pab jpop') fj^Q jjot break the corresponding symmetry — 
the pseudo particle-hole symmetry C. 

Though there are all these constant terms avail- 
able, the perturbed Hamiltonian h{k) = AaF" + 
X^m '^aia2..a„r"i"2-"™ remains gapless because each 
paia2..a,„ commutcs with some F°. If to is even, it com- 
mutes with F*", 6 7^ fls V.S = 1, .., TO. If TO is odd, it com- 
mutes with F°=, Vs = 1, ..; TO. For the Hamiltonian with 
only one constant term maia2...a^^"'^°'^"'°'"^ , one can take 
Pi = for all i except for i — a, where a is chosen such 
that F° commutes with r°i°^- "™. Due to the commu- 
tativity, the Hamiltonian can be diagonalized along the 
a-th axis, with the eigenvalues Pa ± to. Consequently, 
we know the Hamiltonian is gapless. In other words, we 
have shown that no perturbation with the 2"~^ band the- 
ory can open a gap for the 2n dimensional chiral fcrmion, 
which agrees with the topological stability of the (2n-f-l)- 
d bulk system characterized by the n-th Chern number. 

Starting from the (2n -I- l)-d topological insulator, di- 
mensional reduction procedures can be carried out to 



obtain a ((2n — 1) -f l)-d topological insulator. Cor- 
respondingly, the edge theory of the ((2n — 2) -I- l)-d 
topological insulator is given by the dimensional reduc- 
tion of the ((2n — 1) -f l)-d chiral fermion. If a non- 
trivial topological insulator can be defined, it has the 
boundary theory /i(p) = X^a^T^PaT'^- Compared to the 
((2n — 1) -I- l)-d chiral fermion, one F"^ matrix is absent 
in the theory. In the same way, the boundary theory 
of lower dimensional descendants can be obtained by re- 
moving more momenta and F° matrices from the chiral 
fcrmion Hamiltonian. Obviously, if too few F" matri- 
ces are left, a mass term will be available, which anti- 
commutes with all the rest of the F" matrices and thus 
can make gap the whole edge spectrum. The upper crit- 
ical dimension where the edge states become unstable is 
determined by the maximal number of F° matrices that 
anti-commute with some rn^a-.a™ ^ -p^j. ^ even, the max- 
imal number of F^'s that anticommute with p^iaa-.a^ 
TO, while for to odd, the maximal number is 2r7, — 1 — to. 
On the other hand, the available values of to are de- 
fined by the constraint Eq. (|D10p . By studying the cases 
n — 4k — 3,4k ~2, 4k— 1, 4k, k G N separately, we obtain 
that the space-time dimension in which the chiral theory 
becomes unstable is given by d = 2n — 3. For example, in 
the case n = 3 discussed earlier, the mass terms F"'' are 
permitted by the pseudo particle-hole symmetry. When 
the dimension is reduced from (5 + l)-d to (2 + l)-d, only 
two F" matrices are used in the fc-linear terms, so that 
some F"'' can be found which anti-commutes with the 
gapless Hamiltonian and thus can make the surface the- 
ory gapped. Consequently, topologically stable boundary 
theories as descendants of ((2n— 1) -I- l)-d chiral fermions 
can only exist in ((2n — 2) + 1) and {{2n — 3) -|- 1) dimen- 
sions. Correspondingly, the Z2 topological insulators as 
descendants of the {2n+ l)-d topological insulator (with 
n-th Chern number) can only be defined in {(2n — 1) + 1) 
and {{2n — 2) + 1) dimensions. 
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